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HAUPTAUFSÄTZE 


Über eine Schwingungserscheinung in Flüssigkeiten 


mit stabiler Dichteschichtung. 
Von H. Görtler in Göttingen. 
(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) . 


Die Theorie der erzwungenen kleinen Schwingungen in einer der Schwerkraft 
unterworfenen reibungslosen Flüssigkeit mit stabiler Dichteschichtung führt 
zur Voraussage einer neuartigen Schwingungserscheinung: Die Differential- 
gleichung zweiter Ordnung für die Verschiebungsamplituden zeitlich 
harmonischer Schwingungen als Funktionen der Ortskoordinaten ist für 
Frequenzen unterhalb eines von der jeweiligen Dichteschichtung abhängigen 
kritischen Werts vom hyperbolischen Typ. Demgemäf treten Regularitäts- 
defekte dieser Größen an einem Orte zugleich in allen Punkten der durch 
diesen Ort gehenden reellen Charakteristiken der Differentialgleichung auf. 
In einfachen Experimenten läßt sich diese Erscheinung beobachten. Die 
theoretischen Gesetzmäfiigkeiten des Charakteristikenverlaufs im einzelnen 
erfahren dabei ihre Bestätigung. 


8 1. Einleitung. 


Kleine Schwingungen in reibungslosen, der Schwerkraft unterworfenen Flüssigkeiten 
oder Gasen mit in der Ruhelage horizontaler und stabiler Dichteschichtung sind mehrfach 
theoretisch untersucht worden‘). Man interessierte sich vorwiegend für das Problem der 
Wellenausbreitung in solchen Medien mit teilweise festen oder freien Berandungen. Auf 
diese Untersuchungen soll hier nieht näher eingegangen werden. Die vorliegende Betrachtung 
gründet sich auf einen bemerkenswerten mathematischen Sachverhalt bei den Bewegungs- 
gleichungen zeitlich harmonischer Schwingungen, dem bisher merkwürdigerweise nicht die 
verdiente Beachtung geschenkt worden ist. Dies überrascht um so mehr, als bereits durch 
sehr einfache und naheliegende Schlußfolgerungen aus diesem Sachverhalt theoretische Vor- 


1) Genannt seien: Lord Rayleigh: Investigation on the character of the equilibrium of an ineompressible 
heavy liquid of variable density. Proc. Lond. Math. Soc. Bd. 14 (1883), S. 170; W. Burnside: On the small wave-motions 
of a heterogeneous fluid under gravity. Proc. Lond. Math. Soc. Bd. 20 (1889), S. 392; A. E. H. Love: Wave motion in 
a heterogeneous heavy liquid. Proc. Lond. Math. Soc. Bd. 22 (1891), S. 307; H. Lamb: On the theory of waves 
propagated vertically in the atmosphere. Proc. Lond. Math. Soe (:) Bd. 7 (1909), S. 122; On atmospherie oseillations. 
Proc. Roy. Soc. Bd. 84 (1911), S. 551; S. Sano: Some problems on waves propagating in air at an uniform temperature. 
Bull. Central. Meteorol. Observatory, Japan, Bd. 2 (1913), S. 157. Siehe ferner: V. Bjerknes, J. Bjerknes, 
H. Solberg und T. Bergeron: Physikalische Hydrodynamik. Berlin 1923 (Achtes Kapitel) und H. Lamb: Lehrbuch 
der Hydrodynamik. Deutsch von E. Helly. 2. Aufl., Leipzig und Berlin 1931, S. 423 bis 425. 

Den Hinweis auf die Arbeit von S. Sano verdanke ich Herrn H. Ertel. 
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aussagen über eine -Schwingungserscheinung hergeleitet werden können, für welche m. W. 
experimentelle Belege noch nicht vorliegen. In einfachen Versuchen konnte, worüber im 
folgenden ebenfalls berichtet werden soll, diese Erscheinung nachgewiesen werden, und es 
wurden dabei die Voraussagen der Theorie auch quantitativ bestätigt. 

Fragen über die Schwingungsvorgänge im einzelnen sollen in dem vorliegenden Bericht 
weder von theoretischer noch von experimenteller Seite erörtert werden. Das Studium der 
Bewegungsformen soll vielmehr einer späteren ausführlichen Untersuchung vorbehalten bleiben. 


$ 2. Problemstellung. 

Wir beziehen die folgenden Rechnungen auf ein cartesisches Rechtssystem x, y, 2, dessen 
positive z-Achse vertikal nach oben weist. In bezug auf dieses System sei die Dichte o, 
einer reibungslosen, volumbeständigen, der Schwerkraft unterworfeneu und in der Gleich- 
gewichtslage befindlichen Flüssigkeit eine reine Funktion von z°?): 

Von der durch diese Funktion festgelegten Dichteschichtung wir ausdrücklich 

voraussetzen, daß sie nirgends eine instabile Gleichgewichtsanordnung darstellt. Dazu ist 

notwendig und hinreichend, daß o,(z) im Definitionsbereich mit wachsendem 2 nirgends 
zunimmt’), d. h., daß für alle 2 gilt: 

Durch äußere Einwirkung, etwa durch die Bewegung eines in die Flüssigkeit einge- 
tauchten oszillierenden Körpers, möge die Flüssigkeit in zeitlich harmonische Schwingungen 
versetzt werden. Von den dadurch verursachten synchronen periodischen Teilchenver- 
schiebungen £, n, £ in x- bzw. y-, z-Richtung sowie von den zugehörigen Druck- und Dichte- 
änderungen und auch von den räumlichen und zeitlichen Ableitungen dieser Größen setzen 
wir voraus, daß sie hinreichend klein sind, um eine Linearisierung der hydrodynamischen 
Bewegungsgleichungen bezüglich dieser Größen zu gestatten. 

Wenn ein Flüssigkeitsteilchen aus der Gleichgewichtslage in Höhe 2—{£ in die Höhe z 
verlagert wird, hat es, da es voraussetzungsgemäß volumbeständig sein soll, gegenüber der 
dort nach (2.1) definierten Ruhedichte o,(2) die Dichte o(2)=o,(2—{). Für diese können 
wir auf Grund der vorausgesetzten Kleinheit der Verschiebung £ mit hinreichender Näherung 


setzen. 
Die bezüglich der Störungsgrößen linearisierten Bewegungsgleichungen lauten demnach: 


Hierin bezeichnen p den statischen Druck und g die Erdbeschleunigung, während die Indizes 
x,y,2,t partielle Ableitungen nach den Ortsköordinaten bzw. nach der Zeit anzeigen. 


Als vierte Bestimmungsgleichung tritt die Kontinuitätsgleichung hinzu: 


Es bezeichne nunmehr die reelle Größe £ die Kreisfrequenz der kleinen Schwingungen. 
Wir können uns darauf beschränken, £ als nicht negativ anzunehmen. Dann setzen wir an: 


(2.6) 


P(x,y,2) 
20 


(p, = statischer Druck im Gleichgewichtszustand in der Höhe 2=,). 


$& 3. Theoretische Voraussagen für ein einfaches Beispiel. 


Nunmehr soll mit Rücksicht auf die in $4 zu schildernden Experimente zunächst eine 
einfache Schichtung ins Auge gefaßt werden. (Die entsprechende Behandlung des allgemeinen 
Falls folgt erst in $ 5.) Es sei nämlich 


2) Grundsätzlich setzen wir voraus, daß 00 sowie die im weiteren eingeführten Funktionen und ihre sämtlichen 
vorkommenden totalen bzw. partiellen Ableitungen im jeweils zugrunde liegenden Bereich stetig, oder, wo dies ausdrück- 
lich gesagt wird, stückweise stetig sind. 

Lord Rayleigh:a.a.O. 
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(h eine feste Länge, o,=0,(0)). Ferner beschränken wir uns aus demselben Grunde auf 
zweidimensionale Bewegungen, indem wir annehmen, daß Y=0 ist und daß demzufolge 
X, Z und P von y nicht abhängen. Dann führt der Ansatz (2.6) in (2.4) und (2.5) eingesetzt 
auf die drei Gleichungen: 


P®X—P,=0 


X,+2=V0 


Die Koeffizienten sind in diesem Falle durchweg Konstante. (Sie sind dies auch nur in 
dem gewählten Sonderfall des Schichtungsgesetzes (3.1).) Bezeichnet L die Determinante 
des homogenen Gleichungssystems (3.2) für die Unbekannten X, Z und P im Sinne der 
Operatorenrechnung: 


L=|0 , 
, 


so ergibt der Eliminationsprozeß aus den Gleichungen (3.2) 


Der lineare Differentialoperator zweiter Ordnung L, hat nun die bemerkenswerte Eigen- 
schaft, für 


> elliptisch » 
< hyperbolisch 


zu sein. Dies gilt jeweils einheitlich für alle Stellen des flüssigkeitserfüllten Raumes. Die 
Charakteristiken von L sind, da dieser Differentialausdruck konstante Koeffizienten besitzt, 
gerade Linien, und zwar bilden sie die zwei Geradenscharen: 


. . . . 


Sie stellen im hyperbolischen Fall (#?< g/h) zwei reelle Parallelbüschel dar, deren Neigungs- 
winkel @ gegen die «-Achse der Beziehung 


genügt. Hierin ist £,; die für jede Schichtung (3.1) in der ganzen Flüssigkeit konstante 
„kritische Kreisfrequenz* 


des parabolischen Grenzfalls. Es ist dies die Frequenz eines mathematischen Pendels der 
Länge hh wo h nach (3.1) den Höhenabstand zweier im Dichteverhältnis e:1 stehenden 
Schichten darstellt. 

Faßt man nacheinander verschiedene feste, in ihrer Aufeinanderfolge von Null an 
wachsende Frequenzen $< ß; ins Auge, so verlaufen die zugehörigen reellen charakteristischen 
Geraden immer steiler, um schließlich für die kritische Frequenz $=ß, in den Vertikalen 
x=const. paarweise zur Deckung zu kommen. Für 8->f, bilden die charakteristischen 
Geraden zwei imaginäre Geradenbüschel. 

Nun sind die Charakteristiken einer linearen hyperbolischen Differentialgleichung da- 
durch ausgezeichnet, daß längs dieser reellen Linien Regularitätsdefekte der Lösungen 
auftreten können*). In Anwendung auf das vorliegende Problem kann es somit in der ge- 


%) Mathematisch sind die Charakteristiken zunächst als Ort möglicher Unstetiekeiten in den zweiten (oder 
höheren) Ableitungen stetiger und stetig differenzierbarer Lösungsfunktionen der Differentialgleichung ausgezeichnet. 
Über die Frage, unter welchen zusätzlichen Bedingungen die Charakteristiken auch als Ort möglicher Unstetigkeiten 
der Lösungsfunktion selbst oder von deren ersten Ableitungen mathematisch ausgezeichnet sind — wenn keine weiteren 
Forderungen gestellt sind, können solche Unstetigkeiten im Gegensatz zu denen höherer Abieitungen ja längs beliebiger 
Kurven auftreten — findet man gewissen Aufschluß in: R. Courant und D. Hilbert: Methoden der mathematischen 
Physik, Bd. 2. Berliu 1937, S. 859 ff, Es genügt, daß die unstetigen Lösungen Grenzfälle von Lösungen mit stetigem 
Verhalten sind. Es leuchtet daher vom physikalischen Standpunkt ein, daß die Charakteristiken auch in dieser 
Hinsicht eine ausgezeichnete Rolle spielen ‚werden. 
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scehichteten Flüssigkeit gewisse Linien geben, längs denen die physikalischen Größen, welche 
der Differentialgleichung (3.4) genügen, und aus ihnen abgeleitete Größen Unstetigkeiten 
besitzen. Es muß hier jedoch betont werden, daß unsere Theorie insofern eine Idealisierung 
darstellt, als die stets vorhandene Zähigkeit der Flüssigkeit vernachlässigt worden ist. Die 
als möglich erkannten mathematischen Unstetigkeiten der Theorie können ihre physikalische 
Verwirklichung nur ünter einer entsprechenden Modifizierung finden. 

Es erscheint angebracht, vor jeder weitergelienden theoretischen Erörterung, insbesondere 
vor expliziter Lösung einzelner Randwertprobleme, zunächst einmal das Experiment zu be- 
fragen. Wir werden untersuchen, was sich ergibt, wenn die erzwungenen zeitlich harmonischen 
Schwingungen durch einen in die geschichtete Flüssigkeit eingetauchten oszillierenden Körper 
verursacht werden. Wir können nach der Theorie jedenfalls erwarten, daß für Frequenzen 
ß< ßr Schwingungsvorgänge stattfinden werden, welche von dem bei homogenen Flüssigkeiten 
gewohnten Bild völlig abweichen, und zwar wird die Bewegung in einer Form erfolgen, daß 
dabei von der Körpernachbarschaft als Quelle lokaler Störungen (Ausweichströmung) aus- 
gehende Linien das Bewegungsbild aufteilen. In besonders einfacher Weise kann man örtliche 
Unregelmäßigkeiten im Verlauf des Dichtegradienten als Schlieren in Schattenbildern 
sichtbar machen. Treten solche Schlieren auf, so können wir an ihnen ohne irgendwelchen 
größeren experimentellen Aufwand die quantitative Gültigkeit der Gesetze (I) und (II) nach- 
prüfen. Der Bericht über die experimentelle Bestätigung dieser Aussagen bildet den Inhalt 
des folgenden $4. Was jedoch im einzelnen in der Flüssigkeit vor sich geht, soll, wie bereits 
einleitend bemerkt, erst in einer späteren ausführlichen Arbeit theoretisch und experimentell 
untersucht werden. 

Es mögen hier noch einige ergänzende Bemerkungen Platz finden, die auf das Auftreten 
der besonderen Frequenz f;, der kritischen Frequenz unserer erzwungenen Schwingungen, 
bei anderen Fragestellungen und auf eine Verallgemeinerung Bezug nehmen. 


Die Frequenz x =yjg/h tritt bei V. Bjerknes und Mitarbeitern a. a. O. und auch 
schon bei Lord Rayleigh a. a. O. auf mit der Bedeutung der Frequenz möglicher „Trans- 
versalwellen“ (horizontal sich fortpflanzende harmonische ebene Wellen, bei denen die einzelnen 
Massenteilchen vertikal schwingen) in einem inkompressiblen Medium mit der Massenverteilung 
einer isothermen Atmosphäre. Ein solches Medium besitzt gerade eine dem Gesetz (3.1) 
genügende Schichtung, wobei hier h die Höhe der gleichförmigen Atmosphäre, also 8km, ist. 
Als kritische Periode T;=2n/ß, eines solchen Mediums in unserem Sinne ergibt sich 
T=3 Minuten. Alle langsameren Perioden leben demnach die Möglichkeit des Auftretens 
reeller Unstetigkeitslinien. 

In der bereits zitierten Arbeit von S. Sano, in welcher die Wellenausbreitung in einer 
isotherm geschichteten, adiabatisch kompressiblen Atmosphäre untersucht wird, wird die 
fundamentale Rolle einer Fallunterscheidung analog (3.5) zwar erkannt, jedoch wird von der 
von uns oben gefolgerten Möglichkeit kein Gebrauch gemacht. Ohne den Charakteristiken- 
begriff heranzuziehen, gibt Sano zwar partikuläre Integrale seines Schwingungsproblems an, 
welche auf bestimmten Flächen (Ebenen, Kegel) — es sind gerade Charakteristiken- 
- flächen! — unendlich werden. Ich entnehme der Sanoschen Arbeit, daß in diesem Falle 
bei qualitativ gleicher Erscheinung — denn auch hier handelt es sich um eine lineare 
Differentialgleichung mit konstanten Koeffizienten — unsere kritische Frequenz durch 
Pr=YVr—I1gje gegeben ist («=c,lc„, e=yx RT Schallgeschwindigkeit), also, anders ge- 


schrieben, 8; = Yg/ha, wo = h, die Höhe der adiabatischen Atmosphäre (h, Höhe 


der gleichförmigen Atmosphäre) ist, so daß sich also hier mit 7=273° und «=1,41, also 
ha=28km, eine kritische Periode Minuten ergibt. (Die Periode tritt 
natürlich auch wieder als „Iransversalperiode“ bei V. Bjerknes und Mitarbeitern (s. oben) auf.) 


Bemerkenswert ist, daß uns diese Periode auch in einem anderen Zusammenhang 
begegnet, nämlich in der Theorie der Leewellen hinter einem Gebirgskamm in einer isotherm 
geschichteten Atmosphäre nach G. Lyra°). Die Wellenlänge dieger Leewellen ist hinreichend 
weit hinter dem Berg mit guter Näherung durch A=2rüc/yx» —1g (&=ungestörte An- 
strömgeschwindigkeit < c) gegeben. Ein Massenteilchen mit der mittleren Geschwindigkeit & 
führt somit um eine mittlere horizontale Bahn Schwingungen mit der Kreisfrequenz 


Pk —1gje aus. 


5) G. Lyra: Über den Einfluß von Bodenerhebungen auf die Strömung einer stabil geschichteten Atmosphäre. 
Beitr. Phys. d. fr. Atm. Bd. 26 (1940), S. 197; Theorie der stationären Leewellenströmung in freier Atmosphäre. 
Z. angew. Math. Mech. Bd. 23 (1943), S. 1. Für beliebige polytropische Schichtung findet sich die entsprechend allge- 
meinere Frequenzformel der Leewellen übrigens, wie mir Herr G. Lyra mitgeteilt bat, bereits bei P. Raethjen: 
Das zweidimensionale atmosphärische Stromfeld um ein Hindernis. Z. Flugtechn. Motorluftsch. Bd. 17 (1926), S. 185 
und bei F. Linke: Zur Vertikalvewegung isolierter Luftmassen. Meteor. Z. Bd. 45 (1928), S. 255. 
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$ 4. Experimentelle Bestätigung der Voraussagen von $ 3. 

Die gewählte einfache Versuchsanordnung®) zeigt Bild 1 (s. Bildtafel). Der Flüssigkeits- 
behälter (Bodenabmessungen 20,2% 7,6cm?) wurde mit einer Kochsalzlösung gefüllt, wobei 
wegen der Einfachheit der Herstellung eine lineare Dichteabnahme mit der Höhe zwischen 
den Grenzen 0,—=1 und maximal o,=1,19 (konzentrierte Salzlösung) gewählt wurde: 


2 


(2=0 Flüssigkeitsoberfläche bei Ruhe, o, Dichte an dieser Oberfläche, o, Dichte am Boden, 
H Höhe der Flüssigkeitsschicht). Diese kann wegen der Kleinheit der erreichbaren Dichte- 
differenz 0, —o, (=0,19) als hinreichende Approximation des Dichtegesetzes (3.1) 


angesehen werden, wenn hierin 
h=so, Hilo. - - - - 


gesetzt wird. (Diese Approximation ist um z2=0 am besten. Die Ermittlung des Winkels 
wird im folgenden im Gebiet der oberen Schichten vorgenommen werden.) 

Es wurde jeweils zunächst die Schicht geringster Dichte durch ein Zuleitungsröhrchen 
in das Gefäß eingelassen, daraufhin die nächst dichtere Schicht durch das Röhrchen langsam 
auf den Boden geleitet, wo es sich gleichmäßig unter die erste Schicht ausbreitete usf. 
Jeweils wurden etwa 10 Schichten gleicher Höhe eingebracht, und es wurde dann abgewartet, 
bis deren Grenzen sich durch Diffusion hinreichend verwischt hatten. Vor dem Versuch 
wurde dann die Dichteverteilung durch- Einlassen von Tropfen eines Gemisches von Xylol 
und Äthylenbromid in geeigneten Mischungsverhältnissen nachgeprüft, die jeweils in einer 
ihrer eigenen Dichte entsprechenden Höhe zur Ruhe kommen. 

Die zweidimensionalen Schwingungen wurden in der ersten Versuchsreihe durch die in 
Bild 1 sichtbare, vertikal gestellte rechteckige Platte als Störkörper verursacht, welche 
senkrecht von Wand zu Wand quer durch die Breite des Gefäßes reichte und leicht in die 
Flüssigkeit eintauchte (etwa 0,6 cm). Diese Platte wurde durch einen in der Skizze nicht 
wiedergegebenen Antriebsmechanismus in waagrechte harmonische Schwingungen senkrecht 
zu ihrer Ebene (und damit senkrecht zur Zeichenebene von Bild 1) versetzt. Als Amplitude 
dieser Oszillationen wurde zunächst 0,16 em gewählt. 

Bild 2 (s. Bildtafel) zeigt die beobachtete Erscheinung in einer Auswahl von sechs Schlieren- 
aufnahmen (Schattenaufnahmen, s. Bild 1) aus einer Versuchsreihe mit H=3,75 cm, o, = 1,00 
und 0,=1,19gcem” (h=19,7 cm), und zwar mit den dort angegebenen Frequenzen in auf- 
steigender Folge. Dabei ist hier und im folgenden 

Man erkennt die qualitative Übereinstimmung im Verlauf der als Schlieren sichtbar werdenden 
Charakteristiken mit den Aussagen des $ 3. In den letzten Bildern erkennt man auch die 
Reflexion der Charakteristiken am Boden. Entsprechende Reflexionen finden im weiteren 
Verlauf an der freien Oberfläche statt usw. Unmittelbar unter der Plattenkante ist als dunkler 
Kloß ein Wirbelballen sichtbar. Man’ hat sich das Bild der Platte und der Charakteristiken 
hin und her oszillierend zu denken. Im Bereich der Abweichungen der Bewegung von einer 
reibungslosen Bewegung (nächste Plattennähe und Bereich des Wirbels) unterliegen die 
Schlieren nicht den Gesetzen unserer Theorie. 

Für sehr kleine » verlaufen die Charakteristiken fast waagrecht und sind wegen der 
Kleinheit der Störung des Dichtegradienten als Schlieren schlecht zu erkennen. Auch wenn 
sich » andererseits dem kritischen Wert o, des Übergangs vom hyperbolischen zum elliptischen 
Verhalten nähert und die Charakteristiken fast vertikal verlaufen, werden sie undeutlich. 
Wenn sie auch noch für höhere Frequenzen zu verfolgen waren, als dies in der Bildfolge 2 
wiedergegeben ist, so waren sie aber kurz vor dem parabolischen Grenzfall nur in der Nach- 
barschaft des Körpers erkennbar. Wir geben in $5 eine Erklärung dafür. Auf Grund dieses 
Verhaltens ist es zweckmäßig, ©; mit größerer Genauigkeit extrapolatorisch zu ermitteln. 
Zur Auswertung der Beobachtungen wurden die Schlieren in Schattenaufnahmen, wie sie in 
Bild 2 wiedergegeben sind, für jede Schichtung und jede F requenz photographiert. _ 

Bild h (S. Bildtafel) gibt als Beispiel die Auswertung zweier Versuchsreihen (Versuchs- 
reihe Nr. 2 und Nr. 3 in der Übersicht von Zahlentafel 1) zur Nachprüfung des Gesetzes (I). 
Durch das Feld der jeweiligen Versuchspunkte wurde eine mittlere Gerade aus dem Ursprung 
heraus gezogen, welche in ihrer Verlängerung bis sing =1 den experimentellen Wert der 
kritischen Frequenz extrapolatorisch ergibt. 


% Die. Durchführung der Versuche hat Herr W. Müller übernommen. 
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Zahlentafel 1. Versuchsdaten. In allen Versuchen war o, = 1,00 g em”? 


Versuchsreihe Nr. 1 | 2 3 4 5 


| | 
70,0 37,5 50,5 70,0 | wo | | 
1,19 1,19 1,19 1,095 105 | 18 | 198 | 


In dieser Weise wurden mit der oszillierenden Platte als Störkörper insgesamt 7 lineare 
Schichtungen durchgemessen (Versuchsreihen Nr. 1 bis Nr. 7 in Zahlentafel 1). In einer 
achten Versuchsreihe wurde als Störkörper ein waagrecht liegender, bis zur Hälfte eingetauchter 
Kreiszylinder vom Durchmesser 2 cm benutzt, der senkrecht zu seiner Achse in horizontaler 
Richtung oszillierte. Die Ergebnisse sämtlicher acht Versuchsreihen sind in Bild 4 (s. Bild- 
tafel) in der Auftragung sing über w/w; (®; aus den Experimenten, wie oben beschrieben, 
ermittelt) wiedergegeben. Sie bestätigen innerhalb der Versuchsgenauigkeit gut das behauptete 
Gesetz (I). 

Die aus diesen Versuchsreihen ermittelten Werte der kritischen Frequenz »; der 
untersuchten Schichtungen sind in Bild 5 (s. Bildtafel) aufgetragen. Zum Vergleich ist der 
theoretische Verlauf von &; in Abhängigkeit vom Schichtungsparameter h als Kurve einge- 
zeichnet. Auch das Gesetz (II) wird somit befriedigend bestätigt. 

Wir haben uns überzeugt, daß diese Ergebnisse unbeeinflußt bleiben, wenn der Stör- 
körper nicht waagrecht sondern in einer beliebigen anderen Richtung parallel zu den Seiten- 
wänden des Gefäßes oszilliett. Auch die Schwingungsamplitude des oszillierenden Körpers 
kann in weiten Grenzen geändert werden, ohne die obigen Resultate zu beeinflussen. Für 
eine Amplitude von 0,03 bis 0,6 cm — soweit ließ sich bei der benutzten Antriebsvorrichtung 
die Amplitude vergrößern — haben wir dies für die Schichtung des Versuchs Nr. 1 bestätigt. 

Zum Schluß sei noch folgende Beobachtung kurz mitgeteilt, obgleich sie nicht unmittelbar 
in den Rahmen unserer Betrachtungen gehört. In Bild 6 (s. Bildtafel) sieht man neben den 
von der — hier tiefer eingetauchten — oszillierenden Platte ausgehenden und am Boden 
sowie an der freien Oberfläche reflektierten Charakteristiken eine Folge von stark ausgeprägten 
horizontalen Schlieren in Höhe des an der Plattenkante haftenden Wirbelballens. Während 
die Charakteristiken sofort mit dem Einsetzen der Oszillationen als Schlieren sichtbar werden, 
(aber natürlich erst nach einer gewissen Einschwingungszeit den von der Theorie geforderten 
Verlauf zeigen), bildet sich diese Erscheinung erst ganz allmählich vom Wirbel ausgehend 
nach links und rechts symmetrisch aus. Bild 6 gibt das Stadium der vollen Ausbildung 
wieder. Es handelt sich hierbei offenbar um eine Zähigkeitswirkung, und zwar um einen 
durch den hin- und herschlenkernden Wirbel verursachten Austausch der vertikalen Schichtung 
in der Höhenerstreckung dieses Wirbels. Darauf deutet die Folge der Streifen von oben 
nach unten: dunkel, hell, dunkel, hell, in der Schattenaufnahme Bild 6 hin, also: Zunahme 
des Dichtegradienten, dann Abnahme, dann wieder Zunahme und Abnahme. 


85. Verallgemeinerung der Ausführungen des $ 3 für beliebige stabile Schichtungen. 


Es ist bereits bemerkt worden (s. $ 3), dass die Theorie des $ 3 ganz entsprechend für 
ein isotherm geschichtetes, adiabatisch kompressibles Gas entwickelt werden kann. Qualitativ 
bleibt alles unverändert. Wir wollen im folgenden zu qualitativ neuen Ergebnissen fort- 
schreiten, indem wir nunmehr eine beliebige stabile Schichtung 90, (z) zugrundelegen. Das 
Medium sei wieder eine inkompressible Flüssigkeit. ’ 

Gegenüber den Ausführungen des $ 3 ergeben sich folgende Änderungen. An Stelle der 
zweiten Gleichung von (3.2) tritt . 


0... 82a). 


Die Koeffizienten des Systems (3.2) mit (3.2a) sind also im allgemeinen Fall nicht mehr 
durchweg konstant, sondern hängen zum Teil von z ab. Als Folge ergibt sich, daß die 
Grundgrößen X, Z und P nicht mehr alle der gleichen Differentialgleichung genügen. Es 
bezeichne D den Differentialoperator 

Dann ergibt die Elimination von X und P für Z die Differentiälgleichung 

Entsprechend genügen X und P je einer linearen homogenen Differentialgleichung zweiter 
Ordnung, bei der die Glieder der höchsten Ordnung durch D [X] bzw. D [P] gegeben sind, 
und zwar genügen X und P unter sich derselben Differentialgleichung, wie aus der ersten 
Gl. .. unmittelbar folgt. Für unsere Betrachtungen ist allein die Gestalt des Operators D 
wesentlich. 


8 

70,0 

1,19 


Bildtafel 1 zu: 


Görtler, Über eine Schwingungserscheinung in Flüssigkeiten mit stabiler Dichteschichtung. 


oszilherende Platte 


Bild 1. Versuchsanordnung. 


Mattscheibe 


q 


os’ 


Bild 3. Ergebnisse zweier Versuchsreihen als 
Beispiele für die Ermittlung von wg. 


«Versuchsreihe 


Bild 2. Charakteristiken, sichtbar gemacht als Schlieren, welche 

von einer leicht eingetauchten oszillierenden Platte ausgehen, bei 

verschiedenen Frequenzen: von oben nach unten ® = 0,435; 0,538; 

0,625; 0,714; 0,769; 0,848 see-t. Flüssigkeit: Kochsalzlösung mit in der 
Ruhelage linearer Dichteseh ichtung. 

(H =3,75cm, 9% und @=1,19gem-3, daher 19,7 cın.) 


Bild 6. Neben den Charakteristiken, welche sowohl am Boden als 

an der freien Oberfläche reflektiert werden, sieht man hier einen 

waagrechten Schlierenstreifen, der sich von dem an der Plattenkante 
haftenden Wirbel aus allmählich ausgebildet hat. 


| 
| 
05 w 10 
Bild 4. Ergebnisse von 8 Versuchsreihen (Versuchsdaten siehe 
Zahlentafel I, zur Bestätigung des Gesetzes (I) (ausgezogene Gerade). 
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0,5 
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cm 8 


Bild 5. Die experimentell ermittelten 
kritischen Frequenzen in Abhängigkeit vom 
Schiehtungsparameter h. Die ausgezogene 
Kurve stellt das theoretische Gesetz (II) dar- 


ter 
ıng 
nen 
ıng | | | 
ben & | | | | 
für 
der 
2a). 
ehr \ 
die 
Es 
a) 
ter 
nd, | 4 
ten 


| 


Bd. 23 Nr. 2 April 1943 Görtler, Über eine Schwingungserscheinung in Flüssigkeiten 


Es ist o,/o,<0 nach (2.2). Der Operator D ist nun für 
> elliptisch 
< hyperbolisch 


Es kann also nunmehr die grundlegende Differentialgleichung in dem flüssigkeitserfüllten 
Raum bei einer festen Frequenz ß ganz oder auch in einzelnen horizontalen 
Streifen, je nach der Schichtung, zu einem der drei Typen gehören. 
Die Charakteristiken bestehen aus den beiden Kurvenscharen 


Br ......... (86a). 


Sie sind in „hyperbolischen Streifen“ (#°?< g|o,jo,|) reell und verschieden. Die Tangenten 
der einzelnen Kurven bilden mit der positiven x-Achse einen Winkel p, für den 
sin @ (2) = () (Ia) 
gilt. Hierin ist , die im allgemeinen von z abhängige „kritische Kreisfrequenz“: 
Da der Charakteristikenverlauf bei jeweils festen $ nur von dem relativen Dichtegradienten 
0,jo, abhängt, liefern zwei Dichteverteilungen dann (und nur dann) dasselbe Charakteristiken- 
feld, wenn sie sich nur durch einen konstanten Faktor unterscheiden. Wenn in einem ganzen 
endlichen Streifen 2, <z<2, der parabolische Fall bei einer gewissen festen Frequenz 
vorliegt, so muß in diesem Streifen o,/o, = const. sein, d. h. o, muß dort notwendig dem 
speziellen Verteilungsgesetz (3.1) genügen. Im allgemeinen wird der parabolische Fall nur 
als Uebergangsfall auf horizontalen Grenzlinien zwischen hyperbolischen und elliptischen 
Gebieten auftreten. An parabolischen Grenzen hyperbolischer Streifen — immer bezüglich 
einer bestimmten festen Frequenz $ — münden die Charakteristiken senkrecht auf diesen 
Grenzen ein. 
Die Krümmung der Charakteristiken = (2) ist durch 


gegeben. Die Charakteristiken haben also in eventuell vorhandenen parabolischen Grenzlinien 
'hyperbolischer Streifen verschwindenden Krümmungsradius. 


“ elhptisch / Bild 7. Die spezielle Verteilung (3.1) besitzt die konstante kritische 
Vgl Frequenz Vgh. Die experimentelle Aproximation dieser Ver- 
teilung zeigt leichte Abweichungen. Demzufolge schwankt auch 
ihre kritische Frequenz 3% (2)=V—907/0, (ausgezogene Kurve, 
stark überhöht) leieht um den konstanten Wert Yg/h. Daher sind 
die von einer Stelle in Höhe 2=zou ausgehenden Schlieren für 
Frequenzen 3 nahe an V g/h nur für solche z sichtbar, für welche 
die Gerade 3 = const. in diesem Bilde noch nicht das schraffierte 
Gebiet betreten hat. (Siehe jedoch hierzu die Bemerkung im 
unten folgenden Text.) 


—>Zz 


Ich habe bereits in $ 4 erwähnt, daß in unseren Versuchen die Schlieren für 8-Werte 
nahe unter Yg/h nicht gut zu beobachten waren. Das überrascht nun nicht mehr. Denn 
wir hätten in unseren Versuchen die Verteilung (3.1) sehr sorgfältig approximieren müssen, 
um für Frequenzen f, die sehr nahe unter Yg/h liegen, nicht schon durch irgendwelche kleine 
Schwankungen von o, um die approximierte Verteilung (3.1) schmale elliptische Streifen 
verursacht. zu haben. Nach unserer Theorie, welche die Schlieren durch mathematische 
Unstetigkeitslinien wiedergibt, sind solche elliptische Gebiete absolute Sperren für diese 
Linien. In der wirklichen Flüssigkeit, in der die Schlieren örtlich starke Änderungen des 
Dichtegradienten angeben, trifft dies zwar nicht genau zu, aber die*Schlieren werden beim 
Eindringen in ein elliptisches Gebiet sehr rasch diffus (breiter und schwächer). (Es ist jedoch 
möglich, daß sie sich durch sehr schmale elliptische Streifen hindurch, aber nur entsprechend 
geschwächt, fortsetzen.) 

Der Vollständigkeit halber sei zum Schluß ergänzend bemerkt, daß im allgemeinen 
Falldreidimensionaler Schwingungen natürlich ganz entsprechende Verhältnisse vorliegen. 
Da die «-und y- Koordinaten gleichberechtigt auftreten, erscheint an Stelle des Operators D 


lediglich 

und an Stelle der Charakteristikenlinien (zylindrischen Charakteristikenflächen) treten ent- 
sprechend allgemeinere Charakteristikenflächen. 425 
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Strömungen über Bodenerhebungen auf der rotierenden Erde. 
Von W. Rothstein in Göttingen, z. Zt. im Felde'). 


(Aus dem Kaiser-Wilhelm-Institut für Strömungsforschung.) 


Die von H. Görtler aufgestellten Gleichungen für die Strömung einer 
volumbeständigen  Flüssigkeitsschicht konstanter Dichte über Boden- 
erhebungen auf der rotierenden Erde’); werden für den Fall endlich aus- 
gedehnter Bodenerhebungen behandelt. Dabei wird jedoch zunächst der in 
jenen Gleichungen erfaßte Reibungseinfluß vernachlässigt. Als Beispiele 
werden durchgerechnet ein Berg mit kreisförmigen und ein Berg mit 
elliptischen Höhenlinien. 


$ 1. Einleitung. 


H. Görtler hat in einer kürzlich erschienenen Arbeit?) den gemeinsamen Einfluß von 
Erddrehung und Bodentopographie auf den Verlauf einer stationären atmosphärischen 
oder ozeanischen Strömung unter Berücksichtigung der Reibung in Bodennähe untersucht. 
Das Medium wurde der Einfachheit zuliebe als volumbeständig und seine Dichte konstant 
angenommen. Nach oben soll es durch eine freie Oberfläche begrenzt sein. Von weiteren 
für das Zustandekommen solcher Strömungen in der Natur wesentlichen Einflüssen wurde 
bewußt abgesehen, um die mathematischen Schwierigkeiten in erträglichen Grenzen zu halten. 


Wir wollen zunächst die Voraussetzungen und die mathematische Formulierung des 
Problems skizzieren, auf Grund derer Görtler zu jenen Grundgleichungen gelangte, von 
denen wir auszugehen haben. Der Rechnung wurde ein erdfestes kartesisches Koordinaten- 
system x*, y*, 2* zugrunde gelegt, dessen Ursprung o mit einem Punkt der Erdoberfläche 
zusammenfällt, während die 2*-Achse radial nach außen weist. Es seien 


o= die Drehgeschwindigkeit der Erde, 


= {u*, v*, w*} die Geschwindigkeit eines Massenteilchens 


bezüglich der Erde, p* der Druck, o die als konstant angenommene Dichte, H* (x*, y*) die 
Höhe der Flüssigkeitsschicht und B ein Maß für die horizontale Ausdehnung der Boden- 
unebenheit. 


Der Einfluß der Reibung in der bodennahen Grenzschicht wurde in erster Näherung 
folgendermaßen erfaßt. Außerhalb der Grenzschicht, deren Rand durch 2* = {* (x*, y*) ge- 
geben sei, fließt die Hauptströmung reibungslos. In der Grenzschicht entsteht infolge der 
Schwächung der Corioliskräfte eine horizontale Sekundärströmung. Da diese im allgemeinen 
nicht quellenfrei ist, erfolgt aus Kontinuitätsgründen ein vertikaler Massenstrom zwischen 
Hauptströmung und Grenzschichtströmung. Man kann für den horizontalen Massenstrom M 
der Sekundärströmung setzen: M=oh®* Dabei ist die „gedrehte Geschwindigkeit“ 
= t- v*,u*,0% und h ein Proportionalitätsfaktor von der Dimension einer Länge, der auf 
der nördlichen Halbkugel positiv, auf der südlichen negativ ist. Für eine grob-abschätzende 
erste Näherung kann h als konstant angesetzt werden. 


Setzt man nun |Grad {*|<1 voraus?), so wird |w*| <{u*?+v*?}"2 sein. Vernach- 
lässigt man daher die mit w* behafteten Glieder in den Eulerschen Gleichungen und außer- 
dem wie üblich die vertikalen Beschleunigungen gegen die horizontalen, ersetzt man ferner 
die geographische Breite g durch @(0) und die Erdbeschleunigung g durch g(0) = {0,0, 9}, so 
erhält man die Bewegungsgleichungen 


1) Wegen seiner Einberufung zur Wehrmacht ist es dem Verfasser auf unbestimmte Zeit unmöglich geworden, 
gewisse geplante Ergänzungen zu den im vorliegenden Bericht wiedergegebenen Untersuchungen zu vollenden. Diese 
Ergänzungen beziehen sich auf eine Milderung der Voraussetzung bezüglich © in dem Konvergeuzbeweis des $ 4, auf 
eine Verschärfung der Abschätzung in $ 6 und ihre Verallgemeinerung auf beliebige Bergformen, sowie schließlich 
und vor allem auf eine nachträgliche Abschätzung des Reibungseinflusses, der auf Grund der Ergebnisse von 
H. Görtler (a.a. 0.) zunächst einmal vernachlässigt worden ist. 

2)H. Görtler: Einfluß der Bodentopographie auf Strömungen über der rotierenden Erde. Z. angew. Math. 
Mech. Bd. 21 (1941), S. 279 bis 303. Siehe auch den zusammenfassenden Bericht von H.Görtler: Neuere Beiträge zur 
Dynamik atmosphärischer und ozeanischer Strömungen. Naturwiss. Bd. 29 (1941), S. 473 bis 479, ferrer vom selben Ver- 
fasser: Strömungen über Bodenerhebungen auf der rotierenden Erde. Meteor. Z. Bd. 59 (1942), S. 4)5 bis 410. 


3) Wir benutzen mit Görtler die zweidimensionalen Differentialoperatoren Grad /= {fx,fu}. Div {a,b} =ax+by, 
Rot {a,b} =bx —a, (skalar!). 
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* 


mit sin g(o). Hinzu tritt die Kontinuitätsgleichung 


Hat die Hauptströmung beim Eintritt in den betrachteten Bereich eine von 2* unabhängige 
Horizontalgeschwindigkeit «*, »* — und dies soll fernerhin stets angenommen werden —, 
so behält sie diesen Charakter. Denn die Gln. (1.1) und (1.2) haben Lösungen der Form 
u* = u* (a*, y*); = v* (x*, y*); w* = a y*)-2*+b(x*,y). Man setzt daher für die gesuchte 
Geschwindigkeitsverteilung der Hauptströmung: 
— (x*,y*) + iw* (x*, y*, 2*) 
mit 
v*, 0), w* = {0,0, 

und berechnet ®* und w* gesondert. 

Es wird durch die Rechnung bestätigt, das H* bei genügend flachen Erhebungen {* im 


ganzen Strömungsbereich mit genügender Genauigkeit durch einen Mittelwert H ersetzt 
werden kann. Führt man dimensionslose Größen ein: z=«*|B, y=y*|B, z=z*/H, 
H, H=H*/H, wo B eine konstante waagerechte Länge, 
H eine senkrechte Länge (ungestörte Dicke der strömenden Schicht) und B}= {uf, v#, 0} eine 
örtlich und zeitlich konstante Bezugsgeschwindigkeit — im weiteren die ungestörte Anström- 
geschwindigkeit in der Hauptströmung weit vor der Bodenerhebung — bedeutet, und setzt 
man ferner R=|®}|/2»’ (der zu B} gehörige „Trägheitskreisradius“), so lassen sich die ersten 
beiden Gln. (1.1) nach Elimination des Druckes zusammenfassen zu*): 


Die Kontinuitätsgleichung (1.2) liefert bei präziser Durchführung des skizzierten Gedankens 


zur Erfassung des Reibungseinflusses in der Bodenschicht für die Hauptströmung die weitere 
Gleichung): 


wobei ® die gedrehte Geschwindigkeit bedeutet (vgl. S. 72). 
Aus (1.3) und (1.4) läßt sich in Verallgemeinerung einer Bemerkung von Görtler°), 


. in einfacher Weise nach Integration längs Stromlinien die Gleichung 


s 


so 
gewinnen. Das Integral ist über Stromlinien der ®-Strömung zu erstrecken, s ist die Bogen- 
länge dieser Stromlinien, C und s, sind längs Stromlinien konstant. 


Görtler hat der mathematischen Einfachheit halber zunächst den Fall vollständig be- 
handelt, daß £ nur von einer kartesischen Koordinate abhängt (langgestreckte Gebirgskämme). 
Sein wichtigstes Ergebnis war, daß die Reibung in der Bodenschicht bei den praktisch 


' interessierenden atmosphärischen und ozeanischen Strömungen im allgemeinen nur einen 


völlig untergeordneten Einfluß ausübt. Es ist nun von Interesse, wenigstens an einigen Bei- 
spielen das Verhalten der Strömung über Erhebungen z={ von endlicher Erstreckung quer 
zur Anströmungsrichtung zu untersuchen. Mathematisch liegt es nahe, etwa solche Er- 
hebungen zu wählen, deren Höhenlinien einfache geschlossene Kurven, z. B. Ellipsen sind. 
Görtlers Aussage über den Reibungseinfluß läßt demgegenüber die Vereinfachung ange- 
bracht erscheinen, daß man zunächst den Reibungseinfluß vernachlässigt. 


*%) H. Görtler: a.a.0. S. 302. 
5) H. Görtler:a.a.O., S. 302. 
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$ 2. Mathematische Formulierung des Problems. 
Es soll also in (1.4) und (1.5) h=0 gesetzt werden. Wir haben es dann mit dem System 


zu tun. Zur Festlegung der Grenzbedingungen beachten wir, daß unsere Gleichungen auch 
die Strömung über eine sich mit der festen Winkelgeschwindigkeit ©’ drehende unendlich 
ausgedehnte ebene Scheibe darstellen (vgl. Görtler a.a. O. S. 282). Haben wir es, wie wir 
annehmen, mit einer endlich ausgedehnten Bodenunebenheit zu tun, so ist es eine in gleicher 
Weise für unsere Strömung wie für die Modellströmung natürliche Forderung, wenn wir 


verlangen: lim ®=%, und lim Rot®=0. 
Dann folgt C= B/R und wir erhalten 


Auf Grund von (2.2) existiert eine in der ganzen x, y-Ebene eindeutige Funktion A (x, y), so daß 


Ar) 
(2.3) geht dann über in 


Div; Grad A= R 
Unsere Aufgabenstellung ist nun zurückgeführt auf die Bestimmung einer in der ganzen 
Ebene zweimal stetig differenzierbaren Lösung A der Gl. (2.4), von der lediglich noch ver- 
langt wird: 


Grenzbedingung: Im .. 2.2.2... 
>» 


(2.4). 


3. Einzigkeit Lösung. 


Gäbe es zwei Lösungen A" und A“ der Gl. (2.4), die der Grenzbedingung (2.5) ge- 
nügen, so würde für die Differenz A” = A" — 4'” gelten: 


4 
Die Kreisflächke 8={x’+y’< Ri} “mit der Peripherie R möge die Bodenunebenheit ganz 
im Innern enthalten.- Dann ist nach dem Gaußschen Satz: 


o- ; Grad Aw df= 


Im Gebiet gilt: 4A” =0. Wegen (32) ist in die zu A” konjugierte 
Potentialfunktion B eindeutig. Die Funktion F=A"+iB“ ist daher außerhalb 8 
eindeutig und analytisch und insbesondere im Unendlichen infolge der Grenzbedingung 


lim Grad A” —=0 regulär. Daraus folgt: 


(0) 


— 


Nun ist wegen (3.1): 


0= Div 


; Grad Aw = Div 4" i Grad (Grad A'")?. 


Also: 


A” A» 


Grad AU =0 (1-£>0). 


Zufolge (3.3) muß ds=0 sein, also nach (3.4): 


] 


th. Mech. . 
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8 4. Lösung der Gleichung (2.4). 
Gl. (2.4) lautet 


Wir schreiben: 
und machen für den Störungsanteil A’ nun den Ansatz 


Berg 


Dann ergibt (4.1) eine Integralgleichung für y (vgl. hierzu Courant-Hilbert, Methoden 
der mathematischen Physik II, S. 279ff.): 


Berg 


(4.4) - 


1 


Außerhalb des Berges — dort ist {=0 — verschwindet y. Daher ist jedenfalls nach be- 
kannten Sätzen über das logarithmische Flächenpotential unsere Grenzbedingung erfüllt, 
falls eine überall zweimal stetig differenzierbare Lösung y von (4.4) existiert. Es liegt nahe, 
für y die Neumannsche Reihe anzusetzen: 


(im Rahmen unserer Theorie 
der 1. Ordnung) 


|| w(än)d&dn 
Berg 


Gehen wir mit diesem Ansatz in (4.3) ein, so erhalten wir: 


4, mit: || (46). 


Berg 


Bezeichnen wir den Aufpunkt mit P(x,y), so ergibt sich zunächst formal für den uns 
allein interessierenden Gradp A’ die Entwicklung: 


Gradp A,= — || Gradp log 


Berg 


Konvergiert diese Entwicklung gleichmäßig, so liefert sie die Lösung unseres Problems. Die 
Neumannsche Reihe für y (und dann auch die Reihe (4.7)) konvergiert sicher gleichmäßig, 
wenn |Grad£| hinreichend klein ist. Man hat ja: r 


max 
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Diese sehr rohe Abschätzung könnte ohne grundsätzliche Schwierigkeiten wesentlich ver- 
bessert werden®). Vermutlich wird die Konvergenz der Entwicklung auch dann gewähr- 
leistet sein, wenn wesentlich schwächere Forderungen an |Grad {| gestellt werden. 


$ 5. Beispiel: Ein Berg mit elliptischen Höhenlinien. 
Wir führen elliptische Koordinaten ein: 
z==rcos#, y=-Yr—esind; ) 


eco?’ 


Flächenelement do = drdd 


Die Kurvenschar r = const. bestelit aus den konfokalen Ellipsen 
Ei(r): + 
Wir setzen nun voraus, die Bodenerhebung £ sei eine reine Funktion von r: 
Benutzen wir zur Abkürzung: 


so haben wir, wenn wir uns zunächst auf Grad A} als 1. Näherung für Grad A’ beschränken, 
zu berechnen: 


Ro 
\ r? e? cos? ’adar. 


Vr —e 


e 


Partielle Integration nach r ergibt bei Beachtung von k(R,)=0 


eu 


€ 


Hier müssen wir eine Bemerkung über die Anziehung homogen mit Masse belegter 
Ellipsen einschalten. Liegt der Aufpunkt (x, y) innerhalb der Ellipse E(r,) und setzen wir: 
L= \\ log r*do, 
Eirı) 
so gilt: 
I Yr-e 
Liegt dagegen der Aufpunkt (x, y), dessen „Radius“ r sei, außerhalb E(r,), gilt also E(r,)< E(r), 
so folgt aus (+) in Verbindung mit dem Satz von Ivory: 
ryri—e r,Vri —e 
Aus (++) liest man den Satz von Mac Laurin ab: Die Anziehungskräfte zweier homogen 
mit Masse belegter konfokaler Ellipsen auf einen äußeren Punkt sind gleichgerichtet und 
proportional der Gesamtmasse. (Zum Vorstehenden vgl. A. Wangerin, Theorie des Potentials 
und der Kugelfunktionen I, S. 19, 218.) 


Wir kehren zur Gl. (5.2) zurück. Den Aufpunkt bezeichnen wir jetzt mit P,(x,,y,) und 
es möge sein: &,=r,c0s®,, „=Yr; -—esin®,. P, liege in E(R,). Wir zerlegen: 


I,=2n- 


zur 


(++). 


Ro rı Ro 


e ra 
Ist nun F(r) der Flächeninhalt von E(r), so gilt auf Grund von (++) wegen r< r,;: 


6) Man könnte etwa (vgl. Courant-Hilbert a.a.0.) von (4.4) zur iterierten Integralgleichung übergehen, 
deren Kern K2 im Gegensatz zu K quadratintegrierbar ist. Für die Glieder der Neumannschen Reihe dieser 
Gleichung kennt man dann eine Abschätzung (vgl. G. Wiarda: Integralgleiehungen. Leipzig und Berlin 1930, S. 85 ff.). 
Jedoch würde eine solche Abschätzung praktisch recht mühevoll sein. 
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rı 


7 (r,+YVr; —e:) 


Ro 27 
Im Integral kann wegen r>r, nach (+) \ \ do durch x, ersetzt 
e 


ra e 


werden, also: 


Beachten wir schließlich: lim \ —=0, so folgt aus (5.3) und (5.4): 
Tı-ero rı 


dk —e 
e 
und entsprechend: 


Wir haben in diesen Ausdrücken Pet Aufpunkt wieder mit (x, ) bezeichnet. 
Im Sonderfall e=0 gehen die Ellipsen in Kreise über; r,® sind dann gewöhnliche 
Polarkoordinaten. Nach kurzer Rechnung ergibt sich für diesen Fall: 


0A _x 
r 


\r r' k(r')dr’; . . 
v 


und hieraus: 
r 


or 


Die Stromfunktion für die Störströmung lautet also beim Kreise: 
rı 
\ \ r,k(r,)dr,dr, 
r, 
v u 


Im allgemeinen Fall e>0 erhält man zunächst: 


2 2 
r co®®+Yr’ — e’sin 


r+yVr—e 


—24= 


Ro 


dk 
—— dr, + (r? — e?) sin? d \, dr 
r 


Mit Berücksichtigung der Beziehung und der Gleichung | yri 
erhält man aus (5.8): 


dk 
—e’dr, 


13 
\ 
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+Vr-edr, 


r 
r 
r dk 
\ r,Vri—e dr, 
e 
1 dk dk 
un 1 
+ \ r, <edr, dr, +y \ Cojdr, 
r e 
Ist nun ®» = {0,1}, so ergibt sich für die Stromfunktion: 
Die Stromlinien berechnen sich nach dem Schema der folgenden Gleichungen: 


G,(r)= const; 


@G=F,J,+ (r—e)J,— Js; 
G,=F,J, —k(r)+ F,J;; 


dk 
edr, 


Vri — e Ar Cof 


dr, 


e 
Wir haben folgende Zahlenbeispiele durchgerechnet: 
1. Kreissymmetrischer Berg: 


% 

—(1+cos2rr) für rs'],, 

0 für 


Bergdurchmesser B=200 km; Berghöhe [,=0,1. 
Anströmgeschwindigkeit a) |®,|=10 cm sec”' (also R= 1,37 km) 
b) |®,|=%0 cm (alsp R= 2,74 km). 

Das Ergebnis der Rechnung ist in den Bildern 1a und 1b wiedergegeben. Sie zeigen den 
Stromlinienverlauf in der Grundrißprojektion und darunter das Bergprofil im Aufriß. Bei 
der kleineren Geschwindigkeit erhält man einen geschlossenen Wirbel mit relativ geringer 
Geschwindigkeit. Insgesamt hat das Strömungsbild gewisse Ähnlichkeit mit dem Strömungs- 
feld um einen in eine Parallelströmung eingebetteten Einzelpotentialwirbel. Nach (2.3) 
handelt es sich bei unserem Problem um eine Verschmierung von Wirbeln über den Berg. 
Die den Rechnungen zugrunde gelegten Geschwindigkeiten interessieren ihrer Kleinheit ent- 
sprechend im Hinblick auf Meeresströmungen. Bei den wesentlich größeren Geschwindig- 
keiten atmosphärischer Strömungen sind die Ablenkungen entsprechend kleiner. 
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Bild la. 


Bild 1a. Strömung über einen kreissymmetrischen Berg von 200 km 
Durchmesser und der Höhe 0,1 H. Oben die Stromlinien und der Berg 
im Grundriß, unten das Bergprofil im Aufriß. Anströmgeschwindig- 
keit 10 cm seer1, 


>) 


Bild 1b. Entsprechende Darstellung für eine Anströmgeschwindigkeit 
von 20 cm see-1, 


Bild 2 (links). Strömung über einen Berg mit elliptischen Höhenlinien 
(Achsenverhälthis 1:3) Bergbreite 200 km, Höhe 0,1 H. Anströmgeschwin- 
digkeit 20 em see-1. 


2. Elliptischer Berg: 


für esr<r, 


für 


2 
Randellipse E(R,): =1 mit R,=15 und e=y2 (Achsenverhältnis also 1:3); 


B=200km; Berghöhe £,=0,1; |®%| = % cm sec”! (also Trägheitskreisradius R = 2,74). 
Im Bild 2 ist der Berg durch ein Aufrißprofil und durch die Höhenlinien {=0; 0,05; 
0,1 im Grundriß dargestellt. Die Stromlinien weisen wieder einen geschlossenen Wirbel auf. 


8 6. Abschätzung für die Reihe (4.7). 
Für die Gültigkeit der obigen Ergebnisse ist noch der Nachweis zu erbringen, daß wir 
im Rahmen unserer Theorie der ersten Ordnung die Reihe (4.7) nach dem Glied mit A; ab- 
brechen dürfen. Dabei beschränken wir uns auf den Sonderfall e=0. Für e>x gehen 
unsere Ergebnisse in die von Görtler für den zweidimensionalen Fall der Strömungen über 


unendlich lange Gebirgskämme über. Hier spielt in der Gl. (4.1) das Glied — Grad A-Grad{ 


keine ins Gewicht fallende Rolle. Es ist daher schon von vorneherein wahrscheinlich, daß 


auch für andere e-Werte dieses Glied im Rahmen unserer Genauigkeitsansprüche bedeutungs- 
los ist. 


Mit leicht verständlichen Abkürzungen ist 
Gradp A, +, (P})= | Gradp log r* (P, P')- ya+i(P')d op; 
erg 


er 


K(P', 
g 


’ 1 1 * ’ 
also: 
|| Gradp log r* (P’, P") y„ (P')d wpr. 
rg 


, | 

.9) | 
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Somit schließlich 


Gradp (P)J)=— 5, Gradp log r* (P, P’) Grad Grad» A, (P’')dwp. 


In unserm Falle ist Gradp {(P’)- Gradp A, (P) =0. Ferner er- 


halten wir 
N ı 
u 


r 27 
1 diödA, r—r'-cos(p — 

2a dr’ \= 


Ro In 
(Das \ ) verschwindet bekanntlich: Ein homogener Kreisring übt auf innere Punkte keine 


aus!). 


1 
— dr 
or r 


Gilt nun stets s <0, so erhalten wir: 


r 


wo P* einen mittleren Punkt auf dem Berg bezeichnet. Daher wegen: 


r 


max A, 


max), An ‚Berg Or 


Nehmen wir etwa = 0,1 an, so folgt: 


„|_ max 


_ max |dA,| | 
_max 1 

Berg, dr 91— 


_max 04, 1 

Berg! dr | 8° 
Diese Abschätzung ist sehr grob. Sie lehrt jedoch bereits, daß das von uns berücksichtigte 
Glied A, den wesentlichen Störungsanteil wiedergibt. Es ist bemerkenswert, daß bei dieser 
Abschätzung nur £,., und nicht |Grad Z| eine Rolle spielt. Zumindest beim Berge mit kreis- 
förmigen Höhenlinien ist die frühere scharfe Forderung |Grad {| < 1 also überflüssig, Man 
darf erwarten, daß dies auch im allgemeinen Fall gilt. 426 
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Zur Stabilität der elastischen Schalen |. 


Von Ernst-August Deuker in Hannover. 


In der vorliegenden Arbeit werden die Untersuchungen über endliche Ver- 
formungen eines Kontinuums und die aus ihnen abzuleitenden Stabilitäts- 
kriterien des elastischen Gleichgewichts, wie ich sie an anderer Stelle [3]') 
für den allgemeinen dreidimensionalen Körper angestellt habe, auf dünne 
elastische Schalen übertragen. Es werden deshalb zunächst die Grund- 
gleichungen der Schalentheorie für endliche Verformungen entwickelt, als 
notwendige Grundlage für Stabilitätsuntersuchungen. Die Heranziehung 
der Methoden und Ergebnisse der Differentialgeometrie, insbesondere die 
konsequente Anwendung der Tensorrechnung, ermöglicht nicht nur, all- 
gemeine Formeln von großer Durchsichtigkeit zu entwickeln, sondern er- 
leichtert auch die Behandlung spezieller Probleme. In dieser Arbeit wird 
das an dem Beispiel der Kreiszylinderschale, in einer weiteren an dem 
der Kugelschale gezeigt werden. 


I. Der Formänderungstensor. 


Eine Schale wird begrenzt durch zwei (auch im unverformten Zustand) gekrümmte 
Flächen, ihre Laibungen, deren Abstand, die Schalenstärke h, gegenüber den anderen Ab- 
messungen der Schale als klein vorausgesetzt wird. Diejenige zwischen den beiden Laibungen 
liegende Fläche, die überall die Schalenstärke halbiert, heißt Mittelfläche. Die geringe 
Schalenstärke ermöglicht folgende Annahmen (als hinreichende Annäherung) für den -Ver- 
formungszustand der Schale einzuführen, die wir nach’ W. Flügge®) geben: 


a) „Die bei der Verformung der Schale eintretende Verschiebung eines beliebigen Punktes 
sei durch die Verschiebung des zugehörigen Punktes der Mittelfläche derart bestimmt, 
daß die Normalen zur Mittelfläche auch nach der Verformung zur neuen Mittelfläche 
normal sind.“ 

b) „Die Spannungskomponente normal zur Mittelfläche sei überall so klein im Vergleich 
zu den anderen, daß ihr Einfluß auf die Formänderungen vernachlässigt werden 
kann.“ 


Gewöhnlich wird weiter noch gefordert, daß die Verschiebungen klein gegen die Schalen- 
stärke sind, um die .Linearisierung der Schalengleichungen zu ermöglichen. Da wir größere 
Verschiebungen nicht von vornherein ausschließen, sehen wir von dieser Annahme zunächst ab. 

Wir unterscheiden ein raumfestes und ein körperfestes Koordinatensystem. Das körper- 
feste System passen wir den Annahmen über die Verformung der Schale an. Es ist so ge- 
wählt, daß q,=0 die Gleichung der Mittelfläche ist. q, und q, sind geeignete Parameter 
auf der Mittelfläche. g, ist der Abstand eines Punktes von der Mittelfläche und läuft im 


Intervall >: sq,<s+ n Die Wahl des raumfesten Systems ist beliebig. 


Wir ne nun eine einparametrige Schar von Formänderungszuständen der Schale, 
die nach Voraussetzung bereits durch die einparametrige Schar der zugehörigen Mittelflächen 
bestimmt sind. Ist das raumfeste System @in rechtwinkliges kartesisches Koordinatensystem, 
so haben die Raumkoordinaten der Schalenpunkte als Funktionen der körperfesten Koordinaten 
und des Parameters t die besondere Gestalt: 


Die Funktionen 9; (go, t) bestimmen für jeden Wert von £ die Mittelfläche der Schale, Z; (go, t) 
sind die Komponenten des Einheitsvektors in Richtung der Normalen. Der Parameterwert t=0 
kennzeichnet die unverformte Schale. Ihre Koordinaten sind gegeben durch die Funktionen: 


= p; (96,0) Li (9,0) 


= (ge) + 4, &(ge) 


Das Quadrat des Abstandes zweier Schalenpunkte mit den Koordinaten ge, 9, bzw. go +dq®, 
q,+dg? ist gegeben durch das Bogenelement: 


. (1.1a). 


mit 
(14:9, = geo (91:1) — 24, (N + 
1) Die eingeklammerten Ziffern beziehen sich auf das Schrifttumsverzeichnis am Schluß der I. Teils der Arbeit. 
2) Vgl. [6]. 
6 
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"Im unverformten Zustand war die Entfernung derselben Punkte gegeben durch: 


mit 

(91,4) = 905 (9) — 29 #00 (M)- 

In den Gl. (1.2) sind die goo (q;.#) die Koeffizienten der ersten Fundamentalform der Mittel- 
fläche, die Komponenten des Maßtensors, die xo0 (g;,f) die Koeffizienten der zweiten Funda- 
mentalform, die kovarianten Komponenten des Krümmungstensors, und die coo (q4,8) die 
Koeffizienten der dritten Fundamentalform der Mittelfläche. Die coo lassen sich durch die 
kovarianten und gemischten Komponenten des Krümmungstensors darstellen und zwar gilt: 


Der Formänderungstensor hat nur drei von Null verschiedene Komponenten, nämlich: 


1 
€00 (94; t) = 5 (meo mb) 


1 1 
=> (ges — 966) — 4 +5 
£,=0; 2,=0 

Da vorausgesetzt worden war, daß die Schalenstärke h klein gegenüber den Hauptkrümmungs- 
radien der Mittelfläche ist, sind die (q,)” enthaltenden Anteile der moo und &00 klein gegen- 
über den anderen Termen. Sie werden deshalb in der Literatur meistens’ von vornherein 
nicht berücksichtigt. Wir wollen jedoch quadratische Glieder in q, zunächst beibehalten und 
nur an geeigneter Stelle von der Möglichkeit, sie zu vernachlässigen, Gebrauch machen. Wir 
berechnen das algebraische Komplement der moo, d.h. die kontravarianten Komponenten des 
Maßtensors. Sie eıhalten, wenn wir nach Potenzen von q, entwickeln und nach den quadra- 
tischen Gliedern abbrechen, die Werte: 


Die Determinante der Matrix (moo) bestimmt sich mit der gleichen Annäherung zu: 
. (1-44, H+q(4H°+2)K}. 


und daraus 
Dabei bedeuten: 
H= die mittlere Krümmung. 
K=|x|=x1x} das Gaußsche Krümmungsmaß der Mittelfläche. 


Beschränkt man sich nicht auf sehr kleine Formänderungen, so lassen sich dem kovari- 
anten Verzerrungstensor &g0 die zwei gemischten Tensoren zuordnen ‘): ‚ 


mie €)o 


== — ) q, (x4® DV — gie 


1 


II. Der Spannungszustand einer Schale. 


Es seien (90,9); s, (A,u=1,2) die gemischten Komponenten des sym- 
metrischen Spannungstensors, bezogen auf das körperfeste System der q., 9. Sind 
die kovarianten Komponenten der drei Spannungsvektoren auf den drei Ko- 


3) Vgl. [4] Bd. I, Kp. IV, $ 4. 4) vgl. [3]. 


und 
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ordinatenflächen q, =const und q, = const, so bestehen zwischen ihnen und den Komponenten 
des Spannungstensors die Gleichungen’): 


(2.1). 
Bezeichnet man ferner mit 


1 
(2.2) 


die kovarianten Komponenten des Einheitsvektors in Richtung der Normalen der Koordinaten- 
flächen qu.=const, qg,=const mit nu, n, die Komponenten der Normale einer beliebigen 
Schnittfläche, so ist (2.1) in der für beliebige Schnittflächen gültigen Gleichung 


(2.3) 


enthalten. 

Die geringe Schalenstärke erlaubt es, den Spannungszustand der Schale zu vereinfachen. 
Diese Vereinfachung besteht darin, daß man auf die Kenntnis der genauen Verteilung der 
Spannungen über die Schalenstärke verzichten und sich auf ihre Resultierenden und auf die 
resultierenden Momente in bezug auf eine in der Mittelfläche liegende Achse beschränken kann. 

Wir legen einen Schnitt durch die Schale senkrecht zur Mittelfläche und betrachten die 
jedem Punkt des Schnittes zugeordneten Spannungsvektoren. Um die in der ganzen Schnitt- 
fläche übertragene Kraft zu bestimmen, hat man die Summe der Spannungsvektoren, d. h. 
das Integral über die Schnittfläche, zu bilden. 

Wir machen uns von der Länge des Schnittes frei, indem wir die in dem Schnitt über- 
tragene Kraft auf die Längeneinheit der Schnittlinie mit der Mittelfläche beziehen. Die ver- 
bleibende Integration der Spannungsvektoren über die Schalenstärke ist jetzt ohne weiteres 
möglich, da diese sich längs der Normalen ohne Änderung ihrer Komponenten parallel in den 
zugehörigen Punkt der -Mittelfläche verschieben lassen. Man erhält z. B. für die Kompo- 
nenten der in einer Schnittfläche q, =const übertragenen Kraft pro Längeneinheit der g,- 
Linie die Darstellungen: 


Führt man in die Gl. (2.4) die Komponenten des Spannungstensors nach @.1) ein, so gehen 
sie über in: 


s) vgl. [3]. 
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Es ist Ym=yg{1—2q,H+qgK}, enthält also den von q, abhängigen invarianten Faktor 
(1-2, H+4K). 
Wir führen nun folgende Abkürzungen ein: 


h 
+3 


= | 


b) | (1-24, 
| 
In die Integrale (2.6) geht die spezielle Schnittrichtung nicht mehr ein, sie definieren auf der 
Mittelfläche im Falle 
a) einen Tensor, den Längskrafttensor, 
b) einen Vektor, den Querkraftvektor. 


Man beachte, daß N“ die kontravarianten Komponenten eines in jedem Punkt der Mittel- 
fläche definierten Flächenvektors sind, der aber keineswegs eine Kraft in der Fläche bedeutet. 


Durch die Größen T%; N“ sind die in einem Schnitt senkrecht zur Schalenmittelfläche 
übertragenen Kräfte pro Längeneinheit vollständig bestimmt. So erhält man z. B. für die in 
einem Schnitt q, = const wirkenden Kräfte (2.5) die Beziehungen: 


1 
aP;=N' vg" =N* 


. (25a), 


wobei „lu die kovarianten Komponenten der Tangentialnormale der Koordinatenlinie q, = const 
sind. Entsprechend gilt für einen beliebigen Schnitt senkrecht zur Mittelfläche: 


Die Gl. (2.7) sind das Analogon zu den GI. (2.3). 


Die Verschiebung der Spannungsvektoren längs der Normale in den entsprechenden 
Punkt der Mittelfläche führt zu Kräftepaaren mit dem Hebelarm q,. Wir bilden in ganz 
analoger Weise ihre Summen z. B.: 


> 
2 


+ 
h 


Zu dem Längskrafttensor 7% tritt also als zweiter Tensor: 


der Momententensor. 

Die in Richtung der Schalennormale fallenden Komponenten der Spannungsvektoren 
führen zu keinem Moment. 

Über die anschaulich zu deutenden Tensoren T%, M% und PR Vektor N" hinaus, wäre 
es durchaus möglich, weitere Größen zu betrachten, indem man die Integrale über die mit 
höheren Potenzen von q, multiplizierten Spannungskomponenten heranzieht. Diese Mittel- 
werte höherer Ordnung sind jedoch klein gegenüber 7“, M% und N“. Sie können deshalb 
vernachlässigt werden. 

Wir müssen nun noch einiges über die Spannungen in einer Schnittfläche parallel zur 
Mittelfläche, d. h. in den Flächen q,=const, ausführen. Die zugehörigen Spannungs- 


vektoren sind: 


a) 
TE 
(2.6) 
| 
d 
N 
D 
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Wir bilden die Integrale: 


b)N„= s3{1—-24, H+q}K}dg, 


Der Vektor mit den kovarianten Komponenten N, entspricht dem kovarianten Querkraft- 
vektor, weil (4 s3jX...ydq, gleich Null gesetzt worden ist. Man kann also schreiben: 


Der Mittelwert N,,, der Normalspannung s; spielt für die weiteren Rechnungen keine 
Rolle. Es ist jedoch N,, zweckmäßig aufzunehmen, um ein widerspruchsfreies Gleich- 
gewichtssystem für die oben eingeführten Spannungsgrößen der Schale aufzustellen. 

Wir sind bei der Definition der Längskräfte und Momente von den gemischten Kompo- 
nenten des Spannungstensors ausgegangen. Es wäre auch möglich gewesen, die kontra- 
varianten Komponenten des Spannungstensors heranzuziehen. Die mit ihnen gebildeten 
Integrale sind aber nicht identisch mit den kontravarianten Komponenten der Tensoren (2.6) 
bzw. (2.9); denn diese Komponenten sind mit dem Maßtensor g*“ der Mittelfläche zu bilden, 
während die s“* durch Vermittlung des Maßtensors m** aus den gemischten Komponenten s; 
hervorgehen, so daß bei der Integration abweichende Glieder auftreten. Der Längskraft- 
tensor ist deshalb im allgemeinen nicht mehr symmetrisch. Man hat nämlich: 


h 


2 


Tui Tre 2 (x? M} — M%) 


Die Abweichungen sind praktisch ziemlich unbedeutend, insbesondere beim Momententensor, 
für den ähnliches gilt. 


If. Die Gleichungen des Gleichgewichts. 


Ist der durch t=t, gekennzeichnete Formänderungszustand der Schale ein Gleich- 
gewichtszustand, so hat der Spannungstensor s/... den drei Differentialgleichungen 


ds, 
zu genügen®). Von der Aufnahme einer Massenkraft in den Gl. (3.1) sehen wir ab, da man 
sie gegebenenfalls durch eirme Oberflächenbelastung ersetzen kann. Das Symbol d bedeutet 
die kovariante Ableitung in bezug auf das körperfeste System mit der für t=t, maßgebenden 


Metrik. Schreibt man die kovarianten Ableitungen in (3.1) für das gewählte körperfeste 
System aus, so erhält man: 


_ Tiuss ym=0 


ds, 


=0, (3.1) 
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sind die Christoffelschen Dreizeigerzeichen zweiter Art im q,, 9, 9, System, 
deren Abhängigkeit von der Koordinate g, festzustellen ist. Nach Gl. (1.2a) folgt: 


wobei : n die Dreizeigerzeichen erster Art der Mittelfläche sind, und Kiu,e zur Abkürzung für 


gesetzt ist. Wir beschränken uns hier auf die in q, linearen Glieder, da wir in Ziffer II 
bei den Momenten stehen geblieben sind. Ferner gilt: 


Aus (3.3) und (3.3a) bestimmen sich die Dreizeigerzeichen zweiter Art zu: 


Wir setzen die vorstehenden Ausdrücke in die Gl.(3.2) ein und integrieren über die Schalen- 
stärke. Wegen des in jedem Glied auftretenden Faktors Ym führt die Integration auf die in 
(2.6) und (2.9) eingeführten Größen, deren zweckmäßige Definition hier deutlich in Erscheinung 
tritt. Die über die Schalenstärke integrierten Gl. (3.2) ergeben: 
yg) 


(3.4). 


Au +7 


Sie sind die Bedingungen dafür, daß die Resultierende aller auf die Mittelfläche bezogenen 
Kräfte verschwindet. Die noch fehlenden drei Momentengleichungen erhält man, wenn man 
die Gl. (3.2) mit q, multipliziert und wieder über die Schalenstärke integriert. Sie lauten: 


u 


b) — 
In den Gl. (3.5) und (3.6) erkennt man in den Termen von gleichartigem Aufbau die 
verjüngte kovariante Ableitung der Tensoren T{ und M} sowie die des Vektors N“), die 
wir wieder mit dem Symbol d bezeichnen wollen, bezogen aber jetzt auf das q,, q, Koordinaten- 
system der Mittelfläche. Der Klammerausdruck 
98° Kıu,g — 
erhält durch die Einführung der kovarianten Ableitung eine recht einfache Form. Es ist: 


Nach den Formeln von Mainardi-Codazzi gilt ferner®): 


(36). 


r 


?) Man kann auch statt der Tensoren 7%, M% und des Vektors N“ die Tensordichten T7 Vg, M% Vg und die 


Vektordichte N“ Vg einführen. 
‚» Vgl. [4] Bd. I, Kap. V,$ 4. 
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d 


b) MS + (sit...) 


dM% 
) 


d) M5 


Die Bedeutung der einzelnen Gleichungen ist leicht zu erkennen. Die zweı Gl. (3.72) 
sind die Verallgemeinerungen der Gleichgewichtsbedingungen eines ebenen Spannungs- 
zustandes, wobei die Differenzen der Schubspannungen auf den Laibungen der Schale den 
Massenkräften entsprechen. Es treten hier Zusatzglieder auf, welche die Komponenten des 
Momententensors enthalten. Sie haben ihre Ursache in den bei der Verbiegung der Schale 
zwischen den Flächen parallel zur Mittelfläche auftretenden Schubspannungen, deren Summe 
über die Schalenstärke nicht verschwindet, weil dem Ebenbleiben der Schnitte senkrecht zur 
Mittelfläche keine lineare Verteilung der Spannungen über die Schalenstärke entspricht. Die 
Gl. (3.7a) sind jedoch von den anderen Gleichungen auch ohne die Zusatzglieder nicht unab- 
hängig, sondern sie sind durch die Krümmung der Mittelfläche mit der Gl. (3.7b) gekoppelt, 
die das Gleichgewicht aller Kräfte in Richtung der Schalennormalen zum Ausdruck bringt. 
Die beiden Gl. (3.7c) dienen nur zur Definition der Querkräfte als Divergenz des Momenten- 
tensors und sind das Analogon zu der aus der Balkentheorie bekannten Beziehung. 

Die letzte Gl. (3.7d) steht unabhängig neben den übrigen fünf. Sie kann im folgenden 
unberücksichtigt bleiben, da die Normalkraft N, nicht weiter interessiert‘ Die Querkräfte 
lassen sich leicht eliminieren, so daß nur die drei wesentlichen Gleichungen 


h 
dx, +3 
1—2q,H ’K =0 


(3.8) 
+ 


verbleiben. Diese Gleichungen reichen für sich nicht aus, die Längskräfte und .Momente der 
Schale zu bestimmen. Es ist vielmehr. notwendig, das Stoffgesetz zu kennen, das Spannungen 
und Formänderungen miteinander verknüpft und aus dem dann der Zusammenhang zwischen 
den Längskräften, Momenten und den Verschiebungen der Schalenpunkte folgt. Wir betrachten 
deshalb jetzt die Formänderungsarbeit. 


b) 


IV. Die Formänderungsarbeit. 


Die dünne Schale ist ein Beispiel für einen elastischen Körper, dessen Punkte größere 
Verschiebungen erfahren können, ohne daß mit ihnen gleichzeitig die Verzerrungen groß werden. 
Wir wählen deshalb den Ansatz für die Formänderungsarbeit so, daß er dem in der eingangs 


erwähnten Arbeit angegebenen entspricht®). Es ist nur 
ersetzen. Man erhält so den Ausdruck: 


Für die & bzw. ’e$ hat man die Werte (1. 7) und (1.7a) einzusetzen. Der Integrand in (4.1) 
wird ein Polynom in q, mit Koeffizienten, die sich aus dem Maß- und Krümmungstensoren 
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der verformten und unverformten Mittelfläche aufbauen. Die Integration nach 4, läßt sich 
ausführen, wobei böhere als zweite Potenzen von q, im Sinne der vorhergegangenen Fest- 
setzungen zu vernachlässigen sind. Auch die Spannungs - Dehrungsgleichungen, die sich auf 
Grund des Ansatzes (4.1) ergeben, sind unmittelbar zu übertragen und lauten: 


Aus den Gl. (4.2) lassen sich in Verbindung mit den Gl. (1.7a) nach der in Ziffer II 
angegebenen Konstruktion die Beziehungen herleiten, welche die Längskräfte und Momente 
mit den Formänderungen verknüpfen. Die nach Einführung der Verschiebungen in (4.1) nach 
den Sätzen der Variationsrechnung aufzustellenden nielıt linearen Differentialgleichungen sind 
jedoch mit den heutigen Hilfsmitteln der Analysis unangreifbar. Es gelingt leider nicht, alle 
bei vorgegebener Belastung möglichen Gleichgewichtslagen zu überblicken, was das ideale 
Ziel einer Stabilitätstheorie ist, sondern man muß sich darauf beschränken, eine spezielle 
Gleichgewichtslage auf ihre Stabilität zu untersuchen. Nimmt man die Berechnung der Grund- 
lösung als erledigt an, so konzentrieren sich die Stabilitätsbetrachtungen auf die Untersuchung 
der zweiten Variation der potentiellen Energie, die für eine stabile Gleichgewichtslage positiv 
definit sein muß. 

Die Anwendung des Energieprinzips im Sinne der Variationsrechnung setzt ein nicht 
lineares Elastizitätsgesetz, z. B. in der Form (4.2) voraus, und man muß im Formänderungs- 
tensor mindestens noch die quadratischen Glieder in den Verschiebungen und ihren Ableitungen 
berücksichtigen. Es ist nun die Frage, ob es sich durch eine Annahme über die Grundlösung 
nicht vermeiden läßt, daß man die Formänderungen bis zu den quadratischen Gliedern in den 
Verschiebungen und ihren Ableitungen zu berechnen hat. Diese Überlegungen knüpfen zweck- 
mäßig an die Gleichungen des Gleichgewichts (3.8) an. 


V. Die Stabilitätsgleichungen. 


Unter dem kınfluß der Belastung habe sich in der Schale ein Spannungs- und Form- 
änderungszustand ausgebildet, der die Randbedingungen erfüllt. Er wird durch die Ver- 
schiebungen beschrieben, welche die Punkte der Mittelfläche erfahren. Die Zusammenhänge 
zwischen diesen Verschiebungen und den Tensoren T} und M{ (und damit T/ und MÄ{ selbst) 
seien bekannt. Ferner sind dann auch der Krümmungstensor und der Maßtensor der ver- 
formten Mittelfläche bekannt. Innere und äußere Kräfte sind im Gleichgewicht, d. h., die 
Gleichgewiechtsbedingungen (3.8) sind erfüllt. Wir betrachten nun einen benachbarten Spannungs- 
und Formänderungszustand, der durch den Parameterwert f,+ 4t gekennzeichnet werde. Soll 
auch jetzt noch Gleichgewicht vorhanden sein, so müssen die Gleichgewichtsbedingungen erfüllt 
sein für die an der Stelle £,+4At gebildeten Größen. Bildet man nun die Differenz ent- 
sprechender Gleichungen, dividiert durch At und geht zur Grenze At>0 über, so erhält man 
die Ableitungen der Gl. (3.8) nach t an der Stelle t=t,, die also bestehen müssen, wenn es 
zu einem Gleichgewichtszustand infinitesimal benachbarte Gleichgewichtslagen geben soll. 
Die 


(5.1) 


sind ‘drei homogene lineare Differentialgleichungen für die virtuellen Verrückungen, denen 
homogene Randbedingungen zugeordnet sind. Bekanntlich hat ein System homogener Diffe- 
rentialgleichungen mit homogenen Randbedingungen im allgemeinen keine von Null ver- 
schiedene Lösung, wenn die in ihm auftretenden Parameter fest vorgegeben sind. Nur wenn 
‚ die Lastparameter in (5.1) bestimmte Werte annehmen, die den Eigenwerten des Systems bei 
den vorgelegten Randbedingungen entsprechen, existieren von Null verschiedene Lösungen 
und damit benachbarte Gleichgewichtslagen. Wählt man das Stoffgesetz in der Form (4.2), 
so sind bei gegebener Grundlösung die Ableitungen der Längskräfte und Momente nach t als 
Funktionen der virtuellen Geschwindigkeiten festgelegt, und die Gl. (5.1) sind dann im wesent- 
lichen identisch mit den Jacobischen Differentialgleichungen, die der zweiten Variation der 
potentiellen Energie zugeordnet sind. 

Die auf Stabilität zu untersuchenden Grundlösungen haben nun in vielen Fällen die 
Eigenschaft, daß die zugehörigen Verschiebungen der Schalenpunkte und deren Ableitungen 
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so klein bleiben, daß sie mit den Grundspannungen linear nach dem Hookeschen Gesetz 
verknüpft werden können. Man denke z. B. an die geschlossene Kugelschale unter gleich- 
mäßigen äußeren Normaldruck. Hier bleibt die Gestalt der Schalenmittelfläche überhaupt 
erhalten, es ändert sich nur der Kugelradius, wobei aber höhere als erste Potenzen dieser 
Änderung bei der Berechnung der Grundlösung durchaus vernachlässigt werden können. Bei 
den Gleichgewichtsbedingungen (3.5) hat die Vernachlässigung höherer Potenzen der Ver- 
schiebungen zur Folge, daß die Veränderungen, welche die Dieizeigerzeichen und Krümmungen 
der verformten Mittelfläche gegenüber denen der unverformten Mittelfläche erfahren, nicht 
in Erscheinung treten. 

Es sind also die kovarianten Ableitungen auf die unverformte Mittelfläche, d.h, auf die 
für *=0 geltende Metrik zu beziehen. Es ist dann zweckmäßig, den Längskräften und Mo- 
menten der Grundlösung den Parameterwert t=0 zuzuordnen (in Abweichung von den früheren 
Festsetzungen, bei denen Formänderungen und Spannungen für ?=0 verschwanden). Setzt 
man nun in die variierten Gleichgewichtsbedingungen (5.1) t=0 statt ?!=t, ein, so bieten die 
neuen Stabilitätsgleichungen: 
olor; + 


a) 


(5.2) 


h 
b) dt dyrdge T, +% +(st.. 


gegenüber den Gl. (5.1) den Vorteil, daß man sich bei der Berechnung der Formänderungen 
auf die in den Verschiebungen linearen Terme beschränken kann, und die Variationen der 
Längskräfte und Momente mit den virtuellen Verrückungen nach dem Hookeschen Gesetz 
verknüpft werden. Die Stabilitätsgleichungen (5.2) können allerdings nicht mehr die Euler- 
schen Gleichungen eines Variationsproblems sein. Die vorstehenden Betrachtungen gelten 
natürlich nicht, wenn bereits der Grundlösung große Verschiebungen der Schalenpunkte zu- 
geordnet sind. Für Platten und Schalen kommen nun gerade solche Fälle als Stabilitäts- 
probleme wohl kaum in Betracht. Will man dagegen die nach Überschreiten der Stabilitäts- 
grenze mögliche Gleichgewichtslage berechnen, so muß man auf die ursprünglichen nicht- 
linearen Gl. (3.8) in Verbindung mit einem geeigneten Stoffgesetz, bzw. auf die Formänderungs- 
arbeit (4.1) zurückgreifen. Die Erweiterung des Ansatzes der Formänderungsarbeit auf end- 
liche Verschiebungen kann ja in verschiedener Art geschehen. Wieweit die neuen Gleich- 
gewichtslagen davon abhängen, ist eine andere Frage. Die Stabilitätsgrenzen jedenfalls, 
die man nach den verschiedenen Verfahren ermittelt, dürfen nicht wesentlich voneinander 
abweichen, wenn die Grundlösung noch dem gewöhnlichen Elastizitätsgesetz gehorcht, da 
dann sowohl die vereinfachten Stabilitätsgleichungen (5.2) wie auch die verschiedenen Er- 
weiterungen des Energieansatzes als gleichberechtigt anzusehen sind. Wir werden deshalb 
zur Kontrolle die Stabilitätsgrenze in einigen Beispielen einmal mittels der Gl. (5.2) und ein 
zweites Mal auf Grund des Ansatzes (4.1) berechnen. Es ist zu erwarten, daß die Differential- 
gleichungen, die der zweiten Variation der potentiellen Energie entspringen, etwas kompli- 
zierter sein werden, als die Gl. (5.2), da sich der Energieansatz nicht in dem Maße auf die 
besondere Eigenschaft der Grundlösung zuschneiden läßt. Aus den obigen Überlegungen geht 
noch hervor, daß auch der an den Vorschlag von E. Trefftz anschließende Ansatz '°): 


für die Formänderungsarbeit (statt des hier gewählten Ausdruckes) zum Ziele führen muß. 


Der Ansatz (4.1) erweist sich aber in manchen Fällen z. B. besonders bei der Kugelschale 
als vorteilhafter. 


VI. Der Zusammenhang zwischen den Formänderungsgrößen 
und den Verschiebungen. 

Um den Maßtensor und den Krümmungstensor der verformten Mittelfläche in möglichst 
einfacher Weise durch die Verschiebungen auszudrücken, welche die Punkte der Mittelfläche 
bei der Verformung erfahren, ist es zweckmäßig, ein raumfestes Koordinatensystem &;,x, (i=1,2) 
so zu wählen, daß es mit dem körperfesten System im unverformten Zustand zusammenfällt. 
Das Bogenelement hat in einem solchen System entsprechend Gl. (1.2) die Gestalt: 

mi 
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Die Gleichung der unverformten Mittelfläche ist dann. einfach: 


wo de = ist, aber keinen Tensor vorstellt. Die Funktionen bir ent- 
sprechen den Komponenten des Maßtensors bzw. Krümmungstensors der unverformten Mittel- 
fläche nach den Gleichungen: 


(6.3). 
bir(xh)) 6 
Die Raumkoordinaten der verformten Mittelfläche setzen wir in der Form: 
(64) 


an. Die Abhängigkeit des Maßtensors 990(g:,f) von den Verschiebungen ist nun leicht zu 
übersehen. Man erhält: 


dvk) dv, dv, 


(6.5). 


dv, dv, 


Wir entwickeln die Gl. (6.5) nach Potenzen der Verschiebungen und ihren Ableitungen, und 
brechen zunächst bei den linearen Gliedern ab, mit denen wir bei einer Stabilitätsuntersuchung 
nach Ziffer V auskommen. Die Gl. (6, , vereinfachen sich dann erheblich und.man erhält: 


Wir haben durch die ee in \ (6.6) bereits zum Ausdruck BERN daß die Ver- 
schiebungen v* in diesem Falle als kontravariante Komponenten eines Verschiebungsvektors 
im ursprünglichen g0-System aufzufassen sind. Die Normalverschiebung v, ist ein Skalar. 
Die Gl. (6.6) erhalten durch die Einführung der kovarianten Ableitung, die mit d&® bezeichnet 
Sei. um anzudeuten, daß sie sich’ auf das unverformte go-System bezieht, die einfache Gestalt: 


(0) p 
Ieo dp 


= (4) + 


Die Beziehungen zwischen dem Krümmungstensor und den zweiten Ableitungen der Raum- 
koordinaten x;, x, nach den Flächenkoordinaten g. erhält man am einfachsten über die Gauß- 
schen Ableitungsgleichungen der Differentialgeometrie‘'). Ist y;(i=1,2,3) ein beliebiges 
raumfestes System, Yy=Yyi(ge) die Gleichung der Fläche, so erhalten sie die Form: 


In den Gl. (6.8) bedeuten Zi’ die kontravarianten Raumkomponenten der Flächennormale, 
Tr die Dreizeigerzeichen des Raumsystems und er} die der Fläche. Andererseits gilt für 


die kovarianten Komponenten der Normale die Darstellung: 
Der Tensor e;xı verschwindet, wenn zwei Indizes gleich werden tınd ist gleich +1 für gerade 


Permutationen von 1,2,3, dagegen —1 für ungerade Permutationen mit &,,,=+1. Multi- 
pliziert man (6.8) und (6.9) skalar miteinander, so erhält man die Gleichung: 


die für die Berechnung des Krümmungstensors maßgebend ist. Die Gl. (6.10) werden wieder 
bescnders einfach, wenn man das spezielle Raumsystem mit der Metrik (6.1) benutzt. Be- 


1) [4]. 


XUM 


die 


Gli 


Ab 
ant 


deı 
fol; 


"beı 


Bil 


so 


ist, 


Die 
ver! 
ziel 


woh 
dies 
Gle 


Ag Are Deuker, Zur Stabilität der elastischen Schalen I 91 


schränkt man sich ferner auf lineare Glieder in den Verschiebungen, so erhält man für die 
Komponenten des Krümmungstensors der verformten Mittelfläche die Ausdrücke: 


0) 


Wir führen auch hier die kovariante Ableitung (®d ein und erhalten nach einer kleinen Um- 
rechnung: - 


+ 


0 
de om o ou on 


wenn man noch beachtet, daß die kovariante Ableitung eines Skalars mit der gewöhnlichen 


..(6.1la), 


Ableitung übereinstimmt, also auch 4 geschrieben werden kann. Aus (6.11a) ergeben sich 
mit Hilfe der kontravarianten Komponenten des Maßtensors: 


Au io uo 
(0) 7 - q — Io) >» +2 Vz; . . . . . . . . (6.7a) 


die gemischten Komponenten des Krümmungstensors zu: 


(0) 
(0 Id? (0), (0, (0), vr d 


Wir haben nun noch die Dreizeigerzeichen des verformten Systems bis zu den linearen 
Gliedern in den Verschiebungen zu bestimmen. Wir bilden zu diesem Zweck die kovarianten 


. (6.11b). 


Ableitungen an bezogen auf das verformte System, die gemäß der Definition der kovari- 
anten Ableitung alle identisch verschwinden. Setzt man nun für. die 9. die rechten Seiten 
der Gl. (6.7) ein, so hat man nur die > im neuen System kovariant abzuleiten. Bei den 
folgenden Termen dagegen ist die Ableitung d durch Wd zu ersetzen, da nur lineare Glieder 
"berücksichtigt werden sollen. Man erhält so die Gleichungen: 

Bildet man nun die Kombination 


2. 9%, 


so folgt, da 
dg’ 
ist, aus (6.12) die Gleichung: 


d q* d d d g* d 
Die in (6.13) auftretenden zweiten gemischten kovarianten Ableitungen sind nicht ohne weiteres 


vertauschbar, wenn die Mittelfläche keine euklidische M, ist. Es gelten vielmehr die Be- 
ziehungen: 


(6.13). 


wobei RX, die Komponenten des Riemannschen Krümmungstensors sind. Beachtet man 


diese Beziehungen und die Formeln von Mainardi-Codazzi, so lauten die endgültigen 
Gleichungen: 
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VII. Einige Stabilitätsprobleme der Kreiszylinderschale. 


Die in den Ziffern I bis VI entwickelten Grundlagen der Stabilitätstheorie dünner 
elastischer Schalen sollen nun durch einige Beispiele ergänzt werden. Der eigentlichen 
Stabilitätsuntersuchung hat die Bestimmung der Grundlösung, deren Eindeutigkeit festgestellt 
werden soll, vorauszugehen. Eine brauchbare Annäherung der Grundlösung wird in den 
meisten praktischen Fällen durch den in der Literatur als Membranlösung bezeichneten 
Spannungszustand gegeben. Um diese Membranlösung zu bestimmen, werden in den Gleich- 
gewichtsbedingungen, die in Näherung auf die unverformte Mittelfläche bezogen werden, die 
Querkräfte und Momente gleich Null gesetzt. Die restlichen Gleichungen reichen dann aus, 
die Längskräfte ohne Heranziehung des Elastizitätsgesetzes zu berechnen. Das willkürliche 
Nullsetzen der Momente steht allerdings häufig nicht im Einklang mit der sich einstellenden 
Verformung. Die Abweichung ist bei dünnen Schalen aber so gering, daß sie unberücksichtigt 
bleiben kann. Von größerer Bedeutung ist es, daß die Membranlösung sich im allgemeinen 
nicht den vorgeschriebenen Randbedingungen anpassen läßt, so daß sich ein reiner Membran- 
spannungszustand gar nicht ausbilden kann. Diese nicht zu unterschätzende Abweichung 
der Membranlösung von der wirklichen Lösung muß man in Kauf nehmen, um die Behand- 
lung des Stabilitätsproblems nicht außerordentlich zu erschweren. Wir wollen die Frage 
nach den Randbedingungen an einem Stabilitätsproblem der Kreiszylinderschale untersuchen. 

Eine kreiszylindrische Schale sei durch einen äußeren konstanten Normaldruck p [kgem”?] 
und durch eine pro Längeneinheit des Randes konstante Axialkraft P|kgem”'] beansprucht. 
Die Mittelfläche habe den Radius a. Berücksichtigt man zunächst nur den Normaldruck p, 
so wird unter seinem Einfluß die Mittelfläche in einen Kreiszylinder von etwas kleinerem 
Halbmesser übergehen, sofern die, Radialverschiebung überall ungestört erfolgen kann. In 
praktischen Fällen sind freilich die Enden des Zylinders meist durch Böden versteift, so daß 
die Ränder an der Verschiebung nicht teilnehmen. Die Folge davon ist, daß die Mittelfläche 
sich biegt, und außer den Längskräften auch Momente und Querkräfte auftreten. Entsprechend 
wird sich der Zylinder unter der in Richtung der Erzeugenden wirkenden Kraft P verkürzen 
und obendrein etwas ausweiten, was ebenfalls an den Rändern nicht ungestört vonstatten 
geht. Ein reiner Membranspannungszustand kann sich bei der angegebenen Belastung offenbar 
nur einstellen, wenn die Grundlösung keine Randbedingungen für die Verschiebungen zu er- 
füllen braucht. Diese Grundlösung ist leicht anzugeben. Den Zusammenhang zwischen einem 
rechtwinkligen kartesischen «, y, 2-System und dem der Schale angepaßten Zylinderkoordinaten- 
system x,,%,,x, wählen wir in der Form: 


Be 
zie 


fal 
rec 


 z=a(l+x,)cosx, 


Die Gleichung der unverformten Mittelfläche lautet in den Zylinderkoordinaten &,, x,, ®;: 


de: 


bzw. im x, y, 2-System: 


Measing; - . . 0. : (71b). 


Die Komponenten des Maßtensors und des Krümmungstensors erhalten die Werte‘ 


(7.2) 
1 0), 0), 


Die Dreizeigerzeichen verschwinden sämtlich, und die Gleichgewichtsbedingungen (3.8), be- 
zogen auf die unverformte Mittelfläche, lauten einfach: 


10 10 h 
=0 
Der Index (% bei den Längskräften und Momenten soll andeuten, daß es sich um die 
entsprechenden Größen der Grundlösung handelt. Die Membranlösung der Gl. (7.3) lautet in 
unserem Fall: 


u 
wc 
M« 
ve 
“ein 
Ka 
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(7.4), 


wobei a ohne merkbaren Fehler gegenüber 1 vernachlässigt werden kann. 


Wir behandeln das Stabilitätsproblem nach der in Ziffer V auseinandergesetzten 
Methode, und können uns deshalb auf die Variationen des Maß- und Krümmungstensors der 
verformten Mittelfläche beschränken. Nach den Gl. (6.7), (6.11) und (6.15) erhält man: 


69, — ÖV,) 
5 öx.,=aör, 


1 


In den Gl. (7.5) bedeuten, entsprechend einer in der Schalentheorie vielfach üblichen 

Bezeichnungsweise, Punkte Ableitungen nach g, und Striche Ableitungen nach q,. Die Indi- 

zierung der Verschiebungen v, und v, hat hier keine Bedeutung im Sinne der Tensorrechnung. 

Die Variationen der beiden Formen e$ und ’&; des gemischten Formänderungstensors 

fallen für {=0 zusammen und seien mit ö*el bezeichnet. Aus den Gl. (1.7) bzw. (1.7a) be- 
rechnen sich die ö*e% für die Kreiszylinderschale zu: 


1 
+ 


2 


ag, 6%, 
Benutzt man als elastische Konstante den Elastizitätsmodul E und den reziproken Wert 
der Querdehnungszahl „_L, und führt noch die Abkürzungen 
Eh 
ein, so erhält man für die Komponenten von Längskraft- und Momententensor nach der 
Konstruktion von Ziffer II die Ausdrücke: 
T!= Dido, v ka? ö%,} 


öM}=Dka(l —») — — di 


öM;}=Dka(ll—») 


Dka {öv,+röo, — 
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Die Stabilitätsgleichungen, die der Grundlösung (7.4) zugeordnet sind, lauten: 


2 .— > 0 (7.8). 
2 
ı | 


Für die 677; ö8M+= gi) öM} sind die Ausdrücke (7.7) einzusetzen. Man erhält drei 
partielle Dilwentisigleichungen, linear und homogen in öv,, öv, und öv,, die, wenn man 
noch die dimensionslosen Größen 


einführt, folgende Gestalt annehmen: 


1+v 


3 
5 


a) + 5 


1 „ 


— 


++ 9, — 0) 8%, —k 


— +85, — — + — 0,084, —0 

Vergleicht man die Gl. (7.9) mit den entsprechend@n Gleichungen in den Arbeiten von 
W. Flügge und K. v. Sanden-Tölke"), so zeigt sich völlige Bern in den 
von k, o, und 0, freien Gliedern, wie auch nicht anders zu erwarten ist. Die Abweichungen 
in den mit k multiplizierten Termen sind auf die z. T. anders lautenden Definitionen der 
Schnittgrößen zurückzuführen. Die o, und o, enthaltenden Anteile weichen in allen drei 
Systemen mehr oder weniger voneinander ab. Auf die Größe der Knicklasten haben diese 
Abweichungen jedoch, wie wir sehen werden, keinen nennenswerten Einfluß. 

Da man in der technischen Literatur vielfach nicht .scharf genug die ursprünglichen 
Differentialgleichungen von den Variationsgleichungen (Eulersche und Jacobische. Diffe- 
rentialgleichungen im Sprachgebrauch der Variationsrechnung) unterscheidet, so stehen z. B. 
in den Differentialgleichungen, die K. v. Sanden und F. Tölke angeben, auf der rechten 


Seite noch die Belastungen der Schale, so daß die Gleichungen inhomogen sind, während in , 


Wirklichkeit die Variationen der Belastungen einzusetzen sind, die, sofern sie nicht ver- 
schwinden, linear und homogen in den virtuellen Verschiebungen sind. 

Es ist nun die Frage zu entscheiden, welche Randbedingungen den virtuellen Ver- 
schiebungen aufzuerlegen sind. Der einfachste Fall ist der, daß entsprechend den Rand- 
bedingungen für die Grundlösung (7.4) auch die virtuellen Verschiebungen am Rande nicht 
vorgeschrieben sind. Es ist aber auch möglich, die Schar der virtuellen Verrückungen durch 
geeignete Randbedingungen einzuschränken. Wählt man diese Randbedingungen so, wie sie 
gewählt werden müßten, wenn bereits für die Grundlösung bestimmte Randbedingungen vor- 
geschrieben sind, so werden die erhaltenen Kuisklasten ungefähr denen entsprechen, die der 
wirklichen Grundlösung zugehören. 

Die Differentialgleichungen (7.9) gestatten Produktansatz: 


89, Asinng,-sinZq, 
7 


wobei n ganzzahlig sein muß, und der für i=m'T (m ganzzahlig) bereits für sich mecha- 


nisch sinnvolle Randbedingungen erfüllt. Führt man die Funktionen (7.10) in die Differential- 
12) Vgl. [5], [6]. [8]- 
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gleichungen (7.9) ein, so erhält man drei lineare homogene Gleichungen für die Koeffizienten A, 
Bund 


+4) — n’]—o, (m — 1) — 9,4, C0=0 


Die drei Gl]. (7.11) haben nur dann eine von Null verschiedene Lösung, wenn die Deter- 
minante des Systems verschwindet, was einer Bedingungsgleichung zwischen den Last- 
parametern o, und o, gleichkommt. Bei der Berechnung der Determinante können Glieder, 
die von höherer Ordnung klein werden, insbesondere höhere Potenzen von o,, 0,, k und ihre 
Produkte, vernachlässigt werden. Man erhält so: 


In der vorstehenden Gleichung sind offenbar die Glieder 2k (1— »?)A*, o,»A* und o,vi* 
gegenüber (1 —»*)A* bedeutungslos. Vernachlässigt man sie, so stimmt die Gl. (7.12) völlig mit der 
von W. Flügge angegebenen überein. Da die Knickgleichung dort auch eingehend diskutiert 
worden ist, erübrigt es sich, weiter darauf einzugehen. Dem Wert A=0 entspricht eine aus- 
gezeichnete benachbarte Gleichgewichtslage. Sie ist allein in Wirklichkeit diejenige, die der 
Sa (7.4) zugeordnet ist. Die RR reduziert sich für A=0 auf: 


k(n? — 

Die kleinste Knicklast gehört zu dem Wert n=2 mit og, =3k. Mechanisch ist dieser 
Fall dadurch ausgezeichnet, daß für ihn die Variation des Maßtensors der Mittelfläche ver- 
schwindet, d. h. der Übergang in die benachbarte Gleichgewichtslage ist mit keiner Dehnung 
der Mittelfläche verbunden. Der kritische Außendruck nähert sich asymptotisch dem Wert 
0,—=3k, wenn zwar die Böden des Zylinders versteift sind, die Länge des Zylinders aber 
sehr groß ist. Man spricht deshalb im Falle A=0 auch von der Knicklast des unendlich 
langen Zylinders. 

Wir wollen nun dasselbe Problem mit Hilfe des Ansatzes (4.1) der Formänderungsarbeit 
behandeln. In den beiden gemischten Formen ’e und & des Verzerrungstensors kann man 
sich bei den q, enthaltenden Termen auf Anteile beschränken, die in den Verschiebungen 
linear sind. Nur die Verzerrung der Mittelfläche ist bis zu den quadratischen Gliedern zu be- 
rechnen, entsprechend der Tatsache, daß die rn praktisch biegungsfrei ist. 

Wir betrachten diesen Teil der ’e; und &{ und den Beitrag, den er zur Formänderungs- 
arbeit liefert, zunächst für sich. Die Formänderungsgrößen der Mittelfläche sind: 


. (7.12a). 


jeo= (& )aa == (Be — 9%.) (7.13). 


Die neuen Komponenten des Maßtensors der Mittelläche lassen sich, ausgedrückt durch 
die jo, somit schreiben : 


Für die kontravarianten Komponenten ge gilt die Identität 
gi — | 
die man als Iterationsformel für ihre Berechnung benutzen kann. Um in den Formänderungs- 


komponenten j, alle Glieder zu erfassen, die in den Verschiebungen quadratisch sind, ge- 
nügt bereits die erste Iteration, d. h. man kann setzen: 


(7.14), 
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Man gewinnt damit für die zweite gemischte Form des Verzerrungstensors die Darstellung: 


.0 


In die Formänderungsarbeit gehen die Ausdrücke 'j/ j, und ’ee ein, die sich berechnen zu: 


(7.16). 


Wir spezialisieren nun die vorstehenden Formeln auf den Kreiszylinder. Die ’j; erhalten die 
Werte: 


(7.17). 


Bei den Formänderungskomponenten haben wir uns mit den quadratischen Termen begnügt 
und beschränken uns deshalb in dem Ausdruck für die Formänderungsarbeit auf Glieder 
dritter Ordnung. Wir setzen also: 


+ 


-2(,+ +v,) 


(®, + v,)’ + (a? + + 


2 u? 


Nach geeigneter Zusammenfassung und unter Aufnalıme der” Krümmungsänderungen er- 
hält man für die Formänderungsarbeit des Kreiszylinders folgenden Ausdruck: 


1 
A=Gh \ \ 
v v 


m—i ((ö, a? (®, +,,+ v,) a? v,) 


1% 


(7.18). 
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Wir nehmen an, der konstante Außendruck habe stets die Richtung der Flächen- 
normalen, dann sind seine kovarianten Komponenten gegeben durch: 


Er leistet bei einer virtuellen Verrückung öxi die Arbeit: 
54,= | . . (7.20). 
00 00 


Die kovarianten Komponenten der Flächennormale £; sind nach den Gl. (6.9) zu berechnen. 
Für den Kreiszylinder nehmen sie in den hier benutzten Koordinaten die mit Yg multi- 
plizierten Werte: 


v9 a? (1 + v,) (di + | 


v9 ,=a’(l+v,) — v, 
an. Setzt man diese Werte in (7.20) ein, so erhält man den Ausdruck: 


2a 
2.7.2). 


+ 07 + + + (di + v,) ö, dv,] ad d 42 
Er läßt sich als die erste Variation des Integrals: 
27 


00 
ansprechen, da bei der bekannten Umformung der ersten Variation die Randglieder ver- 


schwinden. Dem Außendruck läßt sich also ein Potential U, = — A, zuordnen. 
Der Grundlösung (7.4) entsprechen, unter Vernachlässigung höherer Potenzen von 


und die Verschiebungen: 
Di-—»)’ ° D(1-— »*) 
(7.4a). 
’ 
v\” und vi sind bis auf unbedeutende Glieder höherer Ordnung Lösungen der Euler- 


schen Gleichungen, die dem Variationsproblem 
zugeordnet sind. 

Berechnet man nun die zweite Variation der potentiellen Energie U=A+T, für die 
Grundlösung (7.4a) und bezeichnet die Variationen von ®,, v,, vd, mit u, v, w, so erhält man 
den Ausdruck: 

1 
”u=Gh \ \ kat + m)? 


2 
m—1 
00 


D 


+ +2a? +ö)(2W 4-84) 


a? 
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+0? (am? — 2wov)+ +m+a?w)’ 2(d +5 
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(7.24). 
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Die Jaeobischen Differentialgleichungen, die sich der zweiten Variation zuordnen 
lassen, erhalten hier die Gestalt: 


it», 
2 

+ 


1% +0, 


1—»? 


1+v 
2 
J 

In den Gi. (7.25) fällt zunächst auf, daß die dritte Gleichung auch in den von o, und 
0, freien Gliedern von der Gl. (7.9c) abweicht. Das hat folgenden Grund: Bei der Ein- 
führung der Spannungsresultanten haben wir uns'mit den Resultierenden zweiter Ordnung, 
d. h. mit den Momenten, begnügt. 

Da im Formänderungstensor und bei der Durchführung der Integrationen über die 
Schalenstärke die zweiten Potenzen von q, vollständig berücksichtigt worden sind, so hätte 
man, um diese Glieder von den Spannungsresultanten aus sämtlich zu erfassen, auch noch 
die Resultierenden dritter Ordnung einzuführen. Ihr Einfluß ist jedoch praktisch ohne Be- 
deutung, und es lohnt sich nicht, sie in die Gleichgewichtsbedingungen aufzunehmen. Beim 
Ansatz der Formänderungsarbeit werden alle in q, quadratischen Glieder erfaßt, die Gl. (7.25) 
zeigen daher in den entsprechenden Termen eine größere Symmetrie als die Gl. (7.9). 

Diese Unterschiede, wie auch die erwarteten Unterschiede in den Knickgliedern machen 
sich jedoch bei der Berechnung der Knickdeterminante nicht geltend. Innerhalb der Glieder, 
die in die Gl. (7.12) aufgenommen wurden, zeigt sich kein Unterschied. 

Wir wollen nun noch kurz auf die Frage eingehen, wann ein zylindrisches Rohr unter 
reiner Torsionsbeanspruchung ausknickt, weil es sich zeigt, daß eine spezielle Lösung dieser 
Aufgabe, die in der Literatur gelegentlich als Greenhill-Fall bezeichnet wird, nicht existiert. 

Die Grundlösung für den tordierten Kreiszylinder lautet: 


(0) (0) 
oonst=9. 
Die Momente M ! und M: sind unbedeutend und können vernachlässigt werden. Die 
Stabilitätsgleichungen stellen wir im Anschluß an die Überlegungen in Ziffer V auf. Sie 
lauten für die Grundlösung (7.26): 
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Während bei der vorigen Aufgabe die Grundlösung nicht die Randbedingungen eines 
Zylinders mit ausgesteiften Böden erfüllte, es aber gelang, eine einfache Lösung der Stabilitäts- 
gleichungen mit geeigneten Randbedingungen für die virtuellen Verschiebungen anzugeben, 
liegt beim Torsionsproblem die Sache gerade umgekehrt. Es bleibt nämlich ohne Einfluß auf 
die Grundlösung (7.26), ob an den Rändern die Schubkräfte vorgeschrieben sind, oder ob man 
einen Rand festhält und die Verdrehung des anderen Randes vorgibt. Dagegen macht es 
sehr große Schwierigkeiten, eine Lösung der Stabilitätsgleichungen (7.27) zu konstruieren, bei 
der die virtuellen Verschiebungen die Randbedingungen erfüllen, die sich etwa ergeben, wenn 
bei der Grundlösung die Verdrehung eines Randes vorgeschrieben ist. Die Differential- 
gleichungen (7.27) gestatten nämlich nicht mehr den Produktansatz (7.10), sondern nur noch 
den allgemeineren Ansatz: 


öv, = Asin (na, 


30, = Basin[ng — 4.) 


cos 


Geht man mit diesem Ansatz in die Differentialgleichungen hinein, so erhält man für die 
Koeffizienten A, B und (© die drei linearen homogenen Gleichungen: 


0=0 
B 


k[ + v(m 


2)’ — 
Setzt man wieder die Determinante des Systems gleich Null, so ergibt sich mit den gleichen 
Vernachlässigungen wie vorher die Gleichung: 
wobei die rechte Seite mit der rechten Seite der Gl. (7.12) übereinstimmt. Die Gl. (7.30) ist 
bei gegebenem n als eine Gleichung achten Grades für A anzusehen, deren Wurzeln acht Teil- 
lösungen von der Form (7.28) bestimmen, aus denen dann die Lösung für bestimmte Rand- 
bedingungen aufzubauen ist. Diese Rechnungen sind nur numerisch durchzuführen und wegen 
ihres Umfanges bisher noch nicht in Angriff genommen. Nur wenn keine Randbedingungen 
für die Verschiebungen vorgeschrieben sind, oder diese, ähnlich wie bei der vorigen Aufgabe 
im Falle des reinen Manteldrucks, für einen sehr langen Zylinder außer acht gelassen werden 
können, stellt die Gl. (7.30) bereits eine Bestimmungsgleichung für den Lastparameter o dar, 
in der A und n so zu wählen sind, daß o möglichst klein wird, wobei n ganzzahlig sein muß, 
i aber kontinuierlich veränderlich ist. Wir beschränken uns auf diesen Fall und betrachten 

zunächst den oben erwähnten Sonderfall, der zu dem Wert n=1 gehört. 
Der Koeffizient von o in (7.30) unterscheidet sich von dem von W. Flügge ausgegebenen 
Wert durch das Glied (1+»)4°. Für größere Werte von n hat die Abweichung keine Be- 
deutung, für n=1 ist sie aber wesentlich. Setzt man nämlich n=1, so erhalten wir: 

2(#?—») 

Die weiteren Schlüsse, die an Hand der Gl. (7.3la) von W. Flügge und in ähnlicher 
Weise auch in der Arbeit von K. v. Sanden-F. Tölke gezogen werden, sind folgende '*): 


13) Vgl. [6], [8]. 


(7.31) 
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Die rechte Seite der Gl. (7.31a) erhält offenbar einen kleinen Wert, wenn A klein ist. Dann 
können höhere Potenzen von A vernachlässigt werden, und man erhält einfach: 


Der kleinste Wert, den A annehmen kann, wenn die Lösung periodisch mit der Periode / 
Zylinderlänge) sein soll, ist und damit: 


(7.32). 


3 
Das kritische Torsionsmoment M=2ra?Do=2n a. R hat dann die gleiche Größe, wie 
das kritische Torsionsmoment eines geraden Stabes mit der Steifigkeit Ena*’h. 

Bei der Gl. (7.31) ist dagegen die obige Schlußweise nicht möglich, es wäre also aus 
ihr zu schließen, daß die Greenhill-Lösung nicht existiert. Bevor wir auf diese Frage 
näher eingehen, wollen wir feststellen, wie man bei n >1 vorgeht. Hier verläuft die Be- 
handlung der Gl. (7.30) entspreehend der vorher gemachten Bemerkung, daß nämlich die 
Werte von n und Ä aufzusuchen sind, für die o ein Minimum wird, wobei für n die nächst- 
liegende ganze Zahl zu nehmen ist. Dann zeigt sich, daß nur n=2 in Frage kommt. Setzt 
man n=2 in (7.30) ein, so erhält man: 

Ven den beiden in Frage kommenden Wertbereichen liefert der Bereich A<1 einen 


weitaus kleineren Wert von o als der Bereich A>1. Für A<1 können die höheren Potenzen 
von ./ vernachlässigt werden, so daß man erhält: 


Die erwähnte Abweichung spielt keine Rolle. 
do 1-—»’ 3k 


u’ - pri 
und-die zugehörige Knicklast ist demnach gegeben durch: 


Sie ist unabhängig von der Zylinderlänge, was auch der Fragestellung entspricht. Der Para- 
meter o hat dagegen fürn=1 als Funktion von A kein Minimum an der Stelle A=0, sondern 


nur eine einfache Nullstelle. 


Die abgeleiteten Formeln werden in einer weiteren Arbeit zur Behandlung der Stabilität 
der Kugelschale unter gleichförmigem Außendruck verwandt werden. 
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Zusammenfassender Bericht. 


Die konfluente hypergeometrische Funktion 
mit besonderer Berücksichtigung ihrer Bedeutung für die 
Integration der Wellengleichung in den Koordinaten eines 
Rotationsparaboloides’). 
Von Herbert Buchholz in Berlin. 


(Mitteilung aus dem Zentrallaboratorium für Fernmeldetechnik der AEG.) 


5. Die Entwicklungen im Sonderfall %,;,, @ 2). 


Bei manchen Untersuchungen kommt es darauf an, das allgemeine und auch das nume- 
rische Verhalten der Funktion M,, „,, (2) zu kennen, wenn » und z rein imaginär sind, während 
der Parameter p gleich null oder gleich einer positiven ganzen Zahl ist. Um zu dieser 
Kenntnis zu gelangen, gehen wir von der Integraldarstellung (3.3a) aus. Führt man darin 
für » und z die Größen ir und i{ ein und setzt überdies im Integranden 2x = In (1 +1)/(1—2)), 
so gelangt man nach kurzer Rechnung zu dem Ausdruck: 


\ eos (L/2)- 


Wir wenden hierin auf die cos-Funktion das Additionstheorem an und machen sodann von 
den beiden Gl. (5.2a,b) Gebrauch. Nach einer bekannten Formel von Lommel [4] besteht 


nun aber die für alle A gültige Entwicklung (5.3). Geht man mit ihr in die beiden Teilintegrale 


m: 


von (5.1) ein und vertauscht darnach, was hier zulässig ist, die Reihenfolge von Integration 
und Summation, so stößt man auf Integrale, die sich alle auf ‘das nach [12] auflösbare 


+ 

Bx+r) B-Tin) ® 2 Ren>O für ß,yreell (6.4) 
Fouriersche Integral (5.4) zurückführen lassen. Damit ergibt sich schließlich nach ge- 
eigneter Zusammenfassung der Glieder die in 7 und £ unbeschränkt und absolut konvergente 


ER rtı+2) © riet! 


(T>0,2=0) 


Reihenentwicklung (5.5). Für A=0 ist hierin das unter dem Summenzeichen auftretende 
Produkt einfach gleich 1 zu setzen. Für die erste Ableitung des auf der linken Seite von (5.1) 
stehenden Ausdrucks nach £ und z und für seine zweite Ableitung nach £ und r zugleich, die 
gleichfalls benötigt werden, ergeben sich die ganz ähnlich aufgebauten Beziehungen (5.6), 
(5.7) und (5.8), in denen die innere Summe über r für A=0 gleich Null ist 


*) Die Abschnitte 1 bis 4 des Berichtes sind in dieser Zeitschrift, Bd 23 (1943), S 47 bis 58 erschienen. 
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p+1 
4+1+(pl2) 


1 


(>0,7=0) 


1 


1 ; 


(T>0,2=0) 
i—1 


sgn 27% 


i—1 


; 
.2r 


Aus den obigen Gleichungen geht hervor, daß die Funktion (5.1) für alle reellen Werte 
von r und £ selbst reell ist. Ihrer Integraldarstellung zufolge ändert sie ihren Wert nicht, 
wenn r und £ gleichzeitig ihr Vorzeichen wechseln. Demnach gelten die obigen Entwicklungen 
auch für negative Werte von r, nachdem das wirkliche Vorzeichen von £ in das entgegen- 
gesetzte verwandelt worden ist. Für {>0 gehen die linken Seiten der Gl. (5.5) ... (5.8) 
gegen -die in den Gl. (5.9a, b, c, d) angegebenen Grenzwerte. Geht andererseits 70, so "gelten 


die vier Grenzwertgleichungen (5.10a,b,c,d). Sie ergeben sich aus einer bei N. Nielsen 
angegebenen Formel [13], die in dem Kapitel XIX,'$ 106, dieses Buches zu finden ist. Sie 
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leisten sehr nützliche Dienste sowohl bei der Bestimmung der Lage der Nullstellen von (5.5) 
und (5.6) als zu Kontrollzwecken. 

Für die numerische Berechnung eignen sich die obigen Reihen etwa für Werte von 
t<5 und für Werte von ©&<6...7. Wir benötigen also noch einer Darstellung der obigen 
Funktionen, die sich insbesondere auch für große Werte von r gebrauchen läßt. Um sie zu 
erhalten, schreiben wir das Integral (5.1) als eines mit den Grenzen — und +, das die 
Exponentialfunktion an Stelle der cos-Funktion enthält. Der Faktor vor dem Integralzeichen 
ist dann nur halb so groß. Für ein >0 darf in diesem neuen Integral der Integrationsweg 
längs der reellen Achse ohne Änderung seines Wertes durch den unendlich fernen Halbkreis 
der oberen x-Halbebene, falls er nur nicht gerade durch eine Nullstelle des Nenners geht, zu 
einem geschlossenen Weg ergänzt werden. Innerhalb dieses Weges liegen die unendlich 
vielen, wesentlich singulären Stellen <= (z/2)-i-(24+1) mit A=0,1,2,... Wir denken sie 
uns zunächst durch kleine Kreise aus dem Integrationsbereich ausgeschlossen. Für die Um- 
läufe längs dieser Kreise um die einzelnen singulären Stellen nimmt das Integral bis auf 
einen für diese Zwischenbetrachtung unwesentlichen Faktor die Form an: 


i@r2— - Cotge) 
2ni (Sinz) +1 


Da hierin das Integral von A gar nicht abhängt, so läßt sich die Summation ohne 
weiteres ausführen. Das Umlaufsintegral selbst kann aber mittels der Hilfsfunktion 9 (z) 


i=0 


von Gl. (5.11) auch in der Gestalt des Integrals 


eilzrz— 2.271) 
pri :p(2)-dz 


geschrieben werden. Macht man hierin nach der Substitution 2=(t/2) (£/r)"» von der be- 
kannten Formel (5.12) Gebrauch und löst den ganzen Integranden in eine doppelt unendliche 


ee _ _Ym. „eivrd= Nun. 


Reihe auf, so verschwinden in der Entwicklung nach der Ausführung der Integration alle 
Glieder mit einer anderen Potenz als t'. In Rücksicht auf die durch Gl. (5.13) beschriebene 
Umbenennung der Koeffizienten entsteht auf diese Weise schließlich für die Funktion von 
Gl. (5.1) die neue Darstellung durch Gl. (5.14a). Sie gilt den bisherigen Annahmen gemäß 


M; (id) 


i—0 i=1 


zunächst nur für positive Werte von z und {. Nach der Integraldarstellung (5.1) muß sie 
aber unverändert bestehen bleiben, wenn beide Parameter negative Werte haben, und tat- 
sächlich trifft diese Behauptung auch an der Reihe zu. Treten sie mit den entgegengesetzten 
Vorzeichen auf, so braucht man nur in Gl. (ö.14a) etwa £ durch |{|-e+”i zu ersetzen und die 
bekannte Umlaufsrelation der Besselschen Funktionen anzuwenden. Die ‚dadurch ent- 
stehende Entwicklung lautet: 


(5.6) 
-0) 
| | 
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_ 272:I(p+1) 


| 
(et rer +ir) 


T>0, 
Die in den Gl. (5.14) auftretenden Reihenkoeffizienten lassen sich unabhängig von den 
übrigen Rechnungen aus Gl. (5.14) gesondert bestimmen. Für sie gilt bis zu einem’i<7 


die allgemeine Formel (5.15). Für A=3 z. B. besteht darnach c‘”’ nur aus zwei Gliedern, 
da 1/(—1)! usw. verschwindet. 


35 3.5.7 Tag! 


3481 +89) -5p + 564° — 5164? + 2986-1 — 9111 
* 37.5.7 


5.15) 


Da für A>o sowohl die Funktion J,(2y|r£|) als auch die Funktion I, (2y|r£|) asymp- 
totisch durch den Ausdruck 


angenähert wird und sich dabei nach einer bekannten Abschätzungsformel für die Koeffizienten 
p+1 

die Ungleichung 2 2 berechnet, so konver- 

gieren die Reihen (5.14a,b) absolut und unbeschränkt, und zwar sowohl in bezug auf r als 

auch in bezug auf £. Für größere Werte von £ werden sie allerdings für die numerische 

Berechnung sehr bald ungeeignet. 

Der besondere Aufbau der obigen beiden Reihen gestattet ohne Mühe auch die Her- 
stellung der ersten Ableitung nach {. Wir geben die dafür geltende Beziehung jedoch nur 
für den Fall an, in.dem r und £ gleiche Vorzeichen haben, da ein Blick auf die Gl]. (5.14a, b) 
sofort erkennen läßt, welche Änderungen der entgegengesetzte Fall verlangt. Es ist darnach 


d Mina a-I(p+1) 1 


de + MP ler + (Pete) +in) 


| 

+1 


i=0 


und hierin ist bis zu einem A =6: 
3-5-7 


43102024665) -5p+(562° — 1188 274-3082 — 7035) 
32.5? 


(5.168). 


Auf die Angabe der Formeln, die bei der vorliegenden Entwicklungsart den früheren Gl. (5.8) 
und (5.9) entsprechen, soll an dieser Stelle verzichtet werden. Dagegen mag hier noch der 
Hinweis Platz finden, daß die Kenntnis der Funktior. (5.1) und ihrer ersten Ableitung nach £ 
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“gemäß den Gl. (5.17a,b) auch zugleich über das Verhalten der Wellenfunktion m‘P) (if) und 
deren erster Ableitung nach £ unterrichtet. 


Ist im Falle eines r>0, £<0 der Parameter r sehr groß, so ist nach Gl. (5.14b) in 
erster Nährerung 
Wegen des großen Wertes von r kann man hierin an Stelle der Besselschen Funktion ihre 
asymptotische Näherung verwenden. Wird dann noch wegen {< 0 wie bisher z=|{|- exp(— 
und v=r-exp(+ri/2) gesetzt, so wird |'{|=vz, und es entsteht für M,, 7, (2) selbst die 
folgende asymptotische Darstellung: 


Sie ist in dieser Form, wie man aus Gl. (5.14a) ersieht, auch dann noch gültig, wenn » und 2 
imaginäre Zahlen gleichen Vorzeichens sind. 


Überhaupt ist zu den vorstehenden Entwicklungen zu bemerken, daß ihre Gültigkeit 
keineswegs auf die Annahme rein imaginärer Werte von z und » beschränkt ist. Das folgt 
einmal aus dem Prinzip der analytischen Fortsetzung, da die beiden Reihen (5.5) und (5.14 a) 
ihre absolute und unbeschränkte Konvergenz nicht einbüßen, wenn darin £ durch 2/i und r 
durch v/i ersetzt wird. Es geht dies aber auch aus der Integraldarstellung (5.1) hervor.- 
Schreibt man diese nämlich in der Form des Integrals (5.19) und verfügt darin über den 


2 jarez|<n 


(5.19) 


Verlauf des Integrationsweges, daß er im Endlichen sämtliche Nullstellen von Cojs mit 
Ym(s)>0 zu umschlingen hat, so kann man sich leicht davon überzeugen, daß (5.19) in der 


Tat eine Lösung der Differentialgleichung (5.20) für M,, BR 9 Tu ? darstellt, falls am An- 
(5.20) 

fang und Ende des Weges mit »—=|r|- eiv 
2-|r|-sin(o +(p+ND- |sing|>0 . . . 2.2.2.2... (8.21a) 


ist. Für o=y=n/2 und z=i:.{ geht das Integral (5.19) in das frühere Integral (5.1) über. 
Für die Integralgrenzen © - exp (—ni/2-+iy) als obere und © -exp (— ni/2—iy) als untere 
Grenze gilt von (5.19) das gleiche. Es muß dann nur statt (5.21a) die Bedingung 


erfüllt sein. Der die beiden Enden des Integrationsweges von (5.19) verbindende unend- 
lich ferne Kreisbogen muß natürlich, damit er nur einen verschwindenden Beitrag liefert, 
so geführt werden, daß er die imaginäre Achse der s-Ebene zwischen zwei Polen durch- 
schneidet. Da der Winkel in den beiden Ungleichungen (5.21a, b) beliebig klein gemacht werden 
kann, so darf z.B. in (5.19) o=are (vr) einen Wert haben, der im Bereich ö<o<r—ö mit 
ö>0liegt. Der Übergang von dem Integral (5.19) zu den Reihen (5.5) und (5.14a) läßt sich 
wie früher durchführen. 

Aus der Tatsache der für alle Werte von z2=i{ und v=i-r bestehenden absoluten 
und unbeschränkten Konvergenz der Reihe (5.5) folgt aufs neue, daß die Funktion M,,»,: (2) 
sowohl in 2 als auch in » eine ganze transzendente Funktion ist, die in » zugleich eindeutig 
ist. Das gleiche gilt dann in Rücksicht auf die Gl. (2.5a,b) von der Funktion W,,,ı: (2). 
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a) Der kurvenmäßige 
Verlauf und die Nullstellen 
der Funktionen m!" (i£) und 
Eine anschauliche 
Vorstellung von dem durch (5.5) 
und (5.14) berechenbaren Verlauf 


der Funktion m!"(i£) vermitteln 
das Bild 4 und die Zahlen- 
tafel I. Beide zeigen ihn in 
dem Bereich |r|<3,5 in sogen, 
perspektivischer Darstellung, und 
zwar für ein p=0 und für die 
Parameterwerte {=0,1,2...6. 


Aus dem Bild läßt sich zugleich 
angenähert die Lage der kleinsten 
Nullstellen ersehen. Mit zunehmen- 
den Werten von £ rückt darnach 
die erste Nullstelle von &) 
in bezug auf r bei gegebenem Z schnell von großen negativen Werten gegen Null. Die Stelle 
z=0 selbst bildet nach Gl. (5.104) die kleinste Nullstelle für diejenige Funktion m! (i£), 
die dem Argument {=2j,, = 4,80 %6 entspricht. Für noch größere {-Werte wird die kleinste 
Wurzel in bezug auf r dann sogar positiv. Um eine bequeme Interpolation zu ermöglichen, 
ist in Bild 5 die Abhängigkeit der kleinsten Wurzeln r, von £ graphisch wiedergegeben. Die 
Zahlentafel II bringt die Nullstellen mit einer Genauigkeit bis zu 4 Dezimalstellen. Eine 
formelmäßige Darstellung gerade der kleinsten Wurzel dürfte schwerlich anzugeben sein. Da- 
gegen macht es keine Schwierigkeit, eine Formel für die direkte Berechnung der größeren 


Bild 4. Verlauf der Funktion Y=Mjr, o (id) (i 
in perspektivischer Darstellung. 


Zahlentafel I. Tafel der Funktionswerte von (2/n)"2 m;,‘” (i£). 


1 


2 


3 


| 


4 


5 


6 


— 0,37 563 
— 0,3163 
— 0,20 975 
— 0,04 493 
+0,19 117 
+0,51 333 
+0,93 847 
+ 1,48 581 
+ 2,17715 
+3,03 711 
+ 4,09 346 
+ 5,37 740 
+ 6,92 391 


— 0,16 606 
— 0,28 622 
— 0,36 984 
— 0,37 540 
— 0,24 492 
+.0,10 038 
+ 0,76 520 
+ 1,88 634 
+ 3,64 000 
+ 6,25 044 

+ 10,00 013 

+ 15,24 174 

+ 22,41 207 


+0,20 121 
+ 0,06 538 
— 0,11413 
— 0,28 980 
— 0,36 802 
— 0,18712 
+0,51 183 
+ 2,12 354 
+ 5,22 859 
+ 10,65 797 
+ 19,57 552 
+ 33,58 23 
+ 54,84 81 


+ 0,27 635 
+ 0,27 081 
+ 0,16 083 
+0,06 512 
— 0,27 812 
— 0,32 765 
+0,22 389 
+ 2,14 197 
+ 6,72 328 
+ 16,03 922 
+ 33,38 606 
+ 63,60 201 
113,91 17 


Zahlentafel II. 


+ 0,10 531 
+ 0,23 487 
+ 0,27 646 
+ 0,16918 
— 0,08 766 
— 0,33 052 
— 0,04 384 
+ 1,92 148 
+ 7,87 791 

+ 22,00 278 

+ 51,48 380 

108,16 81 
210.92 53 


— 0,11 433 
+ 0,05 262 
+ 0,21 582 
+ 0,26 647 
+ 0,10 393 
— 0,23 019 
— 0,26 005 
+ 1,48 343 
+ 8,46 297 


 +27,88 341 
+73,21 734 


169,11 05 
357,70 20 


Tabelle der kleinsten Wurzeln 
der Gl. m" 


Tı 


T, 


Bild 5 (links). 


— 
— 1,39 097 
— 0,61138 
— 0,30 959 


— 0,12 190 
+0 

+ 0,02 685 

+0,17 600 


Die $-Abhängigkeit der kleinsten Nullstelle 
der Gleichung Mjr,o dd) in bezug auf r. 


| 
— 25 | | | 
— 2,0 | | | 
—15 | 
| | 
+00 | | | | 
+0,5 | | | | | 
+10 | | | | 
+ 2,5 | 
| | 
06 
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Zahlentafel III. Tafel der Funktionswerte von (i£)/DE. 


T zugleich 


1 


2 


3 


4 


6 


— 2,0 
— 1,0 
—05 
+ 0,0 
+05 
+1,0 
+15 
+2,0 
+25 


+30 


— 0,23 608 
— 0,41.782 
— 0,56 027 
— 0,63 739 
— 0,61 717 
— 0,46 069 
— 0,12 113 
+0,45 727 
+ 1,34 052 
+ 2,60 614 
+4,34 471 
+6,66 146 
+9,67 811 


+0,42 689 
+0,30 946 
+0,13 386 
— 0,07 431 
— 0,26 525 
— 0,35 647 
— 0,22 003 
+ 0,33 049 
+ 1,55 810 
+ 3,82 446 
+ 7,61 989 
+ 13,57 127 
+ 22,52 96 


+ 0,24 085 
+ 0,31 940 
+ 0,31 184 
+ 0,20 166 
+ 0,00 202 
— 0,21 494 
— 0,27 847 
+ 0,13 445 
+ 1,58 179 
+ 4,95 719 
+ 11,63 415 
+ 23,65 25 
+ 43,96 59 


— 0,07 712 
+0,07 915 
+ 0,20 931 
+0,25 251 
+0,15 859 
— 0,06 777 
— 0,28 836 
— 0,10 105 
+ 1,36 585 
+ 5,76 625 

+ 15,99 89 

+ 36,87 62 

+ 76,07 41 


— 0,23 281 
— 0,13 274 
+ 0,02 036 
+0,16 535 
+0,20 545 


+0,05 611° 


— 0,24 855 
— 0,33 634 
+ 0,90 487 
+ 6,01 385 

+ 20,08 99 
+ 52,57 345 
+ 119,99 1 


— 0,18 381 
— 0,20 678 
- — 0,12 828 
+0,04 155 
+ 0,16 601 
+0,18 611 
—0,16953 _ 
— 0,52 939 
+0,23 746 
+5,49 712 
+ 23,13 43 
+ 69,33 19 
+ 175,24 4 


und ganz großen Wurzeln aufzustellen. Beschränkt man sich dabei auf den Fall p=0, so 
errechnet sich dafür aus (5.14b) der folgende Näherungsausdruck : 
2/16 2/32? + 160 — £°11944 
Jon Jon 
Wendet man ihn für die Z-Werte der Tabelle auf die Berechnung der zweiten und dritten 
Nullstellen an, so ergibt sich aus (5.22) 7, bereits mit einer Genauigkeit, die der auf graphi- 


mt .. . (5.2). 


'schem Wege erreichbaren gleichkommt. 


Das Bild 6 bringt eine perspektivische 
Darstellung von dem Verlauf der Funktion 
om;,(i{)/dZ für den gleichen Bereich von r 
und £ wie vordem. Das Bild 7 veranschaulicht 
die Abhängigkeit der ersten drei Nullstellen 
dieser Funktion vom Argument {. Die Zahlen- 
tafeln III und IV ergänzen diese mehr quali- 
tativen Angaben nach der numerischen Seite 
hin. Sie entstammen der. unter [23] erwähnten 
Arbeit des Verfassers. 


' Die numerische Berechnung und die Auf- 
zeichnung der genannten beiden Funktionen hat 
in dankenswerter Weise Herr N. Karsmenkoff 
durchgeführt. 


Bild 6. 
Der Verlauf der Funktion 
> 
in perspektivischer Darstellung. 


6. Eine besondere Entwicklung für W;,,(id). 


Für die Anwendungen ist es außerdem erforderlich, auch für die Funktion W,,7,, (2) bei 
rein imaginären Werten von » und z über eine für die numerische Berechnung geeignete 
Reihenentwicklung zu verfügen. Wir wollen sie an dieser Stelle nur für den wichtigsten 
Sonderfall p=0 herleiten. Setzen wir zu diesem Zweck p zunächst noch als ein wenig von 
Null verschieden voraus und bilden mittels (3.3a) die Gl. (2.5a) nach, so führt der Grenz- 
übergang p>0 unter dem Integralzeichen ohne weiteres zu der Formel: 


| 
| 
| 
| 
| 
| | 
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582 
647 
393 
019 | 
005 
343 
297 
341 
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Zahlentafel IV. 


Tafel der ersten drei Nullstellen 


’ ’ 


| 
T, T; 
| 


— 00 — 
+0,12361 | — 3,5893 | — 12,92 
0,26004 | — 1,6756 | 5,99 
+040715 | — 0,9956 3.84 
+057031 —-0,6310 | 2,75 
 +0,74854 | — 0,3957 2.06 
'+093881 | — 0,2236 | 1,58 


Bild 7 (linke). 1. 2. und 3. Nullstelle von 


in Abhängigkeit von £. 


(6.1). 


Das hierin auftretende Integral läßt sich ganz ebenso behandeln wie das Integral (3.3a) 
zu Beginn des Abschnitts 5. Auch die übrigen Integrale, auf die man bei der weiteren Aus- 
wertung stößt, sind alle wie früher auflösbar, denn da in (5.4) der Parameter n, sofern nur 
Re (n) >0 ist, ganz beliebige Werte haben darf, so kann man ihm auch die Form n=n+e 
geben. Die Differentiation nach e mit dem nachfolgenden Grenzübergang e>0 führt dann 


unmittelbar von Gl. (5.4) zu der noch allgemeineren Formel (6.2). Die weiteren Rechnungen 
reihen sich mühelos aneinander, und für das in (6.1) auftretende Integral entsteht auf diese 
Weise die absolut und unbeschränkt konvergente Reihenentwicklung (6.3), die für die nume- 


+1 


(6.2) 


i=0 


4 


2r 4r 


rische Berechnung unter den gleichen Bedingungen geeignet ist wie die Reihe (5.5). Die 
entsprechende Reihe für die Funktionen W;,,p1), (#£) mit p=1,2,,.. kann nötigenfalls aus 
der Reihe (6.3) und ihrer Ableitung nach £ mit Hilfe der gleich zu erwähnenden Rücklauf- 
formeln hergeleitet werden. 

An allgemeinen Erkenntnissen, die die Gl. (6.1) erschließt, ist noch die Feststellung 


von Wichtigkeit, daß, während sich die Funktion m\"(i£) bei einem gleichzeitigen Vorzeichen- 
wechsel von r und £ überhaupt nicht ändert, die Funktion »‘"(i£) dabei in die konjugiert 
komplexe Größe übergeht. 


7 | 
| 
( 
t 
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7. Die Rücklaufformeln für die Wellenfunktionen m und w. 


Wie für die meisten der höheren -transzendenten Funktionen, die in der mathematischen 
Physik eine Rolle spielen, so bestehen aus bekannten Gründen auch für die beiden konfluenten 
hypergeometrischen Funktionen Beziehungen der Art, die je drei dieser Funktionen in eine 
homogen® und lineare Abhängigkeit voneinander bringen. Jedoch tritt uns hier der Fall ent: 
gegen, daß diese Beziehungen für die m- und w-Funktionen verschieden ausfallen®). Man ge- 
winnt sie für beliebige Zeigerwerte am bequemsten aus den Integraldarstellungen (3.4) oder (4.1), 
indem man darin als Integranden die Ableitungen der wirklichen Integranden nach der je- 
weiligen Integrationsvariablen benutzt. Da diese Integrale matürlich einen verschwindenden 
Wert haben müssen, andererseits die Differentiation aber auch ausgeführt werden kann, so 
liefert dieses Verfahren sofort eine verschwindende dreigliedrige Summe, in der jeder Summand 
hinterher als eine m- oder w-Funktion gedeutet werden kann. Hat man erst einmal einige 
wenige solcher Beziehungen aufgefunden, so lassen sich aus ihnen durch geeignete Kombi- 
nation zahlreiche weitere Formeln gewinnen. Auch die Gl. (2.4) und (2.7) können dazu mit 
Erfolg herangezogen werden. Einige der grundlegenden Beziehungen dieser Art teilen wir im 
folgenden mit: 


2 


1 1» 
2: (2) =(+> _ >) (2) — -5) 


1 
Der Ableitungsstrich bedeutet darin wie immer eine Differentiation nach dem Argument. 


8. Die Asymptotik der M- und W-Funktionen. 


Zur Beurteilung von Konvergenzfragen bei Reihen und Integralen ist es sehr erwünscht, 
das Verhalten der Funktionen M und W zu kennen, wenn entweder das Argument z oder die 
Zeiger » und p in verschiedenen Richtungen ihrer Ebenen ins Unendliche gehen. 


a)Das asymptotische Verhalten für 2>x. Für den Fall, daß es das Argument 
selbst ist, das gegen © strebt, zeigt von den in Rede, stehenden Funktionen die Funktion 
W,,»i,(2) das einfachere Verhalten. Die dafür maßgebende Formel (S.1) läßt sich am be- 
quemsten aus der Integraldarstellung (4.6) dieser Funktion herleiten, und zwar einfach da- 
durch, daß man den Integrationsweg naclı links [3] verschiebt. Eine andere Möglichkeit, die 
asymptotische Entwicklung der Funktion W,,,ı, zu gewinnen, bietet die Gl. (4.1), indem man 
darin v— 1=u/z setzt und den Faktor (1+u/z) mit dem Exponenten (p —1)/2+» in die 
binomische Reihe entwickelt. Die entstehende Reihe hat nach einem bekannten Satz [4] in 
der Tat asymptotischen Charakter. Die auftretenden Integrale lassen sich durchweg durch 
Gammafunktionen ausdrücken. 


6) Diese Unbequemlichkeit ließe sich vermeiden, wenn man anstatt mit den Funktionen m und w mit den ver- 
all-. :seinerten Laguerre-Funktionen arbeiten würde. Das bringt jedoch bei den physikalischen Anwendungen wieder 
den Nachteil mit sich, daß die Schreibweise weit weniger geschlossen ausfällt, da man es hier immer mit den Funk- 
tionen m und w zu tun hat. 
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1 1—p 


Die entsprechende Darstellung für die Funktion M,,»ı, (2) ergibt sich dann ohne weitere 
Rechnung aus der Beziehung (2.9). In der Bezeichnungsweise der Gl. (8.1) lautet sie: 


P(p+1 1 1— 


) 


In Übereinstimmung mit den früheren Festsetzungen ist in (8.2) das +- oder das —-Zeichen 
zu wählen, je nachdem — 37/2 <arc (2) <+n/2 oder —a/2<arc (2) <+3n/2 ist. Die Zeiger 
p und » können in (8.1) und (8.2) jeden beliebigen endlichen Wert haben, der /’(p+ 1) nicht 
unendlich macht. Ist im besonderen v„=n-+(p+ 1)/2 und » eine nichtnegative ganze Zahl, 
so geht in der rechten z-Halbebene für z>x die Funktion M,, »ı2 (2) sr gegen 
Null, weil das erste Glied in (8.2) dann fortfällt. 

b) Das asymptotische Verhalten für »>x. Der andere Grenzfall, bei dem für 
endliche Werte von z und p der Parameter » über alle Grenzen wächst, ist für die Funktion 
M,,»ı» (2) wenigstens unter der Annahme eines reellwertigen v u. a. bereits von OÖ. Perron [10] 
untersucht worden. Er bedient sich dabei der Sattelpunktsmethode in der von ihm verall- 
gemeinerten Form und wendet sie auf die Integraldarstellung (3.2a) an. Den allgemeineren 
Fall eines „>> mit Re» >0 hat nach der gleichen Methode A. Erdelyi [7] an dem 
Integral (3.6) behandelt. Das Ergebnis dieser Untersuchungen läßt sich auf die Form der 

1 
M,„,»ı2(2) r-cos(@yrz— (Rev>0) . . . (83a) 
Gl. (8.3a) bringen. Wie der Vergleich mit (5.18) zeigt, trifft demnach in der Tat die Ver- 
mutung zu, daß der asymptotische Näherungsausdruck (8.3a), der bisher nur unter der An- 
nahme eines Re» >0 bewiesen worden ist, auch noch für rein imaginäre Werte von » gültig 
ist. Aus Gl. (2.4a,b) folgt dann, daß für die Funktion M_,,»1 (2) mit Re» >0 die Dar- 


stellung (8.3b) gelten muß. Dabei dürfte es sich wohl von selbst verstehen, daß, wenn das 
Argument der cos-Funktion in (8.3a,b) eine imaginäre Komponente besitzt, nur diejenige der 
beiden sie zusammensetzenden Exponentialfunktionen beibehalten zu werden braucht, die den 
größeren Absolutbetrag hat. 

Wird p zunächst als nicht ganzzahlig angenommen, so können die entsprechenden 
asymptotischen Ausdrücke für die Funktion W, „,„; und W_,,»ı; mit Hilfe der beiden oben 
angegebenen Formeln aus den Gl. (2.4a,b) hergeleitet werden, Für ein Re» >0 und »>x 
ergeben sich damit die Gl. (8.4a,b). Da in ihnen der Parameter p gar nicht vorkommt, so 

— tr ar — al) . . . . (848), 

gelten sie auch für ganzzahlige Werte von p. Der unendlich kine Punkt der v-Ebene ist 
nach diesen Gleichungen ein wesentlich singulärer Punkt. 

Um diese Verhältnisse in einem einzelnen Falle noch besser zu übersehen, möge an 
Hand der obigen Formeln im besonderen der Fall betrachtet werden, daß z=|z]|- exp(+ zi/2) 
ist, da bei der eingangs erwähnten physikalischen Aufgabe das Argument der Wellenfunktionen 
tatsächlich diesen Aufbau hat, und wir wollen uns vornehmen, die Beziehungen für die 
asymptotischen Absolutbeträge der Funktionen M und W aufzustellen, wenn sich der sehr 
weit entfernte Punkt der »-Ebene mit » = |v|- exp (iy) auf einem Kreis mit dem Radius |v| >» 
bewegt. Die Gl. (8.3a,b) ergeben für diesen Fall die beiden Beziehungen: 
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Es folgt daraus: Die Funktion M,.o.(el -exp (+ i/2)) wächst in allen Richtungen der 
».Ebene für »„>x selbst über alle Grenzen. Sie bleibt nur endlich, falls für ein arc(@)=+n/2 
auf dem Strahl arc (— /2) gegen geht oder für arc (2) auf dem Strahl arc (+ n]2). 
Dieses Resultat ist in voller Übereinstimmung mit den Ergebnissen des Abschnittes 5. 

Ein ganz anderes Verhalten zeigt ‚unter sonst gleichen Voraussetzungen die Funktion 
W,,»ı»(2). Hier führen die beiden Gl. (8.4a,b) zu den folgenden beiden Formeln: 

| 2 

. (8.6), 


|2 


. (8.6b). 

p |v|-cosy-In +7 
In Gl. (8.6a) bezieht sich dabei das eingeklammerte Doppelvorzeichen nicht wie in den anderen 
. auf die Wahl z2=|2|-exp(+i), sondern es entspricht den beiden Möglichkeiten eines 
=0. Da das erste Glied in der geschweiften Klammer der Gl. (8.6a,'b), sei es nun positiv 
alır negativ, für ein |v»|> bereits die Größenordnung der links stehenden Funktion be- 
stimmt, so wächst demnach die Funktion W,,»), in der rechten »-Halbebene mit » über alle 
Grenzen, während sie in der linken Halbebene unter jede endliche Grenze sinkt. Eine Nutz- 
anwendung wird von den letzten beiden Gleichungen bei der Aufstellung der Gl. (10.9) gemacht. 
c) Das asymptotische Verhalten für p>x. Schließlich mögen der Voll- 
ständigkeit halber auch noch die Formeln angeführt werden, aus denen sich das asymptotische 
Verhalten für p>x ersehen läßt. Dabei werden wir uns auf die Funktion M,,»ı, (2) be- 

schränken. Es sind dann immer noch zwei Fälle zu unterscheiden. 

Im ersten Fall ist der vordere Zeiger von p unabhängig. Auf Grund der Integral- 
darstellung (3.2a) findet man dann mit Hilfe der Sattelpunktsmethode leicht heraus, daß unter 
dieser Bedingung das Verhalten von M,,»),(2) für ein are (p)|<x/2 dureh die Gl. (8.7a) be- 


1+» 
2 
schrieben wird. Zu dem gleichen Ergebnis ist auf andere Weise auch A. Erdelyi [10] 


gelangt. 

Im zweiten Falle hänge auch der vordere Zeiger von p ab, indem etwa wie bei den 
Laguerre-Polynomen »=u+(p+1)/2 ist. Die Funktion M,,,ı,(2) zeigt dann das durch 
Gl. (8.7b) wiedergegebene Verhalten. Es folgt dies unmittelbar aus der Erklärungsgleichung 
(2.2), und die Behauptung bleibt gültig im ganzen Winkelraum Jarc (p)|<x. 


9. Die Wronskische Determinante zweier linear unabhängiger Wellenfunktionen. 
Das System der Eigenfunktionen eines Rotationsparaboloides. 


Gemäß ihrer Definition erfüllen die beiden durch die Gl. (1.12a, b) eingeführten Wellen- 
funktionen nach (2.1b) die selbstadjungierte Differentialgleichung (9.1) 


Wir betrachten zwei verschiedene solcher Funktionen, die zu den Parametern »,, p, 
und »,, p, gehören und die wir demgemäß mit y{P» und 4{?? bezeichnen wollen. Eine wechsel- 


seitige Multiplikation der für y|?®) geltenden Differentialgleichung mit y\??’ und umgekehrt 
führt dann mit nachfolgender Subtraktion beider Gleichungen zu der Beziehung (9.2), aus der 


sich eine Reihe wichtiger Folgerungen ziehen lassen. 
Wir prüfen als ersten besonderen Fall die Annahme 


(e). 


| 
(8.2). 
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Dann ist in Gl. 9.2) y„=»,=v und p =p,=p zu setzen, und es muß dann nach dieser 
Gleichung 

2: Im,” (z) - — (z) - Const 
sein. Um hierin den Wert der Konstanten zu bestimmen, betrachten wir den Grenzfall z>x 
an Hand der im Abschnitt 8 mitgeteilten Formeln (8.1) und (8.2). Sie führen nach kurzer 
Rechnung zu der wichtigen Beziehung: 


Mit Hilfe der Rücklaufformeln des Abschnitts 7 lassen sich aus ihr zahlreiche weitere Zu- 
sammenhänge ähnlicher Art auffinden. Als Beispiel bringen wir die Gleichungen: 


PL 


Wird in Gl. (93) »=n+(p+1)/2 mit n als einer nichtnegativ ganzen Zahl, so verschwindet 
die rechte Seite dieser Gleichung, worin sich lediglich auf andere Weise die schon durch 
Gl. (4.3) belegte Tatsache ausdrückt, daß die Funktionen w)(z) und m(z) für den an- 
gegebenen Wert von » bis auf einen von z unabhängigen Faktor übereinstimmen. Für v—=0 
entsteht aus Gl. (9.3) im Hinblick auf die Gl. (4.8a) und (5.10a) eine bekannte Formel aus 
der Theorie der Zylinderfunktionen. 

In einem zweiten Sonderfall von Gl. (9.2) si p, =Pp,=Pp, „=n+(p+1)/2 und 
»,=n'’+(p+1)/2 mit n und n’ als zwei positiven ganzen oder verschwindenden Zahlen. Wird 
dann noch y der Wellenfunktion m gleichgesetzt, so geht nach Abschnitt 8a die Funktion m 
für >» mit are (x)=0 exponentiell gegen Null. Die Funktion + (x) ist daher 
im ganzen Grundgebiet 0... endlich und stetig und besitzt dort in bezug auf x wegen 
ihres polynomialen Faktors n voneinander verschiedene positive Wurzeln. Nach Gl. (9.2) 
ist dann für alle n=#n’: 


0 
Ist aber n=n}’, so setzen wir zunächst n" —=n-+e und lassen hierin e>0 gehen. Das Ver- 
halten von mi), p+n,2(&) für »>& muß dann aber, da n-+e nicht mehr ganzzahlig ist, 
aus der vollständigen Formel (8.2) erschlossen werden. Die Berücksichtigung dieses Um- 
standes führt sofort auf die Beziehung (9.4b). Sie entspricht einer bekannten Formel aus 


r T 1 


0 


der Theorie der Laguerre-Polynome. REN mit der Gl. (9.4a) besagt sie, daß die 
aus diesen Polynomen ableitbaren Wellenfunktionen in dem reellen Grundgebiet 0...o ein 
vollständiges Orthogonalsystem bilden mit den Randbedingungen, für z=0 endlich zu bleiben 
und für 2>% exponentiell zu verschwinden. 

In Verbindung mit der Wellengleichung in den Koordinaten eines Rotationsparaboloides 


ist jedoch von größerer Bedeutung das Orthogonalsystem, das aus den Funktionen m! (2) 
mit rein imaginären Werten von » und z gebildet werden kann. Um zu ihm zu gelangen, 
setzen wir in einem dritten Sonderfall von Gl. (9.2) p, =p,=Pp, v=ir, 2=i{ und identi- 
fizieren wieder die Lösung y von Gl. (9.1) mit der Funktion m. Nach Abschnitt 5 hat dann 
sowohl diese Funktion selbst als auch ihre erste Ableitung nach £ für jeden festen Wert 
von £ unendlich viele einfache Nullstellen in bezug auf r, die sich z. B. für ein &>0 ins 
negativ Unendliche erstrecken. Ihre Einfachheit wird sich sogleich zu erkennen geben. Die 
unendliche Folge dieser Wurzeln möge mit 7, oder r,, bezeichnet werden mit n=1,2,3,.... 
Zuerst beschäftige uns der Fall der Wurzelfolge r„ von Gl. (9.5a): 
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und es mögen r„ und r,- zwei verschiedene unter den unendlich vielen Wurzeln dieser 
Gleichung sein. Nach einer Integration über £ zwischen den Grenzen 0 und £ nimmt dann 
die Gl. (9.2) zunächst die Form an: 


Für ein ©={, von Gl. (9.5a) folgt daraus sofort die Beziehung (9.7a), die den ortho- 
gonalen Charakter der Eigenfunktionen beweist, falls dabei der Annahme gemäß n=#n’ ist. 


Ist aber n=n}’, so setzen wir wiederum n"’—=n-+e und also 7» =1r„+6ö und vollziehen an 


.der Gl. (9.6) den Grenzübergang e,ö>0. Auf diese Weise entsteht dann onne Schwierigkeit 


als Gegenstück zu Gl. (9.7a) die Formel: 
Lo 
Mir, (i£) mir (iL,) 
0 


Die rechte Seite dieser Gleichung kann mit Hilfe der im Abschnitt 5 aufgestellten Formeln 
auch numerisch berechnet werden. 


Wir betrachten darnach den Fall der Gl. (9,5b). Ein ganz ähnlicher Berechnungs- 

vorgang führt hier von Gl. (9.6) geradeswegs zu dem Gleichungspaar (9.8a,b). Das System 
der Wellenfunktionen m®% (i£,) mit den Bedingungsgleichungen (9.5a) oder (9.5b) ist wegen 


Lo 


. (9.7b). 


Lo 


08) 


0 


als das zu dem Rotationsparaboloid £, zugeordnete von orthogonalen 
Eigenfunktionen zu betrachten. Sie liegen über dem Grundgebiet 0...2%k&, und erfüllen an 
seinem rechten Randpunkt die Gl. (9.5a) oder (9.5b) und am linken Randpunkt die Forderung 
des Endlichbleibens. Die unter diesen Bedingungen für sie gültige Eigenwertgleichung nimmt 
für z=i{ und »=-—i:-r gemäß Gl. (9.1) die besondere Form der Gl. (9.9) an. Die Voll- 


y+r-z=0 y=dyldd) ...... 099 


ständigkeit des Orthogonalsystems ist nach dem Gang der Herleitung in Abschnitt 5 evident. 
Ebenso bestätigen die Gl. (9.7b) und (9.8b) durch ihren Aufbau unmittelbar die schon eingangs 
ausgesprochene Behauptung von der Einfachheit der Nullstellen, denn da in diesen Gleichungen 
die linke Seite in Rücksicht auf (5.5) stets von Null verschieden ist, so muß es auch die 
rechte Seite sein. 


10. Die Beziehungen für die einfachsten Wellentypen in den Koordinaten eines 
Rotationsparaboloides. 


Wir betrachten an erster Stelle als die einfachste Form einer Wellenbewegung die 
ebene Welle. Das Bezugssystem der (x, 9, z), dessen z-Achse wie in Bild 1 mit der Ro- 
tationsachse des Paraboloides zusammenfallen möge, sei zu dieser ebenen Welle in der Weise 
orientiert, daß die Wellenbewegung parallel zur «2-Ebene vor sich geht. Das Bezugssystem 
der (x’, y',z’), das mit dem ersten Achsenkreuz die y-Achse gemeinsam habe, und zwar in 
der Art, daß die x’- und z’-Achse in die «z-Ebene fallen, zeige mit seiner 2’- Achse in die 
Richtung des Strahlvektors der Welle. Ist dann y der Winkel zwischen der z- und der 2’ 
Achse, so besteht für die maßgebende skalare Wellenfunktion in leicht erkennbarer Bedeutung 
der Größe o die Beziehung: 


+x- sinn — gikte-cosy (10.1). 
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Auf Grund einer bekannten Formel aus der Theorie der Zylinderfunktionen läßt sich aber 
diesem Ausdruck, wenn ö,, den Kroneckerschen Zahlenfaktor bedeutet, die Form geben: 


Das hierin auftretende Produkt aus einer Exponential- und einer Besselschen Funktion 
gestattet jedoch seinerseits mit 2=°—n und o=2:yE£n die folgende Entwicklung [1] in eine 
Reihe von Produkten zweier Laguerre-Polynome: 


(10.2). 


bp» 
Öpa 


er cos"? (x/2) 


(103). 
( tg? La (#2ink) (F2iEk) 

Setzt man hierin die durch die Gl. (1.12) definierte Wellenfunktion m ein, so entsteht 
in Rücksicht auf Gl. (10.2) für eine ebene Welle in den Koordinaten eines Rotations- 
paraboloides die Reihendarstellung (10.4). An Hand der Formeln des Abschnitts 8 kann man 


n= r= 


2iEk) 


(10.4) 


sich leicht davon überzeugen, daß die darin vorkommende Doppelreihe im Sinne der Konvergenz 
von Doppelreihen absolut konvergiert, solange tg(z/2)<1 ist. Das Auftreten ‘zweier m- 
Funktionen in (10.4) steht in Übereinstimmung mit der physikalischen Forderung, daß die 
Funktion ® sowohl für £=0 als auch für n=0 endlich bleiben muß. Für manche Zwecke 
ist es vorteilhafter, statt der Doppelreihe in (10.4) ein Doppelintegral zu verwenden. Wir 
wollen jedoch auf die Durchführung dieser Umformung hier nicht näher eingehen. 

An zweiter Stelle werde der andere wichtige Fall behandelt, daß das Erregungs- 
zentrum der Welle nicht wie bei der ebenen Welle im Unendlichen liegt, sondern daß es im 
Brennpunkt des Paraboloides angeordnet ist. Es kommt dann darauf an, den Ausdruck 
eikR|R=eik@C+")/(&+n) nach den Wellenfunktionen m und w des Rotationsparaboloides zu 
entwickeln. Die Gl. (10.3) legt es nahe, die Darstellung dieses Ausdrucks in leichter Verall- 


sw +1. . . (105) 


gemeinerung durch die Reihe (10.5) zu versuchen, in der h ein zunächst willkürlicher Faktor 
mit |A]<1 ist. ‘Die Vorzeichen der Funktionsargumente sind in (10.5) entsprechend der 
Gl. (1.13) gewählt worden. 

Wir ersetzen in der obigen Reihe die Funktion w durch die Integraldarstellung (4.1b), 
nachdem zur Sicherung der Konvergenz der Integrationsweg vorher entsprechend geschwenkt 
worden ist. Die nach der Vertauschung der Reihenfolge von Summation und Integration 
unter dem Integralzeichen auftretende Reihe kann mittels der bekannten Formel (10.6) durch 


Zr. (1 (i<i;zbel) . . . (106). 


eine elementare Funktion summiert werden. Für die Summe S$") von Gl. (10.5) entsteht da- 
mit die Darstellung (10.7) durch ein uneigentliches Integral, in dem voraussetzungsgemäß 

(-2iEkpr ?—1P .dt 

immer noch |h|< 1 anzunehmen ist. Für h=1 bleibt zwar das Integral konvergent. In 
Rücksicht auf die Konvergenz der Reihe muß dann aber £>n sein. Überdies läßt sich für 
diesen speziellen Wert von h das Integral (10.7) durch die Hankelsche Funktion darstellen. 
Geht man des weiteren in Gl. (10.5) gemäß der Gl. (4.3) vom Laguerre- Polynom zur Wellen- 
funktion mw über, so kommt für h=1 im ganzegjdje Beziehung zustande: 
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(&>n) 


1,(k(E+ n)) 


Für p=0 wird also im besonderen 


eik6+n) eikR 4i 


n—0 


Ist in den beiden Entwicklungen (10.8) und (10.8a) &<n, so braucht lediglich in den 
beiden zugehörigen Reihen die Stellung der Argumente in den Fuuktionen vertauscht zu 
werden. Die der Reihendarstellung (10.8a) entsprechende Integraldarstellung (10.9), die aus 
der Reihe durch Anwendung des Residuensatzes hergeleitet werden kann, ist dieser insofern 

o+iw 
1 


o—i® 


cosrs 


überlegen, als sie sowohl für &>n als auch für &<n gültig ist. Bei der Herstellung des 
Integrals aus der Reihe leisten die Formeln (8.6a,b) gute Dienste, indem sie den strengen 
Nachweis dafür zu erbringen gestatten, daß der Integrand von (10.9) für &>n auf dem rechten 
und für £&<n auf dem linken unendlich fernen Halbkreis der s-Ebene exponentiell ver- 
schwindet. Mit m“ an Stelle von w‘' wäre das nicht der Fall. 

Mit der Herleitung der bisher noch unbekannten Darstellung (10.9) für die primäre 
Wellenstrahlung eines Dipols in den Koordinaten eines Rotationsparaboloides ist im Verein 
mit den gleichfalls neuen Ergebnissen der Abschnitte 5, 6 und 9 über die Eigenfunktion eines 
solchen Paraboloides der entscheidende Schritt auf dem Wege zur Lösung praktisch wichtiger 
Strahlungsaufgaben getan. 


11. Reihenentwicklungen nach den M-Funktionen und Integraldarstellungen 
für die Produkte von M- und W-Funktionen. 


Um auch erkennen zu lassen, in welcher Richtung sich zur Zeit der weitere Ausbau 
der Theorie in rein mathematischer Hinsicht bewegt, mag zum Schluß noch auf eine Reihe 
von Arbeiten der jüngsten Zeit hingewiesen werden. Alle diese Arbeiten geben deutlich zu 
erkennen, wie stark befruchtend auch auf diesem Gebiet die feste Begründung der Theorie 
der Laplace- Transformierten gewirkt hat. 

Zunächst hat A. Erdelyi [14] mit Hilfe dieser Theorie für die M-Funktion verschiedene 
Arten von Additions- und Multiplikationstheoremen bewiesen. In einer weiteren Arbeit [15] 
werden von ihm systematisch Reihen mit M-Funktionen untersucht, in denen entweder der 
erste oder der zweite oder auch beide Parameter zugleich von Glied zu Glied ansteigen. 
Diese Reihen entsprechen den bekannten Neumannschen Reihen erster Art in der Theorie 
der Zylinderfunktionen. 

Ein wichtiges Beispiel für derartige Entwicklungen läßt sich auf sehr einfachem Wege 
aus der Gl. (3.5) herleiten, wenn: man darin 2 durch 2-4? ersetzt und für die dann rechts 
unter dem Integralzeichen auftretende Funktion J,(A-2u 2”) die von Lommel angegebene 
absolut konvergente Entwicklung (5.3) benutzt. Auf diese Weise entsteht die das Multi- 
plikationstheorem der M-Funktion ausdrückende Beziehung (11.1). 


In zwei weiteren Arbeiten sind von dem gleichen Verfasser, und zwar auch hier auf 

der Grundlage der Transformiertentheorie, Integraldarstellungen für die Produkte zweier M- 
oder W-Funktionen hergeleitet worden. In der ersten dieser Arbeiten [16] beweist er z. B. 
die Beziehung (11.1), in der der Integrationsweg, der von —o kommt und nach dorthin zu- 
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p+1 p+1 


lares|<r 

. larc (2,,2,)|< 


rückläuft, alle Singularitäten des Integranden in den Punkten + (2, + 2,)/2 und + (z, — 2,)/2 
umläuft. In einer zweiten Arbeit [17] wird die entsprechende Beziehung für das Produkt 
zweier W-Funktionen abgeleitet. Sie lautet: 


(+) 


— (21 + 23)/2 


p+1 
2 2 
—$) (2, — 8) 


W,, Di »(2,) Wu, pla (2,) = (2, 2.) e 


larc(2,,2,)|< a, jares|< 
»+u&0, —1,—2,... 
Hier umkreist der Integrationsweg nur den Nullpunkt der s-Ebene. Ist Re(ll—u—»v)>0, so 
läßt sich statt des Umlaufintegrals ein gewöhnliches Linienintegral verwenden, und man erhält 


+ 
u— 


p+1 


(z 2,) 2 


. +1 
W,, p12(21) Wu. = —u—») 


(2,,2,=F0) 
Rell— u—r)>0 


Da in (10.8a) v„=—n— 1/2 und a=+n-+1/2 ist und also in diesem Falle tatsächlich 
1—u—rv>0 ist, so kann die Gl. (11.3a) dazu benutzt werden, auch von dieser Seite her die 
Richtigkeit der Gl. (10.8a) zu bestätigen. Für die angegebenen Werte von » und « ist nämlich 
die hypergeometrische Funktion in (11.3a) darstellbar durch die Kugelfunktion 


22,2, 
2, 
Die in Gl. (10.8a) vorgeschriebene Summation nach n, die wegen der absoluten Konvergenz 


von Reihe und Integral auch unter dem Integralzeichen vorgenommen werden kann, wird 
dann auf Grund der bekannten Formel 


hy". Pu, mit |hl<1 
n—0 
ausführbar und führt damit unmittelbar zu dem früheren Ergebnis. 

Wegen weiterer Arbeiten über Integraldarstellungen von Produkten Whittakerscher 
Funktionen, die zumeist auf anderer Grundlage als die der Transformiertentheorie gewonnen 
worden sind, vergleiche man die unter [18] und [19] genannten Abhandlungen von €. S. Meijer 
und W.N. Bailey. 


12. Besondere M- und W-Funktionen. 


Für besondere Werte der beiden Parameter » und p/2 gehen die beiden konfluenten 
hypergeometrischen Funktionen in eine große Zahl bekannter Funktionen der angewandten 
Mathematik über. Wir bringen zum Schluß eine Zusammenstellung der wichtigsten dieser 
Funktionen. Es sind dies 


1. die unvollständige Gammafunktion 


%e, ? .e:®.W,_ı 


Bat 
(+ 4, 4) 
> 
| 
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Damit sind zugleich erfaßt der Integrallogarithmus, der Integralsinus und der Inte- 
gralcosinus. 


2. die Fehlerfunktion 


3. die Funktion des parabolischen Zylinders 


D, (2)=2? 1.2 2.W, ı (2?1/2), 
-7 


4. die von Schlömilch studierte Funktion 


v 
erw 
0 


5. die Hankelsche Funktion 


1 


6. die Besselsche Funktion 
eriv/2 eriv/2 men 
7. die k-Funktion von Bateman [21] 
W,,ı,(2). 


= 


An Polynomen sind von der konfluenten hypergeometrischen Funktion ableitbar die 
schon im Text erwähnten Laguerre-Polynome und außerdem die Polynome von Hermite 
sowie die Weberschen Polynome des parabolischen Zylinders. 
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KLEINE MITTEILUNGEN 


Zur Ausgleichung des mehrfachen Vor- 
wärtseinschneidens. Die Ausgleichung des 
mehrfachen Vorwärtseinschneidens kann man ent- 
weder nach vermittelnden oder nach bedingten 
Beobachtungen vornehmen. Die Ausgleichung nach 
bedingten Beobachtungen kommt praktisch nur 
dann in Frage, wenn nur eine Bedingungsgleichung 
auftritt, es sich also um ein dreifaches Vorwärts- 
einschneiden handelt). Bei der Ausgleichung nach 
vermittelnden Beobachtungen treten zwei Un- 
bekannte auf, wobei es bei der Wahl der Un- 
bekannten verschiedene Möglichkeiten gibt. Am 
nächstliegendsten, ja fast selbstverständlich und 
auch üblich ist es, als Unbekannte die Koordinaten 
des zu bestimmenden Punkts zu wählen; es hat 
dies die Vorteile, daß man die mittleren Fehler 
der gesuchten Koordinaten ohne besondere Rech- 
nung erhält, und daß die Fehlergleichungen alle 
dieselbe Form haben. Im folgenden soll der Fall 
behandelt werden, daß man als Unbekannte zwei, 
mit den gemessenen Richtungswinkeln zusammen- 
hängende Winkel wählt. Bei der dabei sich er- 
gebenden Lösung sind zwei Fehlergleichungen von 
Haus aus linear, während die anderen linear ge- 
macht werden müssen. Die Lösung würde dem- 
nach praktisch insbesondere dann in Frage kom- 
men, wenn es sich um ein dreifaches Vorwärts- 
einschneiden handelt, so daß man nur eine zu- 
nächst nicht lineare Fehlergleichung hat. Die zu 
lösende Aufgabe lautet folgendermaßen: 

In den n Festpunkten Pı, Pe... P,„ mit den 
gegebenen und fehlerfrei anzunehmenden Koordi- 
naten (2, Yı)» Y2) - - Yn), wurden zur 
Bestimmung der Koordinaten (£, n) eines Neu- 
punkts P die Richtungswinkel 9,, je 
gleich genau gemessen. 


1) Vgl. P. Werkmeister: 
Allg. Vermess.-Nachr. Noch nicht 


5. Jahrgang (1943 
erschienen. 


Dreifaches Vorwärtsein- 
schneiden mit Ausgleichung nach bedingten Beobachtungen. 


Wählt man als Unbekannte x und % die beiden 
Winkel P;,P,P und P,P,;P, und bezeichnet man 
die an %,, 92... anzubringenden Verbesse- 
rungen mit %,, ®...%, so lauten die beiden 
ersten Fehlergleichungen, wenn (P,P,) und (P,P,) 
die Richtungswinkel von P,P, bzw. P, P, sind, so: 


Führt man mit Rücksicht auf die anderen Fehler- 
gleichungen für die Unbekannten x und y Nähe- 
rungswerte z, und %, ein, und setzt man 

z=z +42 und y=y,+4y ....() 
so gehen die beiden ersten Fohlengieichuulgen 
über in 

und 


P)— 2) | 
Yo} 


Die übrigen er sind zunächst von 


der Form 


+7 
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oder mit Berücksichtigung der Unbekannten z 
und % 
siny 
Yi 
sinz 

=F (2, y). 


cos (P,Pı) —y} — 


Führt man die Näherungswerte z, und %, ein, SO 
erhält man 
sin (yo 
sin (u + 42) + 
sin (2, + 42) 


sin 
4Y). 
c08 


Wendet man auf diese Gleichung den Satz von 
Taylor an und bezeichnet man die z, und 4, ent- 


Yı-+tPıP; 


tg = 


sprechenden Werte von $ und n mit &, und 0, 


so findet man bei Vernachlässigung der Glieder 
zweiter und höherer Ordnung 


3), 


wobei 


F (2, %) = 


sin ysint(P,P)+ _,, 
Yıt+PıP, sin Yo) _ 
sinz,cost(P,P)— Yo} _ 
sin (2, + %) 


2%,+P,P, 


Yo) _ 


sin %, sin {(P, P,)— 


Yo} 
sin? (2, + %) 


P,P, 


(P,P,)— Yo} 


(&, — 23)? 


und 
Yo) 
sin z,8in{(P,P)+ 2} 


sin 2, {(P, P,)+ 
— (ne — Yı) Sin? 


— 


Bezeichnet man im Gegensatz zu den „gemesse- 
nen Richtungswinkeln“ 9,, 9» und 9; die dem 
Näherungspunkt P, mit den Koordinaten (£,,%0) 
entsprechenden „genäherten Richtungswinkel“ mit 
Und und beachtet man, daß 


&— 2; —=P;P, 


No —Yi 


ne—Yı= PiP, sin 


=tg 


sin 


’sin (+ Yo) 


1 


sin(+Y) 


so erhält man 


P,P,, 


F (2%0,%) =t8 
IF (zu. Yo) _ p P (£, — Sin (No — Yi) cos 
dx cos? 9,.i8in (zu + Yo) 
und 
_ p P (&—2j) sin (No — Yi) COS 
Pi Po? 8in (2, + Yo) 


Da 


tg = 


COS C08 Pi 


99,i C08 
und 9,.i79; ist, so findet man an Hand der 


Gl. (3) für die zuerst nicht linearen Fehler- 
gleichungen die lineare Form 


- (N, 
wobei 

— 2;) sin — (N — Yi) C08 

sin (zu + 


und 
P,P, (& — 2;) sin 9, — (No — Yi) c08 
sin (2, + 


Die Näherungswerte x, und %, wählt man derart, 
daß in den beiden, von Haus aus linearen Fehler- 
gleichungen (2) die Absolutglieder je gleich null 
sind, womit bei der Aufstellung der Normal- 
gleichungen eine kleine Vereinfachung entsteht. 

Nach Berechnung der Koeffizienten a; und 5; 
der Fehlergleichungen (4) kann man mit diesen 
und den Fehlergleichungen (2) in der üblichen 
Weise die Normalgleichungen aufstellen und damit 
Az und 4, bestimmen. Nach Berechnung von x 
und y an Hand der Gl. (1) erhält man die Koordi- 
naten (£, 7) des Neupunkts P aus 


sin Y 
sin 


Nach Berechnung des mittleren Fehlers « eines 
der gemessenen Richtungswinkel 9,, #2 - - Pn 


aus 
[vv] 


wobei man [vv] in der bekannten Weise zur Probe 
doppelt — einerseits im Zusammenhang mit den 
Normalgleichungen und andererseits unmittelbar 
mit den v-Werten — bestimmt, kann man auch die 
mittleren Fehler «; und «, von € und n berechnen. 
Sind A&, und 4n, die durch eine Veränderung 
von 9ı allein um und 4 die durch eine 
Veränderung von 9, allein um «„, sowie 4£&,„ die 
durch eine Veränderung von 9, allein um « je ge- 
trennt an & und n hervorgerufenen Veränderungen, 
so ist nach dem Fehlerfortpflanzungsgesetz 


und 


n=Yy,+P;P; sin {((Pı P)+ x} 


und 
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Dabei sind die „Teilfehler“ oder „Einzelfehler“ 


die Koordinatenunterschiede zwischen den durch g1;, 
Ps *** 9n bestimmten Koordinaten (£,n7) von P und 


den durch (9, # 4,92 *:* In)» (Pi, 92 Pu) 
u) bestimmten Koordinaten 
(E”, + von n Punkten P’, P”...P”, 
‘Die Bestimmung dieser Koordinaten geschieht in 
derselben Weise wie die von ($, 7), wobei ver- 
schiedene Vereinfachungen bei der Zahlenrech- 
nung — z. B. durch Verwendung der Tafeldiffe- 
renzen — sich ergeben. Wählt man bei der 
Bestimmung der verschiedenen und 
Werten die Näherungswerte 7, und y, jeweils so, 
- daß — wie oben angedeutet — die Absolutglieder 
der Fehlergleichungen (2) je gleich null sind, so 
entsteht dadurch bei den Normalgleichungen inso- 
fern eine Vereinfachung, als alle die neuen Nor- 
malgleichungen von den ursprünglichen sich nur 
in den Absolutgliedern unterscheiden. Bei der 
Aufstellung und der Auflösung der verschiedenen 
Systeme von Normalgleichungen wird weit- 
gehendst vom Rechenschieber Gebrauch gemacht. 
Der Gang der Rechnung soll an dem folgenden 
Zahlenbeispiel gezeigt werden: 


Gegeben sind die Koordinaten 


Punkt Y 
m m 


P, — 20 664,34 + 15 054,58 
P, — 22 013,02 + 12 934,92 
P, — % 874,55 + 11392,02. 


Gemessen wurden in P,. P, und P, je gleich 
genau die drei Richtungswinkel 


9, =253* 00’ 50", 9,=831 32’ 00” 
und 9,= 97° 34’ 17". 
Wählt man zur Bestimmung von (£, n) 
28' 52”,8 
und 
%»=(PR,P)— 59% 57'2, 


so erhält man die drei linearen Fehlergleichungen 


—4y 
= — 1,097 42 + 0,304 


Damit ergeben sich in der üblichen Weise die 
Werte 


578 
sowie 
€ = — 21 150,045 m n=-+ 13 464,448m 


Der mittlere Fehler « eines der drei gleich genau 
gemessenen Richtungswinkel ist «= +4”,1. Mit 


9tu=2353° 0’ 54"1, 32’ 00” 
= 9° 34’ 17" 
erhält man mit 

28 und 59 57”2 


und Jy=—0",8, 
und 2=25* 59 26”4, 


und 


und 


— 1,097,42 + 0,304 Ay + 1” 8 


Az—=+0"9 
15° 29 57",8 
# — — 21 150,027 m 


Ay=— ("2 
y=35* 5 
n'= + 13 464,455 m 

An =0,Tcm. 

9, =233° 00 50", +.=31° 32° 04” 1 

erhält man mit 
— 
— — 1,097 42 + 0,304 Ay + 5”,0 
Ay=—("17 
z=15" 238 5% 52"4 
5” = — 21 150,044 m 7" = + 13 464,469 m 
=0,1cm Ans =2,1cm. 
Y%+u=97° 34' 21”1 
erhält man mit 
38 528 und... 59 
Vs . —4y 
— 1,097 42 40,304 4y + 
Az=+1"0 
z—= 15'238’ 58”,8 
&”" =: — 21 150,068 m 
45, —=23cm 


Ay=— ("2 
v=3°'5% 57”,0 
n""—= + 13 464,430 m 
An, =1,8cm. 


Mit den vorstehenden Werten für (A4&,, 41), 
(4&,, 4n,) und (4&,, 4n,) findet man 
<= +2,9cm und un == + 2,9cm. 
Führt man die Ausgleichung mit den Koordi- 
naten des Neupunkts P als Unbekannte in der 
üblichen Weise durch, so findet man 


= — 21 150,046 + 0,029 m 
und 
+ 15 464,147 + 0,030 m, 


also praktisch dieselben Werte. 
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Eine Vereinfachung bei der numerischen 
Integration gewöhnlicher Differential- 
gleichungen. In dem soeben erschienenen Buche 
„Praxis der Differentialgleichungen“ von H. v. 
Sanden (Berlin 1943) werden auf Seite 18ff."und 
45 ff. numerisch-tabellarische Integrationsverfahren 
für Differentialgleichungen erster und zweiter Ord- 
nung beschrieben. 

Es soll hier gezeigt werden, wie sich die Rech- 
nung noch wesentlich vereinfachen und. abkürzen 
läßt. Die Bezeichnungen sind dabei die gleichen 
wie in dem genannten Buch. 

Bei der Integration einer Gleichung erster 
Ordnung und Systemen von Gleichungen erster 
Ordnung sei die Rechnung bis zum Tabellenwert x; 
fortgeschritten, so daß der Funktionswert %7,, sowie 
/x und die Differenzen 4% bekannt sind. 
Dabei bedeutet » die Ordnung der Ersatzfunktion 
bei der Integration; den sog. Annäherungsgrad. 
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Die erste Näherung für den folgenden Funktions- 
wert, %;+,, erhält man durch die auf Seite 23 
angegebene Formel I: 


th 


Die Beiwerte @,, @,; 4,, a, sind die Integrale 


1 
| 
2). 
0 

(i=1;2;---- (,—1) 


Mit yj,, wird fi}, berechnet, wodurch die Diffe- 


renzen 

Der zweite Näherungswert 

Formel II (Seite 23) berechnet 


bestimmt sind. 
soll nach der 
werden, d. h. 


(3) 
mit 
+D(u+2)- 


d 


Diese Summation (3) kann vermieden 
werden. Man erhält einfacher 


(2) a) 


—h:- 


k+ 

also sofort die an der ersten Näherung anzu- 
bringende Korrektur. Die Differenzen (» + 1) ter 
Ordnung führt man ohnehin der Kontrolle wegen 
in der Rechnung mit. Diese Korrektur läßt sich 
fast immer mit dem Rechenschieber, sehr oft im 
Kopf berechnen, so daß die Zeitersparnis gegenüber 
der Summation erheblich ist. 


Der Beweis für die Richtigkeit der neuen Formel 
ist folgender: 


Aus (1) und (3) folgt 


a) 


— a, —a,- 4 


2 
— 


— dp: 


Nun bestehen in der Bezeichnisweise des Buches 
die Gleichungen 


Setzt man dies in (5) ein, so ergibt sich 


(2) 2 1) 


h 


+») 


ap) 


Die Klammern in dieser Gleichung sind nun alle 
gleich Ap. 


Es ist nämlich nach (2) 


= 
1 
2) (u+p-2)(a—1) . 
— du 
(p—N!p 
0 


Mit einer neuen Integrationsvariabelen z=u—1 
wird daher . (4) 


-@+9-1 


0). 


— Ay „feet: 


Es ist daher auch 


= 
Ap-3 + | 
usw. und 
= 


Setzt man dies in (7) ein, so erhält man 


(2) 
h 


Nach (6) ist die eckige Klammer gleich AR ‚so dab 
in der Tat 

ch) 


— som; 
ap ist. 


Rechnet man z. B. mit p=2, so gilt die Formel 

2 5 h 3 

Oder mit p=4: 

Die letzte der Gleichungen (6) enthält noch eine 
willkommene K der .aus neu ge- 
wonnenen Differenzen 4} A, *!, bei 
deren Bildung erfahrungsgemäß leicht Rechanfehler 
vorkommen. 

Mit 


D+1 


werden alle Differenzen +, kontrolliert. — 
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Z. an 
Bd.: 


Diese Gleichung zeigt, daß man die Differenzen 4x y, 
aus und gar nicht zu bilden braucht, 


sondern unmittelbar 
2) ch 


1 
berechnen kann. Die’ Differenzen A;;, sind erst 
dann zu bilden, wenn die Iteration beendet ist. 
Ist OR ‚ der endgültige Wert, so werden die Diffe- 


renzen aus dem zugehörigen berechnet, 
und durch 
1 2 
kontrolliert. Trotzdem wird man die Differenzen 
Ax+ı bereits mit bilden und 
da) 
rechnen, weil man aus den Differenzen eine Ein- 
sicht in den Verlauf der Lösung und die Brauch- 
barkeit des gewählten Tabellenintervalls A gewinnt. 

Bei der Integration von Systemen gekoppelter 
Gleichungen kommt es zudem vor, daß bei einigen 
Gleichungen keine Iteration mehr nötig ist, so daß 
bereits Y,, ‚ der endgültige Wert bleibt. 

Die weiteren Iterationen werden nach der For- 
mel III des Buches ebenso einfach gewonnen. Der 
darin auftretende Zahlfaktor ist auch a), wie man 
aus den vorstehenden Gleichungen (6) durch Aus- 
rechnen der Differenz leicht nach- 
weisen kann. Die Formel III läßt sich also allge- 
mein schreiben 


Ganz ähnliche Überlegungen gelten für die erste 
Iteration bei Differentialgleichungen zweiter Ord- 
nung. (Vgl. 8. 45 u. 46 des Sandenschen Buches). 


Es wird hier. 
(2) ca) 


— —Bı 4, 


Die Koeffizienten A; und B; sind: 


1 u 

ru 
00 

und 


B= 


-)1-1 
Mit Anwendung der Gl. (6) bekommt man 
(2) 
h? 
— Ir + B, + 
A+B+ 


— Bo) 
— 44 (Ap) 


Wiederum ist: 


N! 


vv 


-/[ 
p! 
und mit (11) wird dann 
4p 


k+ı 


(15), 


wobei man 


1 


1 


k+ı 
als Kontrolle benutzt. 
Wenn man z. B. mit p=2 rechnet, ist 
und entsprechend bei p=4 
135 


(20a) 


(2) ca) 


(20b). 
453 
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Hannover. Karl Stohler. 

Eine vom Vorhaltschießen abgeleitete 
Interpretation einer Näherungslösung 
der gewöhnlichen linearen Differen- 
tialgleichung I. Ordnung. Im Band 22 
(1942), Heft 3, S. 143—152 dieser Zeitschrift habe 
ich eine Näherungslösung für die Differential- 
gleichung 


... (WM 
angegeben. Sie hat die Gestalt 

u=f(z-+u(2)). 
wobei «(z) aus der Funktionalgleichung 


für einen gewissen Teil des Definitionsintervalls 
von A(x) eindeutig zu berechnen ist, wenn A (x) 
stetig und stetig differenzierbar, 
sowie |A (z)| hinreichend klein ist. (Vgl. auch den 
auf S. 148 stehenden Koppelungssatz und den auf 
S. 149 befindlichen Phasenverschiebungssatz 1. Art 
der oben zitierten Arbeit.) 

Diese Näherungslösung läßt sich unter gewissen 
Voraussetzungen über / (z)*und A (x) auf anschau- 
liche Weise interpretieren, wenn man die Theorie 
des Vorhaltschießens zu Hilfe nimmt. Dies sei, 
wie folgt, näher erläutert. 

In Bild 1 stellt B eine Batterie dar, die ein in 
der Ebene BAT längs einer Bahn b sich bewegen- 
des Ziel beschießen soll. Dabei darf im allgemeinen 
ein Geschütz aber nicht auf denjenigen Punkt ge- 
richtet werden, in welchem sich das Ziel im 
Augenblick des Abschusses befindet (Abschuß- 


| 
] u 
(14) 
pu 
bal 
Zie 
sic 
| Tr 
Ge 
zu; 
| | un 
Di 
zw 
wi 
Ni 
al: 
trı 
G 
ab 
Sc 
| 06 
Ap (11). 
D 
T 
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Bild 1. 


punkt), sondern auf einen solchen Punkt der Ziel- 
bahn, der erst nach einer bestimmten Zeit vom 
Ziel eingenommen wird (Treffpunkt). Stellt man 
sich unter A den Abschußpunkt, unter T den 
Treffpunkt vor, so muß der Punkt 7 der Forde- 
rung genügen, daß die Geschoßflugzeit Tr, die das 
Geschoß zum Zurücklegen des Weges BT benötigt, 
zugleich auch die Zeit ist, in welcher das Ziel von 
A nach T gelangt. 

Wir bezeichnen die Uhrzeit des Abschusses mit x 
und setzen 


Die Geschoßflugzeit, die ein Geschoß für den Weg 

zwischen Batterie und Abschußpunkt benötigen 

würde, bezeichnen wir mit 7, und setzen 
Tı=—u(). 

Nun können wir den Abschußpunkt A seinerseits 

als Treffpunkt auffassen, der zum Zeitpunkt x ge- 

troffen wird. Damit dies der Fall ist, müßte das 

Geschoß zum früheren Zeitpunkt T,=u(r)+x 

abgefeuert werden. 

Somit gilt 


Dies aber ist dieselbe Funktionalgleichung (3), auf 
welcher die eingangs genannte Näherungslösung 
beruht. 

Somit können wir den Koeffizienten A (x) der 
Differentialgleichung (1) als Geschoßflugzeit zum 
Treffpunkt, die daraus abgeleitete Funktion — u (x) 
als Geschoßflugzeit zum Abschußpunkt interpre- 
tieren. 

Zur Interpretation der Näherungslösung 

u=/(e+w) 
führen wir irgendeine Ortskoordinate des 
Zieles, beispielweise seinen Seitenwinkel ein, be- 
zeichnen denselben für den Treffpunkt mit o7y, 
für den Abschußpunkt mit o,, und beziehen beide 
Winkel auf den Zeitpunkt z des Abschusses. 
Wir setzen 


or=f(2). 


Beachten wir jetzt wieder, daß auch der Abschuß- 
punkt als ein der Abschußzeit z + « (x) zugeord- 
neter Treffpunkt aufgefaßt werden kann, so gilt 


Demnach können wir die Störfunktion /(z) als 
Ortskoordinate des Treffpunktes, die Näherungs- 
lösung u (x) aber als Ortskoordinate des Abschuß- 
punktes auffassen. 
Zusammenfassend kann man sagen: 


Betrachtet man bei der Dgl. (1) f(x) als Treff- 
punktskoordinate, A(z) als Geschoßflugzeit zum 


oder 


1) Die Kenntnis dieser für die Ge- 
schoßflugzeit verdanke ic'ı Herrn Prof. F. Raab, Karlsruhe. 


Treffpunkt eines bewegten Zieles, so stellt « (x) 
die Geschoßflugzeit zum Abschußpunkt und die 
Näherungslösung u (z) die entsprechende Orts- 
koordinate des Abschußpunktes dar. Beim Vor- 
haltschießen z. B. der Flak sind die Ortskoordi- 
naten des Abschlußpunktes als Funktionen der Uhr- 
zeit gegeben, falls die Ladeverzugszeit und die 
sogenannte tote Zeit gleich Null gesetzt werden 
können. Sieht man auch die Geschoßflugzeit zum 
Treffpunkt als gegeben an, so kann man aus den 
gegebenen Größen näherungsweise die Koordinaten 
zum Treffpunkt berechnen. In Strenge ist die Be- 
rechnung deswegen nicht möglich, weil die Ko- 
ordinaten des Abschußpunktes nur bis zum aktu- 
ellen Zeitpunkt bekannt sind. 

In vielen Fällen wird, wenn u (x) eine Orts- 
koordinate des Abschußpunktes, A (2) die Geschoß- 
flugzeit zum Treffpunkt und x die Uhrzeit be- 
deutet, für die Ortskoordihate vu (x + A) des Treff- 
punktes näherungsweise 


auf Grund der Taylorschen Formel gesetzt. 
Findet das Schießen nach den so berechneten 


Treffpunktskoordinaten statt, so gibt das in 
meiner ersten Arbeit definierte Defizit 


— u(2) 


die sogenannte Sprengpunktsablage für die be- 
treffende Koordinate an. 


Berlin. Hans Bückner. 438 


Ergänzungen zu: Kamke, Differential- 
gleichungen: Lösungsmethoden und Lösungen. 


1. Die Differentialgleichung 


Zu den Lösungen dieser Lagrangeschen Glei- 
chung führen am besten geometrische Überlegun- 
zen. Es werde zunächst der Fall c — 0 betrachtet. 
Dann läßt sich die Gleichung in der Form 
schreiben: 


Die triviale Möglichkeit a — b — 0 werde beiseite 
gelassen und der Fall a<0, < 0 ausgeschlossen, 
weil er zu nicht reellen Lösungen führen würde. 
Ebenso soll, wenn von den Konstanten a und b 
die eine verschwindet, die andere positiv sein. 

Man sieht leicht, daß die Isoklinen der Differen- 
tialgleichung (1) eine doppelte Geradenschar dar- 
stellen. Daher bilden auch die Integralkurven 
eine doppelte Schar; die Differentialgleichung ord- 
net infolge ihres quadratischen Charakters jedem 
Punkt der r-y-Ebene zwei Richtungen y’ zu. Die 
Steigungen der durch denselben Punkt gehenden 


beiden Isoklinen sind — 


Demnach erhält man 


die Integralkurven der Differentialgleichung als 
orthogonale Trajektorien der Isoklinen. 

It a=0, 5>0, so gehen die Isoklinen der 
einen Schar sämtlich durch 2 —=0, y -+YVb, die 


der anderen durch y=—Yb. 
sungskurven der Differentialgleichung 


. (N) 
ergeben sich also die Kreise 


Ist a>0, 5=0, so gehen die Isoklinen durch 


=-+Ya,y=0. Die Lösungskurven der Diffe- 
rentialgleichung 


(a>0) . (1”) 


r 
0% 
(14) 
U | 
— 
te 
ng 
22 
abe 
ial- 
2), 
(3) 
len 
auf 
Art Als Lö- 
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sind die Kreise 
®). 


In dem allgemeinen Fall a$0, b #0 haben die 
beiden Scharen der Isoklinen eine gemeinsame 
Einhüllende, deren Tangenten sie also darstellen. 
Die Evolventen der Einhüllenden sind die ortho- 
gonalen Trajektorien der Tangenten und damit die 
Integralkurven der Differentialgleichung. 


Faßt man (1) als die Gleichung der Isoklinen 
mit dem Scharparameter „’ auf, so findet man 
nach bekannten Methoden als Gleichung der Ein- 
hüllenden 


Danach ergeben sich die Integralkurven als 
Evolventen von Ellipsen oder Hyperbeln. Mit Hilfe 
eines Parameters # schreibt sich die Gleichung der 
Ellipsenevolventen (a>0, 5b >0): 


y=Vbsint — 


mit —=Yasin?t + b cos? t 


die der Hyperbelevolventen (a >0, b<0): 


ei 


y=YV-b -b&int — 


@ 


mit 
— Ya Sin®t_ bCofrt 


Der Fall a<0, 5b >0 liefert ein analoges Er- 
ebnis. Das doppelte Vorzeichen entspricht den 
eiden Hyperbelästen. Für ?=, ergibt sich der 
Ort P, der Evolventenspitze auf der Ellipse bzw. 
Hyperbel. Die Werte £>t, führen zu der einen 
Schar der Evolventen (deren Glieder sich z. B. im 
positiven Sinne um den jeweiligen Anfangspunkt 
P, winden), die Werte 2< t, zu der anderen Schar. 
Mit Hilfe einer Tafel der elliptischen Integrale 
lassen sich die Evolventen punktweise berechnen. 


Für b=-a=—r? ist die Binhüllende der Isoklinen 
gemäß (4) ein Kreis. Pie Lösungsgesamtheit der 
Differentialgleichung wird dann gegeben durch die 
Kreisevolventen 


z=r[eost+ (t —t,) sin 
y=rl[sint — (t — t,) cost] 

Man kann sich die Frage vorlegen, wie die Diffe- 
rentialgleichung aussehen würde, wenn die Ein- 
hüllende geradliniger, zu den Lösungskurven 
orthogonaler Isoklinen die Parabel 


ist. Man findet: 
py’+2yy 


mit den Lösungen (Parabelevolventen) ” 


(ar Sin +c) 


(9). 


Die Differentialgleichung (8) ist für = —2a 
identisch mit einer bei Kamket) angegebenen 
Gleichung, deren Lösung auf dem für Lagrange- 
sche Differentialgleichungen üblichen Wege gefun- 
den wurde und aus (9) zu erhalten ist, wenn man 
t durch den Parameter — pt = 2at ersetzt. 

Es sei nun c$#0. Führt man dann unter Be- 
nutzung von 
(a — b)? 

290 


> 
sin?9=_ c082@ 


c 
—b 
die Drehung 


(10) 


y=!sing+ncosp 


aus, so geht die Differentialgleichung in die 


Form (1) über: 


(11). 


Dabei ist 


(a-+ 5)? 


_a+b[, _ — ab) 


wo man nach Belieben die obaren oder die unteren 
Vorzeichen wählen kann. Man kann nun leicht 
angeben, wie die Lösungen bei c=$0 aussehen 
werden. Ist @ >ab, so haben a und f verschie- 
dene Vorzeichen. Dann entstehen für & und n Lö- 
sungen in der Parameterform (6), d. h. Hyperbel- 
evolventen. Im Falle = ab ist entweder @« — 0 
oder # — 0; setzt man dann die andere Konstante 
als positiv voraus, so erhält man die Kreise (2) 
oder (3) als Integralkurven von (11), andernfalls 
richt reelle Lösungen. Für @<ab sind mit a 
und 5 auch « und B entweder beide positiv oder 
beide negativ. Nur im ersten Falle ergeben sich 
reelle Lösungen, und zwar von der Gestalt (5), 


. also Ellipsenevolventen. Die Lösungen für z und 


wenn man die 
d. h.&£ und 7 


y folgen aus denen für & und n, 
Drehung (10) rückgängig macht, 
in (10) einführt. 


Die hier benutzte Methode ist auf alle Gleichun- 
gen der Form 


anwendbar. Ihre Lösungen sind die orthogonalen 
Trajektorien der geradlinigen Isoklinen, und, falls 
die Isoklinen eine Einhüllende haben, die Evol- 
venten der Einhüllenden. 


Peenemünde, Aer. Inst. W. Heybey. 


1)E „Kamke: I, S. 362, Bspl. 
1.406, Leipzig 19 
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2, Die Riecatische Dgl. 


1 1 


besitzt die Lösung 


2, (2) 
+ Z, (2) 


wo Z, (z) eine Zylinderfunktion vter Ordnung be- 
deutet. Durch den Ansatz 


erhält man die lineare Dgl. 2. Ordnung 
au=0.. 
Sie besitzt als Lösung 


. 


Wählt man für die Zylinderfunktionen Hankel- 
sche Funktionen 2. Art, so geht 4 über in die 
Funktion 7 (z), die in der instationären Tragflügel- 
theorie eine Rolle spielt. Sie wurde von H. G. 
Küßner (Luftf.-Forschg. Bd. 13 (1936), S. 415, 
Gl. (29), (30)) eingeführt, wo auch mit leichter Ab- 
änderung die Dgl. (1) angegeben ist (Vorsicht, 
Druckfehler; falsches Vorzeichen). Verwandte 
Funktionen finden sich bei Theodorsen und Cicala. 
Die Dgl. (1) wurde als Grundlage für Reihenent- 
wieklungen benutzt von L. Schwarz (Luftf.-Forschg. 
Bd. 17 (1940). S. 365, GI. (25) ff.) 


Göttingen. Ludwig Schwarz. 


3. Für die Differentialgleichung 


mr 


findet man 
(Sonderfall: Kamke C 7.7.) 


Göttingen. K. Wieghardt. 458 


BUCHBESPRECHUNGEN 


E. PAWELKA, Übungen aus der Mecha- 
nik. Berlin 1942, Springer-Verlag. 70 S. mit 
59 Abb. Preis 3,60 M. 

Das Büchlein wurde in der Absicht geschrieben, 
bei den Studierenden die Liebe zur Mechanik und 
deren Anwcndung bei der Berufsausübung des 
Ingenieurs zu fördern. Es enthält eine Sammlung 
von 40 Übungsaufgaben, die dem. Verfasser als 
Konstrukteur unterkamen und die der Statik, 
Kinematik, Kinetik, Schwingungslehre und Elasti- 
zitätslehre angehören. Leider sind sowohl die 
Stellung, als auch die Lösung der Aufgaben voll 
von Unklarheiten. Hierzu zwei Beispiele: In der 
1. Übung heißt es: „Ein Differentialgetriebe hat 
3 um eine gemeinsame Achse drehbare Systeme 
oder Zweige“ ... — gibt es ein solches Getriebe? 
Eine Abbildung ist nicht beigefügt. — Und in der 
10. wird bei der Bestimmung des Gravitations- 
feldes bei der Bewegung einer Kugel m um eine 
andere M ein Satz verwendet, der folgende Form 
hat: Absolutkraft — Führungskraft + Relativkraft, 
wobei als Führungskraft jene Kraft bezeichnet 
wird, „welche der Masseneinheit die Beschleuni- 
gung des Koordinatensystems, also die des Schwer- 
punktes S von m erteilt“. Es muß bezweifelt 
werden, ob diese Art, Mechanik zu treiben, der an 
sich begrüßenswerten Absicht des Verfassers ge- 
recht zu werden vermag. 


Karlsruhe i. B. Th. Pöschl. 428. 


Dr. WILHELM MÜLLER, Prof. a. d. Univ. Mün- 
chen, Einführung in die Mechanik des 
Fluges. 2. neubearb. Aufl. 146 S. m. 76 Abb. 
Leipzig 1942, Dr. Max Jänecke Verlagsbuchhand- 
lung. Preis brosch. 3,80 M. 


Das Buch gibt in sechs Abschnitten einen Über- 
blick über die Elemente der Flugmechanik. Es 
wendet sich in erster Linie an Studierende; es ist 
daher im allgemeinen von der Benutzung höherer 
mathematischer Hilfsmittel abgesehen. So wird 
7. B. bei der Behandlung des Auftriebes die Me- 
thode der konformen Abbildung nicht. benutzt. Im 
einleitenden Abschnitt sind die Grundgesetze der 
Aerodynamik wiedergegeben. Die nächsten Ab- 
schnitte behandeln Bewegungszustände und Be- 
wegungen des Flugzeugs. In einfacher Weise wer- 
den die Gleichungen für Steig-, Reise-, Gleit- und 


Segelflug abgeleitet und erläutert, anschließend die 
instationären Vorgänge (Start, Landung, Sturzflug, 
Abfangen) untersucht. Die letzten Kapitel be- 
handeln die statische und dynamische Längs- 
stabilität und die Seitenbewegung des Flugzeugs. 

Die Darstellung ist knapp mit Ausnahme der Be- 
handlung der Probleme, an denen der Verfasser 
selbst forschend gearbeitet hat. Der Leser, der 
sich in die Mechanik des Fluges einführen lassen 
will, kann aus dem Büchlein manche Anregung 
schöpfen. Es wäre zu wünschen, daß im Text 
einige Unklarheiten (z. B. Erklärung der Begriffe 
Potential, Zirkulation usw.) und Ünrichtigkeiten 
ausgemerzt würden. 


Göttingen. J. Stüper. 441 


Dr. phil. BRUNO TIEDEMANN, Über Boden- 
untersuchungen bei Entwurf und 
Ausführung von Ingenieurbauten. 
2. verb. Aufl. 40 S. m. 22 Abb. und 5 Tafeln. Ber- 
lin 1942, Verlag Wilhelm Ernst & Sohn. Preis 
brosch. 2 M. 


Daß das Heft von Tiedemann über Bodenunter- 
suchungen, dessen erste Auflage ich in Bd. 21 
(1941, S. 318) dieser Zeitschrift ausführlicher be- 
sprochen habe, in so kurzer Zeit eine zweite Auf- 
lage gefunden hat, beweist seine Notwendigkeit. 
Da Tiedemann einen ausführlicheren Leitfaden in 
Aussicht gestellt hat, war eine umfangreichere Er- 
eänzung des Heftchens nicht zu erwarten. Wegen 
der Bedeutung für neuzeitliche Erdverdichtungs- 
geräte ist auf die Bedeutung der „Klebegrenze* 
neu hingewiesen worden. Das Heft wird bis zum 
Erscheinen des Leitfadens die früher vorhandene 
Lücke weiter ausfüllen. 


Naumburg (Saale). Dr. Scheidig. 442 


Dr. KARL ULLER, a. o. Prof. a. d. Universität 
Gießen, Idee und Begriff der Welle, 
ein Leitfaden durch die physika- 
lische Grundlagen-Forschung des 
Verfassers. 62 S. m. 1 Abb. Berlin 1942, Dr. 
Georg Lüttke Verlag. Preis brosch. 4,50 M. 


Die Schrift bietet eine Einführung in die bereits 
in einer Reihe früherer Veröffentlichungen dar- 
gestellten Gedanken des Verf. über die Behandlung 
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der Wellenbewegungen. Im ersten Teil unterwirft 
er die’in der theoretischen Physik übliche Darstel- 
lung von Wellenvorgängen einer kritischen Be- 
trachtung und faßt die ihr anhaftenden Mängel zu- 
sammen. Im zweiten Abschnitt begründet er dann 
die von ihm geschaffene „Reine Wellenkinematik“ 
und setzt sich mit den "Wellenbeschreibungen von 
Einstein, de Broglie, Schrödinger und Heisenberg 
auseinander. Im letzten Abschnitt schließlich gibt 
er Hinweise auf die Bedeutung seiner Auffassung 
des Wellenbegriffs für die Deutung der physika- 
lischen Erscheinungen. Im übrigen sei auf die Be- 
sprechung seines früher erschienenen Buches „Die 
Entdeckung des Wellenbegriffs“ in dieser Zeit- 
schrift Bd. 18 (1938), S. 203, hingewiesen. 


Berlin. E. Mosch. 463 


HEINRICH DÖRRIE, Grundriß der Phy- 

sik, mit besonderer Berücksichti- 
gung der Anwendungen. 2. ergänzte Aufl. 
XIT+854 8. m. 464 Fi ig. Breslau 1942, Verlag 
Ferdinand Hirt. Preis geb. 19,50 M. 


Die erste Auflage des Buches wurde in dieser 
Zeitschrift, Bd. 20 (1940), S. 366. gewürdigt. In der 
neuen Auflage hat der Verfasser neben kleinen Ver- 
besserungen und einigen durch die Fortschritte der 
Forschung veranlaßten Änderungen an neuen Stof- 
fen solche behandelt, die in der heutigen Physik 
und Technik erhebliche Bedeutung erlangt haben, 
wie die Erscheinungen des Ultraschalls, die Tem- 
peraturerniedrigung im Bereiche tiefster Tempera- 
turen durch magnetische Abkühlung, das Elek- 
tronenmikroskop, die Frequenzstabilisierung durch 
Schwingquarz u. a. m. Das Buch gewährt damit 
einen recht guten Einblick in den gegenwärtigen 
Stand der physikalischen Forschung. 


Berlin. E. Mosch. 462 


Dr. FRITZ SAUTER, Prof. a. d. Techn. Hoch- 
schule München, Differentialgleichun- 
gender Physik. (Sammlung Göschen Bd. 1070). 
147 S. m. 16 Fig. Berlin 1942, Verlag Walter de 
Gruyter & Co. Preis geb. 1,62 M. 

Das Bändchen will eine Brücke schlagen von 
der mathematischen Theorie der Differentialglei- 
chungen zu den praktischen Anwendungen in der 
Physik. Ausgehend von konkreten physikalischen 
Problemen werden an der Hand von gewöhnlichen 
Differentialgleichungen der Mechanik und von 
partiellen Gleichungen der Akustik, Optik, Wellen- 
mechanik, der Potentialtheorie und der Wärme- 
leitung — also durchgehend linearen Gleichungen 
— eine Reihe von Lösungsmethoden entwickelt. 
Insbesondere wird am Schluß kurz auf die Me- 
thoden der Störungstheorie eingegangen. 

Dem Verf. gelingt es so ausgezeichnet, dem 
Leser die verschiedenen Integrationsmethoden im 
Prinzip verständlich zu machen und ihm an den 
zahlreichen Beispielen ein gewisses Maß an 
Rechentechnik zu vermitteln. Dabei setzt er nur 
geringe Vorkenntnisse voraus, verzichtet z. B. auf 
die Anwendung funktionentheoretischer Methoden. 
stellt allerdings an die Mitarbeit des Lesers bei 
der starken Konzentration der Darstellung recht 
beträchtliche Anforderungen. 
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Prof. Dr. BERNHARD BAULE, Die Mathema- 
tik des Naturforschers und Ingeni- 
eurs, Bd. I Differential- und Inte- 
gralreehnung. VII+158 S. m. 161 Abb. 
Leipzig 1942, Verlag S. Hirzel. Preis kart. 7 M. 

Das Buch will dem Naturforscher und Ingenieur 
das Hilfsmittel der Differential- und Integralrech- 
nung in die Hand geben. Diesen Zweck erfüllt es 
in ausgezeichneter Weise. Der Stoff ist klar und 


übersichtlich angeordnet, die Darstellung schwung- 
voll und leicht verständlich. Immer wieder stößt 
man beim Lesen auf Beweise für das didaktische 
Geschick des Autors. 

Der Inhalt gliedert sich so, daß nach einer aus- 
führlichen Einführung des Funktionsbegriffes zu- 
nächst die Differentialrechnung einer Funktion 
einer Veränderlichen entwickelt wird, dann die 
einer Funktion von mehreren Veränderlichen: daran 
schließt sich die Integralrechnung. Begreiflicher- 
weise sind Existenzbeweise häufig weegelassen 
bzw. durch Hinweis auf den geometrischen Sach- 
verhalt ersetzt. Im übrigen vereinigt das Buch 
mathematische Strenge und Verständlichkeit in 
glücklichster Weise. (Mir ist nur eine Stelle auf- 
gefallen, wo man eine präzise Formulierung ver- 
mißt, beim Funktionsbegriff. Der in Abbildung 5 
gezeichnete Kurvenzug enthält Strecken parallel 
zur Y-Achse, sollte also nicht als Bild einer Funk- 
tion Y—=f(x) genannt werden.) Sehr erfrischend 
wirkt es, daß der Autor nicht vor Improvisationen 
zurückschreckt, z. B. in der Verwendung von 
Vektor und Determinante. (In der Bezeichnungs- 
weise geht das vielleicht gelegentlich zu weit, 
z. B. wenn die Gesamtheit der natürlichen Zahlen 
als Zahlkörper bezeichnet wird oder in der Ver- 
wendung des Wortes Matrix.) 

Nützlich wären gelegentliche Literaturhinweise, 
insbesondere auch auf numerische Methoden, die 
vielleicht etwas zu knapp dargestellt worden sind. 

Das Erscheinen des Buches, gerade in diesem 
Zeitpunkt, wird auch von den Studenten mit größ- 
ter Freude begrüßt werden. 
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Prof. Dr. BERNHARD BAULE, Die Mathe- 
matik des Naturforschers und In- 
genieurs, Bd. Il, Ausgleichs- und Nä- 
herungsrechnung. 54 S. m. 30 Abb. Leip- 
zig 1943, Verlag S. Hirzel. Preis kart. 2,80 M. 


Das zweite Bändchen des Werkes ist in der 
Hauptsache Fragen der angewandten Mathematik 
gewidmet, und zwar gibt es eine Einführung in die 
Ausgleichsrechnung nach der Methode der klein- 
sten Quadrate, behandelt die Gaußsche Häufigkeits- 
verteilung. die Darstellung von Kurven durch ein- 
fachere. Fouriersche Reihen und Integrale, Le- 
zendresche Polynome als weiteres Beispiel für ein 
System orthogonaler Funktionen und schließlich 
kurz die Newtonsche und die Lagrangesche Inter- 
polationsformel. 

Was vom ersten Teil gesagt wurde, gilt auch 
vom zweiten. Man freut sich über die klare, über- 
sichtliche und leicht verständliche Darstellung und 
wird das geschickt geschriebene Heftchen gern 
dem Studenten zur Anschaffung empfehlen. Einige 
Versehen (z. B. die Gleichsetzung von glatt und 
monoton) lassen sich leicht bei einer sicher bald 
nötig werdenden Neuauflage beseitigen. Wün- 
schenswert wäre es, wenn dann auch mehr weiter- 
führendes Schrifttum angegeben würde. 
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Dr.-Ing. FRIEDRICH KLINGER, Einführung 
in die n-dimensionale algebraische 
Geometrie, mit besonderer Berück- 
sichtigung der Physik für Tech- 
nikerund Physiker. Herausgeg. mit Unter- 
stützung der Akademie der Wissenschaften in 
Wien. VII + 152 S. m. 1 Tafel. Wien 1942, Sprin- 
ger-Verlag. Preis brosch. 12 M. 


Der Verfasser will lt. Vorwort mit diesem Buch 
die Kluft überbrücken, die nach seiner Ansicht 
zwischen Technikern und Physikern besteht, in- 
dem die Techniker meistens flinke Rechner sind. 
aber selten den hohen Ansprüchen der theoreti- 


| 
sch 
| Ph 
ist. 
me 
tb 
ma 
Au 
ve! 
Ke 
Mi 
] 
un 
ru 
mi 
die 
Ka 
ist 
ve 
ma 
nil 
Ve 
au 
Hs 
dis 
fäı 
ab 
sc 
na 
Pt 
im 
tri 
Ke 
en 
wi 
de 
| wi 
tis 
al: 
Pt 
| ni 
Bı 
eii 
al 
ne 
ha 
ur 
ge 
T: 
Le 
ge 
lu 
di 
ni 
Fi 
W 
gi 
49 
(F 
is 
V 
sc 


h. Mech. 
pril 1943 


hwung- 
r stößt 
ktische 


er aus- 
es zu- 
ınktion 
ın die 
; daran 
flicher- 
elassen 
Sach- 
Buch 
eit in 
le auf- 
ge ver- 
lune 5 
Jarallel 
Funk- 
schend 
tionen 
r von 
nungsS- 
weit, 
Zahlen 
r Ver- 


|weise, 
n, die 
sind. 
liesem 


w. Math. Mech. 
Nr. 2 April 1943 


Buchbesprechungen 127 


schen Physik gewachsen sind, während bei den 
Physikern fast immer das Umgekehrte der Fall 
ist. Dieser Zweck soll durch eine kurze Zusam- 
menstellung von Definitionen und Lehrsätzen 
über Grundbegriffe aus den Gebieten der Mathe- 
matik, Geometrie und Physik erreicht werden. 
Auf die Ableitungen und Beweise wird vielfach 
verzichtet. Vorausgesetzt werden nur geringe 
Kenntnisse in der Mathematik, wie sie in jeder 
Mittel- und Gewerbeschule gelehrt werden. 

Das erste Kapitel (18 S.) behandelt die Zahlen 
und die verschiedenen Rechnungsarten. Kapitel II 
(27 8.) mit der Überschrift „Die Skalare. Einfüh- 
rung in die Physik“ befaßt sich im wesentlichen 
mit den Maßsystemen. Kapitel III (18 S.) bringt 
die Definitionen und wichtigsten Formeln aus der 
Kombinationslehre. Das Hauptstück des Buches, 
Kapitel IV (69 S.), enthält die Begriffe, Regeln und 
Sätze der Graßmannschen Ausdehnungslehre. Es 
ist sicher sehr verdienstlich, daß der Verfasser 
versucht. diese außerhalb des Kreises der Fach- 
mathematiker fast unbekannte Lehre dem Tech- 
niker und Physiker näherzubringen. Wem die 
Vektorrechnung schon geläufig ist, dem wird es 
auch keine besondere Schwierigkeit machen, an 
Hand des Buches die Graßmannsche Lehre im n- 
dimensionalen Raum zu verstehen. Für einen An- 
fänger ohne solche Vorkenntnisse ist sie reichlich 
abstrakt. Das Schlußkapitel (10 S.) hat die Über- 
schrift „Die extensiven Größen des n-dimensio- 
nalen Raumes. Einführung in die geometrische 
Physik“. Diese Größen unterscheiden sich von den 
im vorhergehenden Kapitel betrachtet®n geome- 
trischen Größen dadurch, daß sie als weiteres 
Kennzeichen eine physikalische Dimension besitzen. 

Wie aus der vorstehenden Übersicht hervorgeht, 
enthält das Buch eine etwas unorganische Aus- 
wahl von Begriffen und Lehrsätzen aus verschie- 
denen Gebieten. Wenn der Verfasser im Schluß- 
wort meint, „damit wären wir mit dem mathema- 
tischen Grundbau des gesamten Lehrgebäudes der 
algebraischen Geometrie bzw. der geometrischen 
Physik fertig“, so wird man dieser Auffassung 
nicht ganz beipflichten. Immerhin darf man dem 
Buch nachrühmen, daß es auf beschränktem Raum 
ein Fülle nützlicher Kenntnisse vermittelt. 
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RÜHLMANN-SCHMIEDEL, V ierstelligeLog- 
arithmen- und Zahlentafel. 17. völlig 
neu bearb. Aufl. 104 S. Leipzig 1943, Julius Klink- 
hardt Verlagsbuchhandlung. Preis kart. 2,20 M. 


Für den Schulgebrauch bestimmte Logarithmen- 
und Zahlentafeln mit umfangreichen astronomischen, 
geographischen, physikalischen und chemischen 
Tabellen. Unter anderem werden vierstellige 
Logarithmen vierstelliger Zahlen gegeben. An- 
gehängt ist eine mathematische Formelsamm- 
lung, die der verbessernden Durchsicht be- 
dürfte: z. B. sollten die trigonometrischen Funk- 
tionen nicht nur für den ersten Quadränten defi- 
niert werden; die Kurve y= m-e“Y', bezeichnet 
man nicht als logarithmische Kurve: der Satz „Die 
Funktion y=f(x) hat für den Wert r=a einen 
Wendepunkt, wenn f”(a)=0” bedürfte einer- er- 
gänzenden Aussage über den Wert von f”” (a), usw. 
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OSKAR BOLZA, Aus meinem Leben. 
49 8. München 1936, Verlag Ernst Reinhardt. 
Preis brosch. 1.50 M. 


Oskar Bolza, geb. 12. Mai 1857 in Bergzabern 
(Rheinpfalz), gest. 5. Juli 1942 in Freiburg i. Br., 
ist allen Mathematikern durch sein Lehrbuch über 
Variationsrechnung bekannt. Von seiner Hand ge- 
schrieben liegt hier eine Skizze seines Lebens vor, 


im wesentlichen unter Beschränkung auf die Dar- 
stellung des mathematischen Entwicklungsganges. 

Bolza hat in Straßburg, Berlin und Göttingen 
studiert. Von entscheidendem Einfluß war für ihn 
das Sommersemester 1879, in dem er bei Weier- 
straß Variationsrechnung hörte. Promoviert hat 
er 1886 in Göttingen bei Klein. Zwei Jahre später 
fuhr er nach Amerika und hat dort bis 1910 ge- 
wirkt, von 1892 ab in Chicago zusammen mit 
H. Maschke und E. H. Moore. Bolza kehrte 1910 
in seine Heimat zurück und las in Freiburg als 
Honorarprofessor bis 1932. 

Die Schrift ist mit liebenswürdiger Bescheiden- 
heit verfaßt. Besonderes Interesse dürfte die 
Schilderung seines Amerikaauferffhaltes haben; sie 
fällt in eine Zeit, die für die Pflege der Mathe- 
matik in den Vereinigten Staaten von großer 
Wichtigkeit war. — Bolza streift auch sein In- 
teresse für religiöse und mystische Fragen, das 
sich in zwei Büchern ausgedrückt hat: „Glaubens- 
lose Religion“ und „Meister Eckehart als Mystiker. 
Eine religionspsychologische Studie.“ 
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Ferner sind bei der Schriftleitung folgende 
Bücher eingegangen (ausführliche Besprechung 
bleibt vorbehalten): 


Dr.-Ing. P. WERKMEISTER, o. Prof. em. a. d. 
Techn. Hochschule Dresden, Vermessungs- 
kunde II, Messung von Horizontal- 
winkeln, Festlegung von Punkten im 
Koordinatensystem, Absteckungen. 
(Sammlung Göschen Bd. 469.) 151 S. m. 98 Fig. 
Berlin 1942, Verlag Walter de Gruyter & Co. Preis 
geb. 1,62 M. 


Dr. phil. Dr. techn. FRITZ CHMELKA, wiss. 
Assistent a. d. Techn. Hochschule Wien, und Dipl.- 
Ing. Dr. techn. ERNST MELAN, o. Prof. a. d. 
Techn. Hochschule Wien, Einführung in die 
Statik. VII+ 132 S. m. 119 Textabb. Wien 1942, 
Springer-Verlag. Preis brosch. 6,60 M. 


Dr. GERHARD KOWALEWSKI, Prof. a d. Deut- 
schen Karls-Universität Prag, Einführung in 
die Determinantentheorie einschließ- 
lich der Fredholmschen Determinanten. 
3., verb. u. erweit. Aufl. 320 S. Berlin 1942, Ver- 
lag Walter de Gruyter & Co. Preis geb. 15 M. 


GUSTAV HÖSSJER, Über die konforme 
Abbildung eines veränderlichen Be- 
reiches (Chalmers Tekniska Högskolas Hand- 
lingar, Transactions of Chalmers University of 
Technology, Gothenburg, Sweden, Nr. 10/1942.) 
15 S. Göteborg 1942, N. J. Gumperts Bokhandel 
A.-B. Preis 2 


Die Tätigkeit der Baltischen Geo- 
dätischen Kommission iin den Jahren 
1938—1941. Berichtet von dem Präsidium. 113 S. 
Helsinki 1942, Ösakeyhtio Weilin & Göös, Aktie- 
bolag. 


Dr.-Ing. Dr. techn. E. RAUSCH, a. pl. Prof. a. d. 
Techn. Hochschule Berlin, Maschinenfun- 
damente und andere dynamische 
Bauaufgaben, 3. Teil Rahmenfunda- 
mente beihoher Maschinendrehzahl 
(insbesondere Dampfturbinenfunda- 
mente), Aufstellung von Maschinen 
in Gebäuden, Werkzeugmaschinen, 
Schäden an Maschinenfundamenten 
undErschütterungsschädenanderer 
Art, andere dynamische Aufgaben 
im Bauwesen nebst Ergänzung zum 1. und 
2, Teil des Buches. 365 S. m. 377 Bildern. Berlin 
1942, VDI-Verlag. Preis brosch. 35 M. 
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NACHRICHTEN 
David Hilbert 7. 


Mit Hilbert verlor die Welt den größten Mathematiker, den sie seit Gauß gehabt hat. Hilbert 
war groß durch den Umfang und die Vielseitigkeit seiner Entdeckungen, die von der Algebra und der 
Zahlentheorie über Geometrie und Analysis bis zu den Anwendungen hin einerseits und bis zu den tiefsten 
Forschungen über die Grundlagen seiner Wissenschaft andererseits überall wesentlich Neues gegeben 
haben, so noch einmal die Mathematik als ein Ganzes erscheinen lassend; Hilbert war groß durch die 
tiefschürfende Kraft und revolutionierende Wirkung seiner Forschungsmethoden; Hilbert war groß 
durch die unanfechtbare Strenge seiner Beweisführung, die klassisch geworden ist. 

Hilbert wurde am 23. Januar 1862 zu Königsberg aus einer Juristenfamilie geboren. In Königs- 
berg studierte er auch, war Schüler von H. Weber und Lindemann, promovierte und habilitierte 
sich, wurde außerordentlicher und bald auch ordentlicher Professor. Früh erwarb er Ruhm durch den 
Beweis der Endlichkeit des Invariantensystems und andere algebraische Arbeiten. Der weitblickende 
und großzügige Felix Klein holte den 33-jährigen nach Göttingen, wo nun Hilbert bis zu seinem 
Tode (am 14. Februar 1943) eine erstaunliche Tätigkeit entfaltete, die ihn nicht etwa nur als stillen für 
sich schaffenden Gelehrten zeigte, sondern auch als den großen Anreger, der unbekümmert schenkte 
und eine große Schülerzahl aus der ganzen Welt um sich sammelte. Durch die „Probleme der Mathe- 
matik“, die er auf dem internationalen Mathematikerkongreß in Paris 1900 vortrug, wurde er Lehr- 
meister der ganzen mathematischen Welt. 

Sein zweites großes Arbeitsgebiet wurde die Zahlentheorie. Der „Bericht über die Theorie der 
algebraischen Zahlkörper“ vom Jahre 1897 stellt sich Gauß’ Disquisitiones arithmeticae würdig zur 
Seite. Andere Arbeiten über den Gegenstand folgten; die vollständige Lösung des Waringschen 
Problems ist eine weitere schöne Leistung. 

1899 begannen dann nach einer vorhergehenden kleineren aber sehr hübschen Arbeit über die 
gerade Linie als kürzeste Verbindung zweier Punkte die Untersuchungen über die „Grundlagen der 
Geometrie“. Sie erregten mit Recht das größte Aufsehen. Denn hier wurde mit aller Schärfe wohl zum 
ersten Male die axiomatische Methode in der Geometrie und damit in der Wissenschaft überhaupt ganz 
klar herausgestellt: für den mathematischen Aufbau einer Wissenschaft kommt es nur auf die logischen 
Verknüpfungen der vorangestellten Axiome an, nicht auf deren reale oder anschauliche Bedeutung. 
Kleine Geister haben daraufhin Hilbert dahin mißverstanden, als vertrete er einen extrem nomina- 
listischen Standpunkt: man könne irgendwelche Spielregeln vereinbaren, ganz gleichgültig welche, 
müsse nur streng logisch weiterschließen und die Widerspruchslosigkeit beweisen. Hilbert hat es 
besser gewußt; seine axiomatische Methode war auch für ihn nur eine Seite der Sache; das Kantazitat: 
„so fängt denn alle menschliche Erkentnis mit Anschauungen an, geht von da zu Begriffen und endigt 
mit Ideen“, das er seinen „Grundlagen der Geometrie“ voranstellte, und der Geometer Hilbert, der 
die prachtvolle Vorlesung über anschauliche Geometrie gehalten hat, hätte die Dogmatiker des Konven- 
tionalismus eines Besseren belehren können. 

Kurz nach der Jahrhundertwende beginnen nun Hilberts große Untersuchungen zur Analysis. 
1901 bringt seine Arbeit über das „Dirichletsche Prinzip“, in der es ihm gelingt, mit einer unver- 
gleichlichen Kraft die Schwierigkeiten, die sich einer strengen Begründung dieses an sich so schönen 
und einfachen Prinzips entgegenstellen, zu überwinden. 1904 bis 1910 kommen dann die ebenso 
bedeutungsvollen Untersuchungen zur Theorie der Integralgleichungen, in denen Hilbert das Werkzeug 
der unendlich vielen Variabeln und damit den Hilbertschen Raum schafft. Unvergleichlich ist, wie 
er hier mit aller Strenge Probleme meistert, au die sich wohl sonst niemand herangewagt hätte. Er 
hat damit methodisch der Mathematik die schärfsten Hilfsmittel geliefert. 

Gerade auch die angewandte Mathematik ist hier Hilbert zu größtem Dank verpflichtet, nicht 
etwa nur wegen einiger Probleme, die er gelöst hat, sondern wegen der weittragenden Methoden. Die 
sogenannten direkten Methoden. der Variationsrechnung nehmen von Hilbert ihren Ausgängspunkt 
und bilden heute ein Hilfsmittel, das die Praxis nicht mehr entbehren kann. 

Die Untersuchungen über die Grundlagen der Geometrie mußten einen Geist von der umfassenden 
Kraft Hilberts notwendig weiterführen. Die Widerspruchslosigkeit der Euklidischen und anderer 
Geometrien konnte durch Arithmetisierung nachgewiesen werden. Aber sind die Sätze der Arithmetik 
selbst widerspruchslos? Diese Frage mußte sich erheben, zumal die sogenannten Paradoxien der 
Mengenlehre darauf aufmerksam gemacht hatten, daß man nicht ganz willkürlich ins Blaue hinein 
Begriffssysteme aufbauen kann. Ist die Arithmetik ein unbedingt zuverlässiges Begriffssystem ? Diese 
Frage führte Hilbert zu der kühnen Idee, die Arithmetik selbst als widerspruchsfrei dadurch 
nachzuweisen, daß sie formalisiert wird, d.h. daß eine Kerr: geschaffen wird, deren Symbole 


Begriffe, auch so bedenkliche. wie den Begriff Unendlich, durch eine ndliche überschaubare Anzahl von 
Zeichen mit endlicher Verknüpfungsmöglichkeit darstellen, so daß’die Widerspruchslosigkeit direkt zu 
erkennen wäre. Auch der Widerspruch muß eine bestimmte, Formel erhalten, und es bleibt zu zeigen, 
daß diese Widerspruchsformel zwischen den Zeichen nie auftreten kann. Die Durchführung dieses 
gewaltigen Gedankens hat uns Hilbert, trotz weitgehender Ergebnisse seiner selbst und seiner Schüler 
unvollendet hinterlassen müssen. Doch liegt vielleicht gerade in diesem fast Übermenschlichen, das er 
gewollt und mit Erfolg angegriffen hat, das stärkste Zeugnis seiner Größe. Hamel. 477 
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